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SYNOPSIS 


The subject matter of this dissertation is a study of the 
dynamics and equilibrium properties of a spin 3/2 pure nuclear 
quadrupole resonance subjected to a scheme of 

radiofrequency Crf^ radiation pulse p>ertur bat ions. In magnetic 
resonance literature, the various schemes of rf radiation pulses 
are known collectively as multiple— pulse sequences. 

A central theme that has guided our study is the behavior 
of the system’s response right from zero-time to times of the 
order of where is the spin-CquadrupoleI> lattice 

relaxation time of the sample under study. The establishment 
of quasi -equil ibri um state under varying conditions of the 
external parameters has been studied theoretically. Furthermore, 
wherever possible, appropriate experiments have been p>er formed 
on model quadrupolar spin 3/2 solids to seek support for our 
theoretical findings. Very interestingly, it turns out that the 
mechanism of this quasi -equilibrium state leads to a 
characteristic prolongation of the quadrupolar relaxation time 
CT^Z>, similar to previous findings on the prolongation of the 
spin-spin relaxation time in the case of pure dipolar 

solids. In addition to the abo Yemen tioned focal theme, our 
studies have enabled us to adapt several multiple-pulse cycles 
which have been in use for studies of dipolar solids; the effects 
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of t-hese cycles on the various major interactions in quadrupolar 
solids* namely* electric field gradient Cefg> i nhomogenei ty» 
torsional oscillations and dipole couplings* are studied 

theor eti cal 1 y. 

The thesis consists of seven chapters. In Chapter I we have 
surve]/ed the essentials of the phenomenon of NQR, This has been 
presented vis—' -vis the phenomenon of nuclear magnetic resonance 
CNMR>* and similarities and differences have been pointed out. 
This Chapter also discusses the importance and usefulness of 
multiple-pulse studies in N<3R. A review of the work already done 
in multiple-pulse N<3R is also presented in this Chapter* with a 
preview of the present work outlined towards the end of the 
Chapter . 

In Chapter II we have presented the quantum mechanical 
background for the work presented in later chapters. After 

introducing all the Hand 1 toni ans needed in the usual Cartesian 
basis of spin angular momentum operators* we present these 
Hamiltonians in terms of the single transition operators CSTO> in 
a basis where the quadrupole Hamiltonian is diagonal. The choice 
of STO or fictitious spin 1x2 operators CA. Wokaun and R. R, Ernst* 
J. Chem. Phys. 67* 1752 C1©77D; S. Vega* J. Chem. Phys. 

66 * BSl SC 1 07S!> 3 turns out to be very suitable for describing the 
dynamics of the evolution of spin 3x2 nuclear assembly subjercted 
to an rf field perturbation. We also introduce in this Chapter 
the quadrupole interaction frame CQIFD; the high frequency 
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approximation for first order perturbation theory calculations; 
the observables of NQR experiments* namely, the magnetization 
operators for spin 3/^; the density operator; the high temperature 
approximation for the equilibrium density operator; the 
Liouvilie-von Neumann equation; and the pulse propagator. Also we 
describe in this Chapter the roles of the internal Hamiltonians in 
modifying N<^ line shapes , 

In Chapter III are presented the single pulse response and 
spin echo calculations for spin 3/^2 NC^ in a system of 
non-interacting nuclei. The purpose is to demonstrate, perhaps 
for the first time, the effectiveness of the STO formalism in 
typical spin 3/^ N<SR calculations. We also demonstrate in this 
Chapter the applicability of STO formalism to analyse 
theoretically the zero-time resolution CZTfO technique for the 
recovery of free induction decay CFID) completely in a spin 3^2 
NQR system with non -axially s 2 >Tnmetric efg. This is an extension 
of the recent work of Singh and Armstrong CM. Singh and R. L. 
Armstrong, J. Magn. Reson. 78, S38 Ci988>] v/ho have applied the 

same to study an axially symmetric spin N<^. 

In Chapter IV, we go to the next level of complexity in spin 
3^ N<3^. Here we present a theoretical study of the multi pie -pulse 
sequences, namely. Phase Alternated Pulse Sequence CPAPSI> and 
WAugh-HUber -HAeber len CWAHUHAD' using Average Hamiltonian Theory 
CAHT) within the Magnus expansion framework. The cjrcllclty 
criterion, as it applies to a polycrystalline spin 3/^ NC^ solid. 
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is explicitly stated; this enables us to formulate several 
multiple-pulse cycles which are adapted to the study of the major 
interactions in quadrupolar solids. The nuclear assembly is now 
considered as being no longer a non-interacting one and* 
therefore, heteronuclear dipolar interactions are taken into 
account. The homonuclear dipolar interactions are very weak in 
spin 3/^ NQR and are discarded Also considered are the effects 
of impurity and torsional oscillations on the efg. A noteworthy 
result of this Chapter is the dependence of the effectiveness of 
the pulse sequences on crystal orientation and asymmetry 

parameter. Finally our experimental results on the effects of 

9SP 

PAPS on spin 3.*^ in powdered KCIO^ C Cl 5 are presented in this 
Chapter . 

The generalized AHT within the framework of Floquet theorem 
has been found to be most appropriate to analyse theoretically 
the results of Ostrof f-Waugh COW> or Mansf ield-WareCMW> 
multiple-pulse sequence experiments in spin 3./2 NC35. Chapter V* 
therefore, starts with a detailed exposition of the Floquet 
theorem, the generalized AHT, and the derivation of equilibrium 
properties of spin 3^3 NCS? therefrom, in the presence of 
periodic rf pulse perturbations. We also describe in this Chapter 
the advantages of the generalized AHT within the Floquet framework 
over AHT within the Magnus expansion framework and the conditions 
under which they yield similar results. The results of AHT in 
three formalisms, namely, Magnus expansion, Floquet theorem and 



xii 


Krylov-Bogoliubov-MitropolskyCKBMD method » are compared in this 
Chapter . 

Chapter VI gives the detailed results of our calculations 
performed for both si ngle--cr ystal and piowder samples of spin 3/^ 
NO? under the OW or MW pulse sequence. A Monte Carlo procedure 
based on '‘systematic sampling* technique was used for powder 
averaging, and the Fortran program which performs the powder 
average is listed in Append!.. A of the thesis. 

Theoretical plots of the calculations are compared with the 
corresponding experimental results performed on powdered HgCl^ and 
NaClO^C Cl resonance!) samples. Attainment of a quasi -“equilibrium 
state, viz. , a spin-locked state, is demonstrated, both 
theoretically and experimentally. The dominant mechanism 

responsible for achieving this quasi -equilibrium state is the 
coupling of quadrupolar spins to the lattice via the oscillating 
efg due to the torsional oscillations. Unlike in NMR, the 
homonuclear dipolar interaction does not play any significant role 
in spin 3 ^^ NQR. The consequence of the above quasi -equilibrium 
is the persistence of the signal at times much longer than ST^. 
This finding not only considerably generalizes and adapts the 
earlier work on spin 1 NQJ? by Maricq CM. M. Maricq. Phys.Rev, B33t 
4S01 C 1 0©6!> 3 , but also enhances our understanding of NQR phenomena 
f ur ther . 

In Chapter VII, we sununarixe the present work and the 
important conclusions drawn therefrom. We also point out the 
scope for future work in this Chapter. 



Thm experiments carried oat in this thesis have value from 
the standpoint of demonstration of theoretical predictions. In 
particular, the long-time behavior of the spin 3X2 
quasi -equilibrium magnetization has been unequi vocally- 
demonstrated. 
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CHAPTER I 


INTRODUCTION 


The subject matter of this dissertation is a study of the 
dynamics and equilibrium properties of a spin 3x2 pure nuclear 
quadrupole resonance CN<3RI> subjected to a scheme of radio- 
frequency CrtD radiation pulse perturbations. In magnetic 

resonance literature, the various schemes of rf radiation pulses 
are known collectively as multiple-pulse sequences. 

I. A. A Survey of the Essentials of Pulsed Experiments 

In most pulsed NOR experiments, a multiple-pulse sequence is 
applied to a sample containing quadrupolar nuclei Cspin 5: 1 in 

units of tD initially at equilibrium in the inhomogeneous electric 
field created by the surrounding electronic charge cloud. The 
typical experimental arrangement is shown schematically in Figure 
I.l. The rf pulses excite the nuclear spins to a non -equilibrium 
state, from which much weaker rf signals are emitted and recorded. 
That state may be , specified in detail by 


»*From here onwards, NQR will mean pure NQR, i.e. NOR with no 
Zeeman field, unless specified otherwise. 



Figure 1*1: Schematic experimental set up for pulsed NQR. A 
sample is placed in the coil of a tuned rf 
circuit. The i nhomogeneous electric field of the 
surrounding electronic charge cloud interacts with 
the quadrupolar nuclei to give rise to quantized 
energy levels. The nuclei in the sample distribute 
themselves among these energy levels according to 
the Maxwell -Bol tzmann distribution law, rf pulses 
create an oscillating field along the x-ax±s. 

Oscillating rf signalst emitted by the sample in 
the x—axis direction, are detected. The Cartesian 
coordinate system shown defines the laboratory 
f name. 




NQR Sample. 



Figure M 
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a density operator. However* there is often one property of this 
perturbed spin state that is of particular interest. This 
property may then be considered to be the response of the spin 
system to the sequence of the applied rf pulses. Figure 1.2 gives 
a schematic representation of the excitation process. 

The rf pulses used to excite the sample are* most commonly* 
well approximated by rectangular shapes. The height of these 
rectangles represents the amplitude of the rf field. A 
multiple-pulse consisting of n pulses is schematically written 
as: 


cd -T -ce 

t (p 1 z 


-T -ce > . -T 

Z B tp B 

a 


n <p 


where, e^^and .re the flip .ngle .hd ph.se of the ith poise. 

respectively. is the time that elapses between the ith and 

<;i-^l>th pulses- The flip angle is proportional to the duration of 
the rf pulse* and the name is reminiscent of the fact that* in the 
classical vector model of spin 1/^2 nuclear magnetic resonance 
CNMR>* it represents the angle by which the spin vector is rotated 
by application of the rf pulse for a certain duration [i-43 . 
Ideally, a certain pulse in the sequence is considered to be an 
impulse or Dirac’s <5— function. 


»*A11 the references appear at the end of the last Chapter. 
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Figure I^iSx An NQjR experiment, viewed as a measurement of the 
response of a sample consisting of quadrupolar 
nuclei C spins S: 1 in units of h> to a given excita- 
tion. The excitation is typically a sequence of rf 
pul ses wi t h wel 1 -def i ned 1 engt hs and phases » 

possibly separated by delays. 
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The response is determined by the set of parameters involved 
in the experiment. Indeed* the response is modified not only by 
the parameters that can be manipulated externally, but also by the 
ones that are intrinsic to the system that is 
studied. When the rf amplitude is kept constant, the external 
parameters that are involved in a mul ti pie -pulse experiment are 
the individual pulse lengths, rf -phases, inter' pulse time-durations 
and the resonance offset, a quantity defined as the difference 
between the r adi of requency and the resonance frequency of the 
spin. The parameters intrinsic to the NQR sample under study can 
all be put in two classes - Ci) the strengths of spin-couplings 
and Cii> the distribution in the electric field gradient CefgD due 
to impurities, lattice defects and torsional motion of the units 
constituting the lattice. Some of the spin— coupling terms give 
rise to fine-structure in the response signal, thus making it rich 
in information content, while others make it structureless and 

hence less informative. This masking of information content of 
the response is known in magnetic resonance literature as 

line-broadening. Here it is worthwhile pointing out that, by 
manipulating the external parameters available to the 
experimenter, it is possible to enhance the information content of 
the response. These experiments in magnetic resonance are known 
collectively as selective-averaging or line-narrowing techniques. 
Later on, we shall deal with the nr^chanlsms of 1 ine -broadening in 

N<3R in general terms and also the techniques of 
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selective-averaging in detail. The distribution in efg invariably 
gives rise to broadening of lines. 

The focus of our attention in this thesis is the MC2i? of spin 
3/^ systems. A central theme that has guided our study is the 
behavior of the response of the system right from zero-time to 
times of the order of - 8T^ . where is the epi n-Cqoadr epole5 

i ait ice relaxation time of the sample under si jdy. The 

establishment of quasi -equil ibrium state under varying conditions 
of the external parameters has been studied theoretically. 
Furthermore, vherever possible, appropriate experiments have been 
performed on model quadrupolar spin 3^^ solids to seek support 
for our theoretical findings. Very interestingly, it turns out 
that the mechanism of this quasi -equilibrium state leads to a 
characteristic prolongation of the quadrupolar relaxation time 
, similar to previous findings on the prolongation of the 
spin-spin relaxation time CT^> in the case of pure dipolar 
solids. In add! t ion to the above mentioned focal theme » our 

studies have enabled us to adapt several multiple-pulse cycles 
which have been in use for studies of dipolar solids ; the 
effects of these cycles on the various major interactions in 
quadrupolar solids are studied theoretically. 

Before we go into a discussion of the importance and 
usefulness of the present work,, it may be worthwhile to point out 
the salient features of NQR vis-A-vis NMR in general, and how the 
spin 1/^ NMR, spin 1 NQR and spin 3/^ compare with each other 

in particular. For the ensuing discussion, we have heavily drawn 
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upon the relevant material from standard monographs Cl '-43 and 
research papers C5-*-73. 

The quantized energy levels in NMR have their origin in the 
magnetism v^hich the nuclear spins possess inasmuch as they are 
created by the interaction of the nuclear magnetic moments with 
the applied static magnetic field. On the other hand, the 
energy levels owe their origin to the electrostatic interaction of 
the electric quadrupole moment C which is due to the non'-spher ical 
charge distribution in the nucleus> with the ef g due to the 
surrounding electronic charge cloud- This is a tensorial 

interaction C5-73 . 

The principal axis system of the efg tensor is chosen such 
that its diagonal components, V , V and V satisfy the 
condition [83, 



ci.iy 


Here V is the electrostatic potential at the nucleus due to the 

surrounding electronic charges and V = ; 1, j = x, y, 

la a ^ 

It is assumed that the electric field at the site of the nucleus 

is produced entirely by charges external to it, in which case the 

Laplace equation, V •♦‘V -+•¥ «0, is satisfied- A measure of 

XX yy 22 I 2 Z: 

the deviation of the efg from cylindrical symmetry is given by a 
parameter r?- In the NQR literature, t) is known as the asymmetry 
parameter and is defined as 
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V -.V 

1 XX VY J 

y . CI.2> 

The two constraints* Eqn I.l and the Laplace equation* lead to the 
result that the value of t) can vary between O and 1 only. n is 
zero for the axial symmetry of the efg tensor around z-axis of the 
principal axis system. In most chemical systems* the principal 
axis system is fixed in the molecule or crystal lattice which 
contains the quadrupolar nucleus or nuclei and is known in NCS^ 
literature as the quadrupole principal axis system tQPASO. 

The non -spherical charge distribution is possible only in 
nuclei with spin > 1. We therefore have NCy? of spin 5: 1 only. 
This is unlike NMR where spin is the most studied case. 

In NMR the external static magnetic field defines the 
quantization axis for all the nuclei constituting the system. The 
situation is quite different in NQR. The quantization axis in this 
case is intrinsic to the system and an experimenter has no control 
over it- It is defined by the main component of the efg tensor, 
namely* the QpAS. Wnetx we have a powder sample* which 

happens to be the case most commonly in NCS?* each crystallite 
constituting the powder will have its own quantization axis and 
an external rf field perturbation will not have the same effect on 
nuclei in different crystaLli tes. This one point makes the 
theoretical treatment of an NCS? problem quite different from that 


of NMR. 
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Once the energy levels are created* the external rf field 
perturbation in a plane perpendicular to the quantization axis 
induces resonant transitions among them. In both NMR and NQR these 
transitions are magnetic dipole induced tr ansi ti ons * and hence the 
selection rule in both the cases is in general* Am - ±1 * where m 
is the magnetic quantum number for spin I and has values ranging 
from --I through *^1 in steps of unity. Seen from this point of 
view the term Nuclear Quardrupole Resonace is a misnomer f 73 . It 
is interesting to mention here that * recently * using a dc 
superconducting quantum interference device CSQUID7 detector Hahn 
et ai [tO~-lS3 have been able to observe the effect of electric 
quadrupole induced transitions in spin 3/^2 nuclear assembly a 
phenomenon that is due to nuclear quadrupole resonance in the 
literal sense of the phrase. 

The NMR energy levels for any spin I are determined by 

E = mxB * Cl . 3> 

m ^ o 

where y 1^ the magnetogyric ratio and B the strength of the 

o 

static magnetic field applied externally. The number of energy 
levels are* therefore* 21 1 and they are ail non-degenerate. It 

is also clear from Eqn 1.3 that the energy difference between 
any two adjacent levels is determined only by y and Therefore* 

the levels are equidistant* Ip the particular case of spin IxS 
NMR* there are just two energy levels. 

In the most simple case Cr) 07* the NQR energy levels for 
any spin I >1 are given by 



wher e , A 


la 

E ■= A [3m^ - ICI + i:>3 » Cl. 45 

m. 

*= e qQ^'^41 C SI ~1 !> , Q is the quadrupole moment of the 
nucleus, eq is the magnitude of V and e is the electronic 
charge. e^qQ is known as the quadrupole coupling constant CQCO, 
The above expression holds both for hal f -^int -grai and integral 
spins. Here the energy levels are doubly degenerate in m, because 
the tvso states }-+-ni> and j-m> have the same energy as is evident 
from Eqn 1.4 , Also, the energy difference between the adjacent 
pairs of levels will not be the same for various pairs of m 
values. TTierefore, NC^ energy levels are non-equi distant. For 
half -integral spins there are non-equi distant energy levels 

all doubly degenerate, while for integral spins there are I •+■1 
non-equi distant energy levels where I of thena are doubly 
degenerate and only one of them with m = O is non-degenerate. It 
may be mentioned here that, for integral spins, when O, the 

degeneracy is lifted completely, whereas for half -i ntegr al spins 
it is not. This follows from a special case of a general theorem 
due to Kramers Cl 31, Figures 1.3 and 1,4 illustrate the points we 
have discussed above. 

In the general case, therefore, the spin 1 NQR S 3 rstem will 
have a non-degenerate three energy level pattern, la±>eiled by the 
states |0> and j-l>. Accordingly, a rf field applied in any 

arbitrary direction gives rise to three transitions, namely, 
l+l>-fr>fO>, |0>-f-j-l> and f+l> f-|-l>. The first two form a doublet 
wi th f requenci es C i n rad sec 



13 


3e*qQ 

Cl ± r)/3D . Cl . SZ) 

The third transition ^ gives rise to a line at a very 

low frequency and is not of much interest. The intensities of the 
two members of" the doublet are proportional to \<^i j0> { and 

iO> j* .respective! V. Here Sf ^ is the Handltonian 

* ‘ r I ‘ " r 1 

representing the radiof r equency field perturbation. The higher 
frequency line has a maximum intensity when the rf field lies in 
the X““direction of the QPAS» and the lower frequency line then has 
zero intensity. The apposite is true with the rf field in the 
y-“di recti on of the QPAS. In other words > the spin 1 HQR resonant 
transitions are polarized C Figure I.4Ca>>. 

The situation in spin 3/^lS NQf? is quite di fferent. A1 though the 
presence of rj leads to mixing of states jm> differing by Am 
±2* the energy eigenvalues are still doubly degenerate and are 
gi ven as 


; ci.6a:> 

1X2 

^ - CI.6b5 

If an rf field is now applied, it is generally possible for four 
transitions to occur between eigenstates jni> and lm*> because of 
mixing of m states. The frequency for all the four transitions. 
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Figure 1^3: CaD and CbZ>: Sinvple energy level diagram for NMR of 
spin I in a static magnetic field CBJ), Ca> is for 
half -integral spin and Cb> is for integral spin. 
All the adjacent levels are equidistant, as seen. 
Cc]) and CdD: Simple energy level diagram for NC^ 
of spin I in a cyl i ndri cally symmetric electric 
field gradient CefgD. Cc> is for hal f — i ntegral 

spin and Cd!> is for integral spin. The quantization 
axis in this case is defined by V component of 

2525 

the efg in quadrupole principal axis system CQPAS>. 
All the adjacent levels are not at equal difference 
in these cases as seen. Thereforejin these cases we 
have a degenerate C except m = 0> and non~equi dis- 
tant energy levels. Energy differences in all the 
diagrams are in arbitrary units. 
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Figure 1.4: Ca> Spin 1 NQR. For t) — O, we have a single resonant 
transition. For ?} ^ O, there are three resonant 
transitions possiblev as explained in the text. These 
transitions are linearly polarized. 

Cb!> Spin 3/^ NQR. For both -0=0 and t) O only a 
single resonant transition is possible as shown. This 
transition involves all the energy levels and is not 
linearly polarized. 
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namely, l+l/^> <- j •+3/^> , 1 +1 /^a> |-l/<2> <-|+3/2> and 

1 -1 ^ ( -3/^> , is given by CFigure I.4CbZ>I>. 

6A i/z 

^ * — — Cl + . ci.7:> 

cat 

The transitions other than those corresponding to Am - ±1 are 

forbidden for axial symmetry and are expected to be weak except 
for large values of rj. The transition frequencies in Eqn 1.7 do 
not depend on the direction in the x-y plane in which the rf field 
is applied. Unlike spin 1 NQR^ therefore » the transitions are 
unpolarized in this case. 

We take the example of a sample of polycrystalline spin 1 NOR 
system. In a typical NOR experiment, it will be irradiated with an 
rf field of fr&a^u&nay that matches with one of the resonant 
frequencies of the spin. Let us also assume that the frequency of 
irradiation, in this example, corresponds to as given in Eqn 

1.5 and Figure I.4Ca5. If the rf field is applied such that the 
oscillating ma.gnetic field direction of x-axls of 

laboratory frame, the system under consideration can closely, if 
not exactly, be made to reseaxble spin NMR. This is because, of 

all the crystallites in the polycrystalline sample, only those 
that have their OP AS coinciding with the laboratory frame will 
favor their constituent nuclei to undergo resonant transitions. 
Only crystallites having this orientation in the powder are, 
therefore, excited. Also the rf field will connect only one pair 



of energy levels out of three, namely, the one that has energy 
difference proportional to the angular frequency Under this 

circumstance, the spin 1 NC2R system is Virell approximated by 
assuming it to be a two-level quantum system. The two energy 
levels of this system are the ones that are coru'iected by resonant 
rf field irradiation. This also follows from the fictitious spin 
1/2 formulation of the problem [143. According to this 
formulation, in a multi-level quantum system, when we connect only 
two of the levels by an external resonant perturbation, we can 
model the system as a two-level system to an excellent 
approximation, disregarding non-resonant levels completely. 

This however, cannot be the case in spin 3/2 NQR, where a 
single frequency radiation will connect all the four levels. 


I.B. Ii^ortance of Pulsed NQR of Spin 3/2 Nuclei 

Each of the quantum systems described above is representative 
of a class. An independent study of dynamics and equilibrium 
properties of each one of them, namely, spin 1/2 NMR, spin 1 NOS 
and spin 3/''2 NQR under a scheme of rf perturbation pulses deserves 
merit in its own right. As pointed out earlier, spin 1/2 NMR is 
the most widely studied case [15-193. While there also exist a few 
reports of the multiple-pulse :^tudies in spin 1 NQP [20-333, there 
is a need to study mul tipi e— pulse spin 3/2 NQR in more detail both 
theoretically and experimentally. The few studies that have been 
carried out in spin 3/2 systems will be described in section I . E. 
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From the point of viev<r of applications of NQR, the spin 3^ NQf^ 
has an additional importance. Only three nuclei out of 

eighty-seven naturally occurring NQR-active isotopes have spin 
l.In contrast » spin 3/^ nuclei are thirty-two in number, A 
theoretical and experimental investigation of the kind presented 
here» therefore, may extend the benefits of a vital domain of 
techniques to a large variety of NOR systems. 


!•€# The Usef«jJLness of Multiple-Pulse Experiments 

The multiple-pulse experiments are very useful for measuring 
transient properties of spin systems both in NMR and In 

particular! they are extensively used for the measurement of 
relaxation times [34-301, for the study of difCusLan processes 
C3Q3 and for the investigation of chemical reactions {40-4S1. 
Multiple-pulse experiments in solids, carried out on time scales 
of the order of T , show how the various factors contributing to 
the width of the resonance line respond to strong rf field 
perturbations. Also, these experiments reveal pulse sequence 
characteristics that lead to the narrowing of resonance line 
widths, thus yielding valuable information on the interactions 
leading to fine structure of the resonance lines [15-191. Some of 
the multipie-pulse sequences, such as spin-locking sequences 
of which the Ostroff -Waugh COVD C433 or Mansfield-Ware CMVD [443 
are examples, improve substantially the effective sensitivity of 
the NMR and the signals if the individual echoes within each 
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spi n-1 ocked-spi n-echo CSLSEI> sequence [20,43-46} are coherently 
added. This sensitivity improvement accrues even without our 
taking recourse to the complicated double resonance schemes 
C 47-503. So versatile are the multi pi e-pulse experiments that 
they can be used nowadays for solid state NMR [15-193, high 
resolution HMR in liquids [51-543 and for HQR C553. Indeed , such 
techniques, often called spin-gymnastics, have substantiality 
extended the scope of magnetic resonance. 


A Review of Multi pie -Pulse NQR 
The OW C433 or the MW [443 sequence, which we shall 
henceforth call as spin-locking sequence, is in general : 


<p {p'^n^z 


It represents one of the first uses of pulsed rf fields for 
studying a powder system of spin 1 in poiycrystalline 

at 77 K;> [203. This sequence also represents one of the first 

uses of pulsed rf fields for selectively suppressing internal 
spin-interactions that cause homogeneous line-broadening of the 
resonance lines, thereby increasing the resolution of NMR lines in 
solids [43,443. This particular sequence causes a refocusing of 
the transverse magnetization for periods much longer than T^. In 
some respects, the effect is similar to the more familiar 
spin-locking by a constant amplitude continuous wave CCWD rf field 



[563. In the case of NQR also, the effects are very similar to 
those observed in the NMR experiment, namely, the persistence of 
the magnetization for long times compared with T^. Initially, this 
was not conceptually reconcilable in light of differences between 
the origin of the interactions that give rise to the phenomena of 
NMR and of N<^ in nuclear assemblies. The average Hamiltonian 
theory CAHTD within the framework of Magnus expansion 
C 15-19, 57,503 has been found to be a powerful theoretical tool to 
deal wi th anal ogous si tuati ons i n spi n 1 ^ NMR. However , the AHT 
technique was thought to become intractable in spin 1 NQR case 
[£1,223. This is because the overall ti me -dependent Hamiltonian 
in the appropriate interaction representation in an NQR system 
is, in general, not cyclic in the sense of the system returning 
to its original state after the pulse sequence has been completed 
Csee Chapter IV>. This situation arises in a powder sample, as 
pKDinted out earlier, because the rf field has random orientations 
with respect to the individual QPAS of the crystallites 
constituting the powder. An explicit time evolution calculation 
[21,223 on a model system of dipolar coupled spin 1 nuclei, 
however, confirmed the refocusing of the spin magnetization in the 
multiple-pulse NQR experiment of *^N nuclei in NaNO^ powder 
sample, in the short-time regime * Later on, Maricq C333 

gave a detailed theoretical analysis based on the generalized AHT 
within the framework of Floquet theorem C 59-64 3 . A1 though the 
spin-locking sequence does not remain cyclic independent of the 
orientation of the rf field, Maricq has shown that the average 
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Hamiltonian approach can still be taken if one considers a 
single-crystal sample. In fact» this also has the advantage of 
revealing effects which are obscured in powdered samples. The 
result of single— crystal calculation can then be subjected to a 
powder -aver agi ng procedure, taking into account the proper 
weightages for various orientations of the crystallites in the 
powder. Mar icq also took into account the many-body nature of the 
interactions in the real system CNaNO^J in his theoretical 
analysis and successfully accounted for both the short C-^ T^Z> and 
1 ong-t i nieC ~ behavior of the overall magnetization. 

Incidentally, we mention here that earlier, generalized AHT was 
used by Maricq to settle the issue of the validity of the AHT for 
theoretical anai 2 rsis of the mul tipi e -pulse experiments in general 
C64] ; furthermore, generalized AHT has also been used for the 
theoretical analysis of spin 1-^ NMR multiple-pulse experiments, 
namely, OW or MW and WAugh-HUber -HAeberlen CWAHUHAD [64-663. 

In mul tipi e -pulse experiments we confront not only a 
many-body system, but one which is driven by time-dependent 
external forces. Because the Hamiltonian corresponding to such a 
system is energy non-conservative, one cannot apply the normal 
methods of statistical thermodynamics to determine the equilibrium 
properties of the system. Maricq £653 showed how the Floquet 
solution for the time -dependent Schrddinger equation can be used 
to extend statistical thermodynamical ideas to systems evolving 
under a periodic ti me -dependent Hamiltonian. This was used by him 
to analyse theoretically the equilibrium properties of the overall 



magnetization in spin 1/S NMR C0S3 and spin 1 NCS5 C33] driven by a 
scheme of rf field. perturbation sequences. Provotorov and 
co-workers [24,26,292 have also given a parallel theoretical 
analysis of the response of spin -locking sequence in a spin 1 
Certain pulse sequences such as phase-alternated pulse sequence 


CPAPS5 [67, 

683 

and 

WAHUHA 

remain 

cyclic irrespective of 

the 

orientation 

of 

the 

r adi of r equenc y 

field in the QPAS of 

each 

crystal li te 

in 

a 

powder . 

Osokin 

[23,25,27,30—323 has 

taken 


advantage of this special property of these multiple-pulse 
sequences, and has used the AHT within the framework of Magnus 
expansion to discuss the effects of these pulse sequences in spin 
1 NOR experiments of in NaNO^ polycrystalline sample. 


I.E. Motlvatiori for the Present Work 

All of the foregoing discussion has been about an NO? system 
of an assembly of spin 1 nuclei, where the multiple-pulse 
experiments have been carried out on one of the three possible 
resonances. Subsequent to the above work, there have been reports 
of experiments, from our laboratory [69-713 as well as elsewhere 
C723, on spins 3/2, S/2 and 7/2. The pulse sequences employed by 
these authors can be broadly classified into four categories; 

CiD 4f±+nyz'~^^ri ~ spin-locking sequence, 

CiiD Crr/2> , -2 t -C , -t} - 

<p n 


the PAPS 
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-ir-cey ^ -ar-cey 

^ ' ~<^ 4f-i-Tzyz -<4>+nyz} 

-T-C6>y^~r} - the WAHUHA , 

& ci^y cn^y-ir-c&y-T-c&y ^ -ar-c&y.^ ^ -r-c&y -ar-c&y . 

' <P - < ip-i-n/'z ) -<p —<p 

-T-C&y . ^ Sr-COy . ^ . —r-coy.-r} ~ the Mansf i ei d-Rhi m- 

-<^+n :/2 > <p+n^z <i) n 

Elleman-Vaughan CMREV-8:> [73-773, 

where ^ =0 when rf pulses are applied in the x-direction of 
the laboratory frame. In the above, we have intentionally kept 
the phases Cify of the pulses arbitrary. Unlike spin t^a NMR, as 
will be shown in Chapter II, in NQ2 ^ loses its usual meaning of 
representing directions in the rotating frame . Ainbinder et al 
[723 have examined the spin-locking effect under the influence of 
Carr— Pur cell -Mei boom— Gill CCPMGI> sequence [783 both theoretically 
and experimentally in spin 3ya N<^ C Cl in polycrystalline 
samples of SbCl , C Cl and CCl COOH at 77 KD and also in spin 7/^8 

9 Z 3 

NQR t^'^Sb in polycrystalline SbCl^ sample at 77 K!!) . Narsimha 

Reddy [69,703 has performed expert n;^ents using the pulse sequences 
of the categories Ci> and Cii'^ listed above, with varying values 
of the pulse parameters 8, <fi, r and offset from resonance 
conditions. His most important observation was that in the 
powdered sample of a spin NQR system C Cl in NaC3 0^ 


and 
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KCIO^ polycrystalline sample at room temperature?, spin-locJcing 
effect is observed for any value of i>y the phase of the pulses in 
the sequence with respect to that of the preparatory first pulse. 
His experiments included well-known sequences such as OW or MW, 
PAPS, Carr-Purcell CCP? £35] and CPMG. Ravinder Reddy £711 
reported experiments both on spin 3/^ and spin 5-^ NQR systems. 
His multiple-pulse sequences were from the first three categories 
listed above. His spin 3/^ experiments were on Cl in the 
samples of NaClO^, KCIO^ and a^Cl^. While experiments on NaClO^ 
were done both on a single crystal and a powder, only powder 
samples were used in case of KCIO^ and SbCl^. For spin 5/^ 
experiments, a powder sample of KReO^ was used and the nuclei 
studied were **^Re and ^*^Re. In all the spin-locking 

experiments, there is strong evidence for establishment of a 
quasi -equilibrium state for the overall magnetization under the 
spin-locked state. It means that the transverse magnetization, 
under mui ti pi e-pul se sequences of the type mentioned above, 
remains in the x-y plane of the laboratory frame for times much 
longer than the characteristic relaxation time of the system . 
The experiments in £711 using WAHUHA also show lengthening of the 
characteristic relaxation rate constant. The theoretical analysis 
of these experiments as well as of MREV-8 £713 is a preliminary 
one and needs further investigation. 

The other kind of multiple-pulse sequence that has been used 
to narrow the NQR lines is Lee-Goldburg CLG? £703 sequence. The 
possibility of application of this sequence has been discussed by 
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Pratt et al [80,81]. Mefed and Pavlov [82] used this sequence on 
I in a Cdl^ single crystal Cl = Their experiment was 

performed on the resonance line corresponding to the 

ILZT 

ir ansi t i on 1 ±i/ 2 > 4 -^ 1 ±3/'E> in I nucl ei . They obser ved 

line-nar rowing of the* NQR line in a rotating reference frame that 
rotates with one of the components of the oscillating rf field. 
The free induction decay CFID!> in this case was found to be -- 25 
times slower than that in the laboratory frame. In his 

theoretical analysis^ Pratt [81] considered dipolar couplings and 
distribution of electric field gradients due to static 
imperfections in the crystal lattice, to be the line-broadening 
interactions. His conclusion was that a multiple-pulse technique 
for partial suppression of the dipolar couplings followed by a 
spin-echo method for eliminating the i nhomogeneous broadening 
would lead to enhanced sensitivity and resolution in N<^. 

An alternative theoretical method of analysis of the response 
of LG sequence in spin nuclear assembly has been used by Zueva 

[S3]. In this procedure, the Hamiltonian is transformed to a 
generalized interaction frame, where a second moment calculation 


of the 

NQR line shape in 

the presence of 

the 

LG 

sequence 

is 

carried 

out . The second 

moment i s shown 

to 

be 

smaller 

in 


magnitude when the LG sequence is used. The line -broadening 
interactions considered in this calculation are also the dipolar 

couplings and the distribution in efg. 

Zueva and Kessel [84—861 have also reported a theoretical 
analysis of the response of a multi pie -pulse sequence from an NC^ 
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system of any half -integral nuclear spin assembly. In this 
analysis* the axially symmetric systen.\s Ciq«=OI) were considered. The 
theoretical analysis used was the same as has been described above 
for spin 3/lS systems* namely* the second moment calculations in 
generalized interaction representation for the line shape* in the 
presence of line-broadening interactions. These authors have also 
considered the following multiple-pulse sequences: 


Ci5 Cn^:) , -aT-C6^:> -ar-COD -tI , 

o o o o n 


& CiiZ> Cn/^y , -{r-C^y -SaT-C<9> -8o(t-C0> 

rr /2 ‘ o o o o rrxz 


-aotr-coy ^ -ar-cey ^ -aar-cey ^ - tj - 
rr/'z rr /^2 n/'z n 


Ca — a positive integer!). 

They have been able to show that these pulse sequences ajre also 
very effective in line-narrowing of NQR lines in hal f -i ntegral 
nuclear spin S 3 rs terns. We may also mention here that Ainbinder and 
Furma.n EB7] have given an independent theoretical framework based 
on Kryiov-Bogoiiubov-Mi tropolsky CKBM> C 18*88,893 averaging 
procedure and the canonical transf ormation technique to analyse 
multiple-pulse experiments for a general spin I with non- 
equidistant spectra. This theory offers an average Hamiltonian 
which is similar to the ones based on Magnus expansion 
[15-19,57*583 and on Floquet theory [643 in the lowest order 
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approximation. For higher order approximations, the differences 
start appearing. This theory, however, remains largely 

unexploited till date in magnetic resonance. 

More recently, there have been reports on both experimental 
and theoretical studies of the kind where a modified 
Jeener-Broekaert CJBl sequence C903 has been used to transfer 
quadrupolar order into dipolar order in spin 3/2 pure NQR with 
line-broadening due to heteronucl ear dipolar interaction 1913 of 
®^Cl-^^Na in NaClO^ and due to efg 1 nhomogenel ty C923 in K^OsCl^. 

guy 

In both the cases the nucleus studied was Cl. Theoretically, 

the introduction of the concept of spin temperature, cor respondi ng 

to the Haxni 1 tonians for internal interactions, has been very 

useful in predicting the correct experimental behavior of the 

overall magnetization. Also, the zero-time resolution CZTR> C931 

technique has been applied recently to recover the FID completely 

in an NQR experiment of spin 3/2 nuclei C Cl in K OsCl 5 in. 

z <s 

axially S 30 iiittetric case [941. 


A Preview of th.e Present Work 

This study starts with the choice of a suitable operator 
formalism to describe the dynamics of the evolution of spin 3/^ 

nuclear assembly subjected to an rf field perturbation. The 
generators of the SUC4> Lie algebra [95,963, known in magnetic 
resonance as single transition operators [973 or fictitious spin 
1/2 operators [983 , provide a suitable choice. After introducing 
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all the Hamiltonians needed in the usual Cartesian basis of 
angular momentum operators, we present, in Chapter II, these 
Hamiltonians in terms of the single transition operators CST02> in 
a basis where the guadrupole Hamiltonian is diagonal. This is done 
for the mathematical convenience of the calculations in later 
Chapters. We also introduce in this Chapter the quadrupole 
interaction frame CQIF;>; the high-frequency approximation for 
first o^der perturbation theory calculations; the observables of 
NQR experiments, namely, the magnetization operators for spin 3/2; 
the density operator; the high temperature approximation for the 
equilibrium density operator; the Li ouville-von Neumann equation; 
and the pulse propagator . In brief, this Chapter provides the 
quantum mechanical background for the work presented in later 
Chapters. Also we describe in this Chapter the roles of the 
internal Hamiltonians in modifying N<^ line shapes. 

In Chapter III are illustrated the single pulse response and 
spin echo calculations for spin 3/2 NC3jR in a system of non- 
interacting nuclei. The purpose is to demonstrate, perhaps for 
the first time, the effectiveness of the STO formalism in typical 
spin 3/2 NCy^ calculations. We also demonstrate in this Chapter 
the applicability of STO formalism to analyse theoretically the 
zero-time resolution technique for the recovery of FID completely 

in a spin 3/2 NQR system with non-axially symmetric efg. This Is 

■% 

an extension of the recent work of Singh and Armstrong [941 who 
have applied the same to study an axially symmetric spin 3/2 


NQR. 
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In Chapter IV, we go to the next level of complexity in spin 
3 ^ NC2?- Here we present a theoretical study of the 

multiple-pulse sequences, namely, PAPS and WAHUHA using AHT within 
the Magnus expansion framework. The cyclicity criterion , as it 
applies to a p>olycrystalline spin 3x2 NC3? solid, is explicitly 
stated; this enables us to formulate several multiple-pulse 
cycles which are adapted for studying the major interactions in 
quadrupolar solids. The nuclear assembly is now considered as 
being no longer a non-i nteracti ng one and, therefore, 
heteronuciear dipolar interactions are taken into account. The 
homonuclear dipolar interactions are very weak in spin 3X2 NO? and 
are discarded. Also considered are the effects of impurity and 
torsional oscillations on the efg. A noteworthy result of this 
Chapter is the dependence of the effectiveness of the pulse 
sequences on crystal orientation and asymmetry parameter. 
Finally, our experimental results on the effects of PAPS on spin 
3x2 in KCIO^ <I^^C1!> are presented. 

The generalized AHT in the framework of Floquet theorem has 
been found to be most appropriate to analyse theoretically the 
results of OW or MW pulse sequence experiments in spin 3x2 also, 
as in spin 1X2 NMR and spin 1 NQR. Chapter V therefore starts 
with a detailed exposition of the Floquet theorem, the 
generalized AHT, and the derivation of equilibrium properties of 
spin 3X2 NCJjB therefrom, in the presence of periodic rf pulse 
perturbations. Ihe theory presented in this Chapter is in the 
same spirit as that of Maricq £33,64-66] but for the spin 3x2 N<3R 
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experiment itself. We also describe in this Chapter the 

advantages of the generalized AHT within the Floquet framework 
over AHT within the Magnus expansion framework and the conditions 
in which they yield similar results. The results of AHT in three 
formalisms, namely, Magnus expansion, Floquet theorem and 
Krylov-Bogol iubov-Mi tropoisky CKBM> method, are compared in this 
Chapter . 

Chapter VI gi \/es the detailed results of our calculations 
performed for both single-crystal and powder samples of spin 3/B 
under the OW or MW pulse sequence. A Monte Carlo procedure 
based on the * systematic sampling* technique was used for powder 
averaging, and the Fortran program which performs the powder 
average is listed in Appendix A. 

The theoretical plots are compared with the experimental 

3 S 

results performed on powdered HgCl^ and NaClO^ C Cl 

resonance^ samples. Attainment of a quasi -equil i bri um state, viz. , 
a spin-locked state, is demonstrated , both theoretically and 
experimentally. The dominant mechanism responsible for achieving 
this quasi -equilibrium state is the coupling of quadrupolar spln.s 
to the lattice via the oscillating efg due to the torsional 
oscillations. Unlike in NMR, the homonuclear dipolar interaction 
does not play any significant role in spin 3/^ NOR C 9©. 1003. The 
consequence of the above quasi -equilibrium is persistence of the 
signal even at times much longer than 5T^- This finding not only 
considerably generalizes and adapts the earlier work on spin 1 NQR 
by Maricq [333 , but also enhances our understanding of NC3? 
phenomena further. 
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In Chapter VII, we summarize the present work and the 
important conclusions dravm therefrom. We also point out the 
scope for future work in this Chapter. 

Overall, the thesis has been written up in a style that 
blends the quantum mechanical and magnetic resonance languages, 
with the Jargon of the latter predominating somewhat. This is 
me nt to be in keeping with our belief that current problems of 
NCS$ dynamics and equilibrium properties are topics of research 
with an intrinsic interest. On the other hand, certain details 
which may seem somewhat elementary have been treated at length in 
this thesis, and this has been done intentionally to make up for 
the absence of such treatment in currently available monographs, 
textbooks and other dissertations. 

The experiments carried out in this thesis have value from 
the standpoint of demonstration of theoretical predictions. In 
particular, the long-time behavior of the spin 3/^2 
quasi -equilibrium magnetization has been unequivocally 


demonstrated. 



CHAPTER II 


THE QUANTUM MECHANICAL BACKGROUND 

II, A. The Hamiltonian of the Sjrstem 

The interactions of nuclear spins that occur in a quadrupole 
nuclear assembly can be described by a Hamiltonian of the 
following general form Cl -43: 

ffC = X ffC, , CII. 

Q r f t rit ^ 

where we have also Included the interaction of the quadrupolar 

system with the extremal radi of requeue y field perturbation. 

has units of radians par second or energy divided by h. the 

quadrupole Hamiltonian, corresponds to the interaction of the 

electric quadrupole moment of the nucleus with the surrounding 

electric field gradient Cefg3> due to electronic charge cloud. 

^ , known as the radiofrequency CrfD Hamiltonian, describes the 
rf 

interaction of the quadrupolar nuclei with a linearly oscillating 

rf field. 9 € stands for internal Hamiltonians, and it 

ir.1 

represents all other Internal interactions of the nuclei that 
constitute the NQR system. The internal interactions include 
coupling between pairs of spins and between spins and internal 
oscillating fields and field gradients. 



The Hamiltonians take the most simple form in their own 
principal a>d.s system CPASD. To begin with therefore, it is 
convenient to define each of the three terms on the right hand 
side of Eqn II. 1 in a different coordinate system. We adopt 
the following notation for this purpose: 

Coordinate system 
Quadr upol e PAS C CSP’AS^ 

Laboratory PAS CLABPASC? 

Internal Interaction PAS CIPAS> 


Notati on 
X, Vj z 
X, y, X 
x’ , y ’ ,2* 


The Quadruoole Hamiltonian: 96 
^ 

The OPAS is defined by the principal components of the efg 
tensor as described in Chapter I Cp.9>. The single-particle 
quadr upol e Hamiltonian in this coordinate system is given by 
[1-43 ; 


96 


Ol 


£31* - 5 tl* + !*:>] 

6 *• z £ + - ^ 


cii . s:> 


c»> 


V6 


r T 


20 




t'^ cs:»] 

ZZ ^ 


CII . aa> 


Is in frequency units and is given by: 


Qt 


3e*qQ 


6A 


£icai-i:>h 
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A is the same as defined in Eqn 1.4. and CS^ in Eqn II. £a 

20 ZZ 

are the rank two irreducible spherical tensors: 

cs:) - ] 

ZZ 22 2-Z-* 

For spin 3/^ HQR, co^ is the one and only resonance frequency. T) 
is the asymmetry parameter defined in Chapter I Cp.9»10Z>. I'‘s are 
the components* or linear combinations thereof* of the total 
angular momentum operator of the nuclear spin. 

The Radi of reguenc v Crf!> Harni 1 toni an : 

T1 

Here we assume the classical field form of radiofrequency 
radiation. This form is quite appropriate for describing most of 
the magnetic resonance phenomena. The interaction of the 

radxof requency radiation field with a nuclear spin is described by 
PC . We place the axis of the rf coil along the y-axls of LABPAS 

ri 

defined earlier. The explicit form of the single particle in 

LABPAS Is then ; 

PC ^ = “--CO C t]>cosCa>t 
rf 1 y 

= --'ico C t>cosCoo«t-^^]>T*'^CSI> . CII.3a> 

1 li 

oy^CtD equals ^B^CtO, where y is the magnetogyric ratio of 
the nucleus that is being irradiated. B^CtD is the rf magnetic 
field amplitude. ci> and tp are^ the frequency and phase of the rf^ 



respectively. In multiple-pulse experiments, typically the rf 
field is applied in pulses that are ideally rectangles, sq that 
t]) is piecewise constant, taking on the values O and only. 
0 l'° will be referred to as the nominal rf aniplitude. Also, 
typically co is constant and 4> is piecewise constant. The first 
pulse in these experiments, known as the prepaj'atory pulse, has 4> 
= O. Other pulses in the sequence will have different values for 0 
depending on the kind of expert n>ent being performed. Some of the 
examples of the form of multiple-pulse sequences have already been 
given in Chapter I Cp.£i,£4 and ££> . = Cl.xys> £ ] 

in Eqn 11.3a is a first rank tensor. 

The Internal Hamiltonians : Sf 

The internal interactions that play a role in the dynamics of 
the NQR are of different physical origins. We describe belov/ the 
Hami 1 toni axis cor respond! ng to each one of them. 

1. The efa i nhomoaexiei tv Hamiltonian: D SC 
^ ^ 

The single particle efg i nhomoge?iei ty Hamiltonian in the QPAS 

i s r 




D se 
<i} a 


D 


ijS 


a 


6 


[31 


SK 








ci^ + 




[ T 


20 


-rr CSOl 
VS 22 ^ 


CII . 4) 


= D -rr 
w Vo 


CII . 4&J 
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It causes inhomogeneous broadening of NCS? lines. is a function 

that is introduced phenomenologically to take into account the 
distribution of resonance frequencies around the centre frequency 
This distribution of frequencies is the result of the efg 
i nhomogenei ty C 102-1053. Dislocations and strains in the crystal 
or powder grain cause random distortions in the i nter mol ecul ar 
interactions and a corresponding random cistrlbution of electric 
field gradients at nuclear sites. Strains may result from the 
crystal growth process, cold work or from cooling the sample 
under study. Disorder in the crystal also causes a similar 
distribution. The most appropriate distribution function in this 
case is the Gaussian distribution function: 

1/2 

gCw-w^:) = [ exp-Ca>-0J.^:>^/C26*:>] /C2 »tc5*> . CII.5> 


where 


J-oo 


gC D d<i> 
a. 




and 

i 



2 


C A(a.O 

ty'Z 

. ^ 

C£ln2:> 


Cd 


a 


CII. 7> 


is the root-mean'-square deviation in A<a> and CA€v>5 is the half 

width of the NQR line at half intensity. T** is the effective 

** 2 

spin-^spin C transv^rse> relaxation time. o> is a frequency away 

from the centre frequency co • 

" a 
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£. 7>>e Torsional Oscillation Hamiltonian: ffC 

— — ^ 

In the NC95 sample, molecules containing the quadrupolar 
nucleus undergo torsional motions about the three principal axes 
oT the field gradient. Let us denote the axes fixed in space 
Cfixed CSPASD by Cx.y.zD, and the axes fixed in the molecule 
Coscillating Cf’ASD by Cx*,y’, 2 *>. For small rotations & 

H y 

a±}out the three axes x* , y* aiid z' , respectively, it may then 
be shown that the components of the field gradient tensor in the 
two coordinate systems are related by the following equations 
C8, 103,1063 : 


V = c±-e^ , , 

KK y 2 : XX z yy y zz 


V = , , 

yy r.* X X* 


X z y y 


+ ^ .V , . 

X z z 


ztz y x*x x* yy 


X y T X 


V = V 

XY YX 


-8 .V , , + .8 .-*-8 , 5 V . , 

z XX X y 2 y y 


X y XX 


25y y z xx x yy y x* ^x 


2SK KSI5 y*x*x* x*x*y*y* x*z* y* 

4 


CII . 8> 
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eq eq 

V , , = - "o 

X ' X ’ d: y y 2 

^ eq* the relations in CII.8> can be rewi'-itten as; 


Substituting V ^ V » and 


eq 

XX d y 2 * y’ 


eq 

V = — f 7lC^*, +ad^. -ID - Cl-3^*,D3 

vy S X z X 


eq 

V = ^ [ -7)C.9^. -<9^.D - C 3©*. +3©^, -2D 3 

ZZ £ x y* X y 


eq 

V = V = — r -»C2© .W-© .© .D - 3© .© .3 . 

XV VX 2 z X y X y 


eq 

VZ ZV 2 y x X* x’ y'z’-’ 


eq 

y ^ y „ [ - 77 C© .-© .© .D 3C© .+© ,© ,D3 

ZX XZ 2*- 'y*x*z* y*x’z‘^ 


CII.9D 

Eqns II. 9 lead to the following form for the single-particle 
torsional oscillation Hamiltonian in the fixed QP^AS: 


ffe = A CI*-I*D + B CI*-I*D + C CI^-I*D 

T TXY TVZ TZX 


+ D Cl I +I I D + E Cl I +I I D + F Cl I +I I D , 

TXV YX TYZZY TZXXZ 


CII. tOD 
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where 






O) 

_ [l-Cn+ 3 >©^,] 


<D ^ 

S * 6“ ti-xn-35e^.J 


CO 

Ck. ^ ^ ^ ^ ^ 

D == - — 130 ,e ,^T)Cc^» ,:)] 

T 6 ^ yc* y* X* y* z* ■* 


CO 

Ct _ _ _ _ _ _ 

E = — lace ,-o ,o ,:>’^r)ce ,o ,>] , 

T O X y X X y z 


I3t _ _ _ _ 

P ^ [ 3 ca ,^ 0,0 

T 6 y X z y X z 


V i 

of 6^*. 

V 


O = 0^ -- <i O^ > I i = X* » y* » z*. < &^ > is the thermal average 


i 


For our later discussions, the torsional Hamiltonian also needs to 
be expressed in terms of irreducible spherical tensors* The fixed 
QPAS Cx,T,5K> and the oscillating <^AS Cx%y*,zO can be related 
by the Euler angle rotations O sCo/^j^O* The torsional oscillation 
Hamiltonian in fixed Cff^AS can then also be written as: 


CO 

a , 

^ = ~7r CCD + CCD + 5>* CfDJl-T'^ 

T 0 ‘ OSS '^ 22 -22 2-2 


+ ifCSi^C CD + ^ 



•i-iVesf CfD + T}C^>* CCD + sf CfD)}-T“ 

* CSO '^20 “20 20 

-iVesf CCD + 

0-1 2-1 - 2—1 Z -1 

+{yei5>* CCD + r}CS>^ CCD + sf CCD:»t“ ] , Cll.lOai 

0~2 2—2 — 2—2 2—2 

where *s are tl s VI gner rotation matrices for Euler angle 

m’ w 

rotation D = Caftyli . They are related to the rotation 

anales described above as © , © and & . The assumption that 

X- y' Z- 

rotations are infinitesimal, has also been used in Eqn II. 10a. 

3. The Dipole Hamiltonian; SC 


Another contribution to Is tbe dipole coupling. m^is 

coupling can be of two kinds. Vhen both the nuclei being coupled 
by the int erection are of the same species, we have a 

homonuclear dipole coupling. The coupling of the resonant nucleus 
with other non-resonant nuclei in the sample is considered as 
heteronuclear dipole coupling. This Hamiltonian is a two-particle 
Hamiltonian. We describe here the explicit form for both the 
abovementioned cases ; 

Expressed in its own principal axis system CEtf^AS> , the dipole 
coupling between i and j nuclei is given by [i63 ; 

SC = o>. . Cl/'Vfe^CEI .1 .-I .1 .-I .1 

» vj z't Z-J y'\ y-% X ' x x - j 


5 ; 
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Here cu = -Er .■hyC3-''£Dr7^ . r. . is the distance between nuclei 

i-i l J ‘■J '■i 

i. and i- y. and are the magnetog 2 /ric ratios of the i. and j 

nuclei, respectively. The principal axis system is defined such 

that the z — axis is along r.. , the internuclear vector. T** in 

VJ zo 

Eqn Il.iia is the rank two irreducible spherical tensor operator. 

A point to be noted at this stage is that when horaonuclear 

dipole coupling is significant in the sample under study, we have 

to take into account the two-particle nature of the full 

Hamiltonian. The total Hamiltonian is now a sum of 

single-particle Hamiltonians pertaining to each nucleus; while for 

the dipole term, the sum will be over each distinct pair of 

nuclei. In all our calculations. where homonuciear dipolar 

interaction is considered to be present, this summation will be 

implied, even when not mentioned explicitly. In the particular 

case of the quadrupole Hamiltonian we will assume that the QPAS‘s 

of the two nuclear spins coincide. Also = w = a> ; n. =n = n. 

^ ai a.\ a t j 


II. B. The Description of the Various Terms in the Hamiltonian 
in the same Principal Axis System C<1PAS3 
In the present study, the magnitude of the quadrupole 
Hamiltonian will be the largest. It is, therefore^ convenient to 
describe the system when all the terms in the total Hamiltonian 
are written in C^AS. The efg 1 nhomogenei ty <;EqnII.4D and the 
torsional oscillation CEqn II.IOIJ Hamiltonians are already in 
<2^ AS. In this section we transform X and X also to the <^AS, 

rf o 
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The rf Hamiltonian , i n <yAS 

The S)C has been defined in Eqn II. 3 along the y-axls of 
LABPAS. The relative orientation of this y-axis and the QPAS 
depends on the crystal orientation. In a pxDwder sample, 
crystallites yd.ll be randomly oriented in all possible directions. 
In general, we can represent the relation between y-axis of LABPAS 
and the axes of the CPAS by two angles C0 , 4> ^ as shown in Figure 
III. Cl earl y* 

y = xsin^ sin& + vcos^ + sccost^ sin^ , CII.12> 

JL IL. L Lm L 

and al so 

I = sln0 sin^ I 4- cos^ I cos0 sin^ I . CII.132> 

y L L UK JLi y JLi Xji Sj 

Thus, in QPAS we get from Eqn II. 3: 

= “CO C t])cosCa>t'^<i>Z)Csin<^ sin<9 I + cos© I + cosc^ sin© I D . 

rf 4 Iw L X L Y L t* ZS 

CII. 143 

The Dipole Hami 1 toni an 96^ , i n QPAS 

In this case also, we shall assume for generality that the 
E>PAS and QPAS are randomly oriented with respect to each other. 

We can then express in the QPAS by using Cin II.llaD the 

appropriate Wi gner matrix for that transformation t 16, 21, 22] 
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Figure II. 1; Relation between LABPAS y-ax’S and CPAS Cx.y.tD 
of one of the crystallites in the sample, n is 
the projection of y on the x-z plane. O is the 
angle between y and y. ip is the angle between 

Hi 

n and z. Projection of v on y is therefore 
vcos© . Projection of z on n is zcos0 . Projec- 
tion of ZCOS0 on y is zcost^ cosCny'B-O ^ = 

I-* Li !.« 

xcos<p sin© . Similarly, projection of x on y 

is xsin0 cosC tt.^' 2-© ) = xsin<^ sin© . Therefore, 

l l l 

y = xsin^ sin© + ycos© + z casd> sin© . 

)L L. L JU L 
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X - S C-lD*" 

a> o-w 2m 

I j m =“2 


Cl I. is:) 


yc3^:)sln^(9 ■{T'^ e>cpC-£i<;^ !) + T*‘_ expC+2i^ ? 

22 2“2 


yC3/8:>sinSd {T^ expC-i^ 5 - T® expCi0 !)> 

D* 21 ^ ^0 2“1 ^ ' 


+ Cl/^£5C3cos*d -1:)T“ 

D 20 


CII. 16> 


yc3/e:>sin O cosS^ Cl I - I I 5 

D 0 Xv Xj VI vj 


3cos © -1 1 

D 


~&7e 


C£I . I .-I .1 .-I .1 > 

Zv Zj Xv Xj Yv Yj 


+ yC3^:)-{sin 6^ sin20 Cl .1 +11:? 

D D Xv Yj Vv Xj 


+ sin2<9 sin^ Cl . I + I .1 5 

D Yv Zj Zv Yj 


+ sin2^ COS0 Cl .1 +I I 2>)- 

0 ^0 Zv Xj Xi. Zj 


cii. i7:> 


From Eqn 11.17, after rearrangement , we obtain 


Se = oi. .{& 1 ,I .+ al .1 + & 1 .1 } 

0 vj‘ X Xv XJ y Yv Yj Z Zv Zj' 


bCI .1 .+ I .1 .> + cCI .1 .+ I .1 .:> + dCI .1 .+ I .1 

XV Yj YV Zj Xv Zj Zv Xj YV ZJ ZV Yj ' 


CII. i8:> 
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where 

a == yC3/^I>sin^8 cosSi;^ —yCl /^>C3cos^8 — 12> » 

X 0 I> 0 

a = -yc3/a:>sin*^ cos£<^ -yci^DCScos*© -1> . 

Y 0 0 0 


a = ycS^^])C3cos^6^ —1!) 

ZJ 0 

b = ycS/^Dsin^^ sinS^ 

0 0 

c = yC3/^2)sin88 cos<^ 

r\ ^ rk 


d = yC3/^I)sinS8 sin0 

0 0 


Eqn 11.18 is the dipole Hamiltonian in OF* AS. 


II. C* The Cartesian Single Transition Operator CSTO^ Formalism 

Wokaun and Ernst [97,54] and Vega [98] have developed an 
operator formalism which has been found to be very useful in 
describing spin 3^^ NQR. We shall presently express the 

Hamiltonian terms described in Sections II, A and II. B in terms of 
STO*s. In this section we describe their definitions, 

commutation and other propenties. We mention in passing that 
STO*s are generators of SUC4> Lie Algebra [05,96]. Its use in NQR 
is not new. A very similar operator formalism constructed out of 
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Pauli matrices was used by Leppelmeier and Hahn [107] to study 
spin-lattice relaxation phenomenon in spin 3/^ NQjR. Also, the 
operators used by Singh and Arn\strong [94] and Ramachandran and 
Mur thy [1083 are very closely related to STO*s. 

The single transition operators *and 

* refer to the transition between two energy levels |r> and 
ls>. All other energy levels are disregarded* and the subsystem is 
treated as a fictitious two— level system. The single transition 
operators associated with the transition between jr> and ls>, 
which may represent a zero-* single-* or multiple-quantum 
transition are defined in the following way: 

< > = i/iac6 6 + 6. 6 

K i ^ i» jr 

< > = iy^C-6. 6. * S 6 

V xr js jr 

< > = 1^6. 6. ~ 6. S :>. CII-i93 

tzs \r jr ts Js 

The matrix representations of these three operators CFigure II. SZ)* 
lead us after elimination of rows and columns containing only 
zeroes, to the Pauli matrices. 

It is easily seen that reversing the order of Indices leads 
to the relations 

- s *”r - r -s 



V 
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I®"*" = -I*""® . CII.S05 

z z 

These sign changes have t^o be taken into account when 
constructing single transition operators for arbitrary eigenstates 
I r> and 1 s> . 

The three operators belonging to a paj-t^cular transition 
between |r> and ls> obey the standard commutation relations for 
Paul i spi n n\at r i ces : 

= il^'"® . CII.£i:> 

•■a ft Y 

where a., ft, y ~ x, v, z and their cylcic permutations. 

For operators describing two connected transitions that 
involve three distinct states lr> , |s> and |t> , the following 

commutation rules hold : 


X X 


r-i js t 

*■ y V ■’ 


= cL^:)i 


I®-"] = 0 

25 25 


j-jr-t t 

X Y ^ 




X 


r-t ja-t 
^ X 25 




r-t 


j, jr-t jS ^ ^ 

V Z ^ X 


CII. ££> 
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Figixre 11*2: Matrix representation of the three single 

transition operators ^ and They 

yc hr 525 

are represented in the eigenbasis of the correspond- 
ing Hamiltonian C54}. 


CE'MS-^ ■ 'RARv 

4CC. ^O, ' * 



r-s 
















54 


I 


r~« 


I 


r-B 


K 


•f il 


r-m 

y 


1 r>< s I , 


r-» t*><r| = Cl*" “i)'*' . CII.e7:> 

X V » * 


Because the y component is reversed In sign if we permute the 
indices, we have formally : ^ 


= I®'" . cii.ae:) 

- -t- 

However it is desirable to use ordered Indices such that m > m * 

r « 

to ensure that the raising and lowering operators increase and 
decrease the magnetic quantum numbers respectively : 

= lr> . 

= is> . CII.aQD 


II. D. The Diagonalized 96 Representation of the Hamiltonian 

o 

Terms 

In order to follow conveniently the evolution of the assembly 
of quadrupolar nuclei in the single transition operator formalism, 
we diagonalize 96 ^^ the main , term in the overall Hamiltonian. 
Also, all other Hamiltonian terms are expressed in this 
representation. This representation is denoted throughout by a 
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superscript prime C * > . In general, this di agonal izati on does not 
represent a pure rotation. However, since the quadrupole 

Hamiltonian is Her mi ti an, it can be di agonal i zed by a unitary 
transformation and such a transformation need not always 
represent a rotation. In such a case we can no longer assign 
probability weights P to states }»n> with respect to any choice of 
the z-axis in space, but only to a different set of orthonormal 
states which would, in general^ be linear superpositions of states 
characterized by magnetic quantum nun±>er [1093. 

The unitary transformation that diagonalizes the quadrupole 
Hamiltonian is 


R 


cos^'/''S O sln?^/^ 0 


O cosf/^ O -sinf/^ 


-sin?/2 O cos?/S O 


O sin^/S O cos?/^ 


CII. 30> 

where tan^ = . Also R * is the transpose of Eqn 11.30. 
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The exponential operator form of this unitary 
transformation Is 


R = exp [ *:>K} 


CII. 30a5 


where tan? = rjyys. It is interesting to point out here that since 
Cfor spin 3/^2) has degenerate eigenvalues, the unitary matrix R 
that diagonalizes it is not unique [1103. For instance, in place 
of Eqn II.S9a we can construct another unitary transformation that 
achieves the aforesaid di agonal izati on; 



cos^xS O sinfxS O 


O -cosif.xa O sin?/-'£ 


sin^xS O -cos?.'''£ O 


O sin^fy''£ O cosfx£ 


CII . 30b:) 


where tanf = i^xys. 

When we apply a Zeeman field, the degeneracy of the energy 
levels is lifted. The unitary transformation that will now 
diagonalize the Hamiltonian will be unique. 
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Using Eqn 11.30 in Eqn II. i we obtain 

ge' ^ ge ^ ge\ ^ ge\ , cii. 3 i:> 

a r f V n t 

* - 1 

where* for example* ^ ^ represents R86 R . We give here the form 
of various terms in the right hand side of Eqn 11 . 3 i in terms of 
the single transition operators: 


d(' = KCl^ = KL 

a Z 25 

L = ci 

25 25 

vhere K = Cl + 


cii. -aaD 


C 1 1 . 32 a:> 


Also from Eqn II. 4 , the ef g inhomogeneity Hamiltonian is 


1} X = D KCl‘ ^-I® *J = D KL 

lO O Ci> Z Z CO 


From Eqn 11.14 we obtain 


CII. 33 :> 


Sc = -CO CtOcosCcot+01) [ sin0 sing!> -{-C -/Scos^" +sinif 5CI * ^+1® "*> 

r S 4 L* Li XX 

+c I ‘ ■*+! ^“®:> +c yssi n? -cost" :> c I ^ "®:> y 

K H ^ ^ X X 

+cos(9 -{C si n?’ - VScosf 5 C I ‘ ~^+I ® '“l -C I * “*-1 *"®5 

L* YV YY 
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+C V3si nK +COS? D C I ‘ “‘‘+1 i 

+sin<9 cos^ ■{Cl‘~‘‘-I*“b+ScosfCl‘"‘‘+I*"®:»+2sln?CI*"®-I*"‘‘:)}3 ; 

JL !Z!Z ^ 52 55 ^XK -*^ 

CII.343 

From Eqn 11,10 we obtain 

T T5KZ; TXK TYV TYY 

» cii.35:i 

T X X 

where 

A* = {8-/3A sinf-ySB C y3cos?^si nf > Ct^cos? ~si nf 

T T T T 

= 3-ay3A^cos?-y3B^cy3sinf-cos?f:>+y3C^cy3sin?-^cos?>}-. 

c‘ = Sy3D 

T T 

d' = aVsE 

T T 

E* = aySF ; 

T T 

and from Eqn 11.18 we obtain 

» 

se' = to. [a <-cy3cos?+sinf:)ci*‘*+i®"‘‘:>-^ci‘7*+i*''®:> 

t> t.y X ' ’ Xt XV Xv Xv 
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+C -r^si nf -cos? 5CI *■*-1 * . ®:>> 

Xi Kl 

x<-cy^cos?-t-sin?:)Ci *■‘‘+1*"'^:) 

XjXji XjXj 

+Cy3si n?-cos?5CI * .*+1* 

xj xj 

+a <c-y3cos?+sin?>ci ‘~*+i®~*:)-cr*7‘‘-i^7®5 

V Vv Vi Yl V\- 

+Ci^si n? +COS? > CI ‘ " *-I * ■ ^ > > 

x< C - yScos ?+sin?:)Cl‘7*'+ 

Y^l Yj Y^ Yj 

+C ■/3sin?+cos?:>Cl‘7‘‘+l’'7*5> 

vj Vj 

<Cl^-*-I^~^y * acos?c;i‘-*-^I='~®3+SsineCl‘“®-I^'‘‘>> 

SC Zi 251. SCv set Xv Xt 

+ ecos?CI^“*+I*“®)-*-2sin?CI*~®-I*7*>> 

zj zj ^ zj zj xj xa 


+bC'C-CyScos?+sln?5Cl‘.*+I® *>+Cl‘ . *+1*7^3 

Xt Xt Xt Xt 

+<: y 3 sin?-cos?:)C I ‘ 7 *-! 

’ Xt Xt 

x< c -yscos? +si rv? > ci ‘7*-^ I -c r ^“*-1 

VJ VJ Ya Yj 

+C yssi n? +cosf :>ci y~*+I > 
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■K C -yscos? +si n? > Cl ‘ +1^"* 5 -C f~* -f"® 5 

^ ^ Yi Yl Yl Yl 

xC -C yscosf +si n? > C I ‘ “^+1 ® “ *:> +C I * ~ •♦•I * ■ ® > 

XjXj Xj Xj 

+c T^si nf -cos?:' > Cl ^ T^+i * T®:> > 5 

Xj XJ 

+cC< -cyScos;^+siu?:>ci*"^+i®~*:)+ci‘7*+i*~®:) 

Xi Xt Xi Xt 

+c yssi n?: -cosc :) c I ‘ 7*-i ^7®:) > 

^ Xv XL 

+ acos?:ci‘~*+i*7®5-+-2sirLfci*7®-i^~*>> 

SEJ 55J ^ 25j 25i Xj Xj 

+<ci*~*-i*7®:) + acos?:ci‘7‘‘+i*~®>+asin?fci*"®-i^"*>> 

Zl 251 ^ 25i 251 Xl Xl 

xC-cy3cos?:+sin?:>CI*7*+I®7‘‘o +CI ‘7*+! ^7®:> 

XJ XJ XJ X] 

+cvssin?-cos?:jci ‘7*-i^"®:)>) 

XJ XJ 

+d(c c -yscosf +s 1 nf 5 c I * 7^+i‘® 7‘‘:> -c I ‘ 7 * -I * " ® 5 

Yi Vt VV Yt 

+C ySsi n? +cosf 3C I ‘■“‘+1 J'®::) > 
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+ ScosfCI* 

zj zj ^ Zi zj ^ xj xj 


+<ci‘~*-i^.®:> + scosf c I ‘.*+i*.®:)+ssineci*7®-i *.'*:>> 

5Zi tZx. Zt 2Sv " XL XL 


x-CC-yScosf +sinf:>CI*~*+ 

Vj Y) Yj vj 

H-C ySsi nf +eos? :> CI * ^*+1 * > >] 

Cl I. 36 :> 


II. E. The Quadrupoie Interaction Frame CQIFI and Effective 
Hamiltonian Description of the Spin System 
The Interaction representation provides a relatively simple 
and convenient means for calculating the pulse response of an 
ensemble of nuclear spins experiencing electric quadrupole 
coupling [£1 ,S£,107, 111 3 . The procedure vsrhich we adopt here is 
completely analogous to the rotating frame concept in NMR, 

where the laboratory frame Hamiltonian consisting of a Biain 
ti me -independent part and a much smaller but time-dependent 
perturbation Cthe coupling of the magnetic moments of the spins to 
the rf magnetic fieldD.is transformed in the interaction picture 
into a new Hamiltonian consisting of a low-field 

ti me -Independent part PC and tiij^ -dependent parts oscillating at 
frequencies <j> and Sw. Since these oscillating terms are far from 
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resonance In the interaction picture, they exert a negligible 
torque on the spin system, the evolution of which may then, to an 
excellent approximation, be described by alone. This procedure 
of dropping the oscillating terms in the interaction 

representation is completely equivalent to the standard 
first-order perturbation theory used to tackle ti n\e-dependent 
quantum mechanical problems til S3. 

In the quadrupolar case, the transformation involved in the 
interaction representation does not correspond to a rotation and 
has no simple geometrical interpretation. Nevertheless, In what 
follows, we will see how the interaction representation can be 
used in our case in a way completely analogous to its better knovm 
use in NMR £3,4,1133. Rather than dealing with the quadrupole 
interaction frame transformation specifically. It is useful to 
describe the general quantum mechanical procedure for a change of 

9 >€ 

reference frames Eli 43. 

A change of reference frames is defined by a unitary 
transforn\ation U^^CtD, which may be time -dependent. If jy'CtI>> is 
the state of a system at time t in the original frame, the state 
at time t in the new frame is |v^tZ>>: 

jyCt:)> . CII.37> 

% 

Given that St is the Hamiltonian in the original frame, the 

We are not considering boosts here. 
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Hamiltonin ^ in the new frame may be found by the following 

-V 

argument. In their respective frames, jv'tt5> and jv' Ct>> evolve 
according to the Schrddinger equation : 

'EiT ^ jv<t5> . CII.38> 


jvk:t?> = |v<t>> 


dl. 3 ©> 


Substituting Eqn 11.37 into Eqn 11.39: 

dU 4 .. 

i — jV'<t>> + iU ^ jv^'t5> = ^ jv'<t:>> 


dt 


CII . 40'J 


Using Eqn 11.38 : 

dU 


[iC |v^Ct)> = ^ j^t:)> 




CII. 4i:> 


Eqn 1 1 . 41 implies that 


dU 

se = iC-qi-^*') u-‘ u se u'* 

dt IF IF IF 


CII. 4a:) 


•V 

This expression for ^ holds even if ^ is time dependent. 

The specific example of the quadrupole interaction frame; 
CQIFD trasf or mation employs a unitary transformation given by ; 

% "* I 

U Ct> = exp CiSif’t> . 

o o 


CII . 431) 
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where 

9?’ = af’ + Lse' =ci+A> a? 

o Q o Q 

= C 1 +A3KC I I ““*3 

ss ss 

= k‘ci‘-^^ !*-♦:> 

as z 

= K ’ L ; CK*= CI + ADK:) . CII.44> 

is the quadrupole Hamiltonian plus the resonance offset. To 
distinguish the special on-resonance case A = O, we will denote 
the coefficient by K° . 

Applying Eqn II. 4£ to Eqn 11.31 : 

% 

d expCiaf t> , 

^ = i-{ 3 T: expC-iPf t5 

dt o 

+ expCi«f ' X exp 

o O r f 1 nl o 


= -sf* +u sr*u*+u a?’u‘ + 

OoOo orfo 


u se* u~ 

o I n t o 


wher e 


=-A a*' + ^ ^ ' 

a r f i n t 


CII. 45 :> 



expCiaf’t:> a?' 

o Oi 


exp c-ise t:>= a? 

o a 


CII. 46> 
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because 

i ge* , 96^ 1 = O . CII.475 

o <k 

To evaluate the expressions of the kind in Eqn II.46» we note 
that the operators involved satisfy the cyclic conunutation 
relations in the sense that^if 
rP,Q} = LhR 
CQ^R} = inP 
and [ R P 3 = i nQ > 

with ?e = 1 ♦ or -1^ and with P,Q»R representing the 

single transition operators or their linear combinations, we get 
the following transfer Jiiat ion relations : 

expC -i^P2)OexpCi©P3) = Qpa^nO Rsinae© 

expC -i ^PI> RexpC i <9PZ> = Rcos;^^ - Qsin?€^ . CII-493 

Also, of course, if 
fP,Q3 = 0, then 

expC-iOP3QexpCi6^P> = Q . CII.50> 

Using Eqns 11.49 and 11.50 to evaluate 96 and , we 

ri i Tit 

arrive at the following expressions for the terms in the 
Hamiltonian in QIF: ^ 

== -CO C tZ>cosCa>t"»-(;t') [ si sinep -f-CyScos^* -^sin^O 

rit L, L. 
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X X r V 

4-c;i *"'*+1 +c yssi nf -cosf :>ci *~*-i 

+aose ■{Csin?-yScosf:)(Cl‘~*+I^“*:»cosK‘t + 

z-« Y y XX 

-c I ‘ ~ *-i ^ +c yssi n^: +cos^ >c i 

V V Y Y 

+sin© cos<p 

L L Z 55 55 5 B 

+2sinfCtl‘~®-I^'*:>cosK‘t-CI*~'‘+ I t»] , 

XX y y 

cii. 5 i:> 

x' = D «?* + ^' + . CII. 525 

i nt Ci> a T D 

Hie expression for each of the Hamiltonian terms in Eqn 
1 1 . 52 i s as f ol 1 ovirs : 

D = D Se' = D K CI*"*-I®"*> , CII. 53 > 

a> €fe a Ci> z 55 

= a' 

T T 55 *25 

+b’ CCI‘~®+I^"*:jcosK:’ t-Cl‘"^-I^~*)sinK’ t> 

T K K y y 

i 

-C* <CI*~®+I^~‘‘:>cosK* t+Cl‘“®-I*~‘‘>slnK* t> 

T y y XX 
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-D 

T 


» 


“E 


T 


Y y XX 

XX Y y 


CII. 54> 


The dipole Hamiltonian will take different forms in QIF 

depending upon the kind of the nuclei that are dipole -coupled to 

each other. We have two distinct cases, namely* the 

heteronuclear dipole and the homonuciear dipole. We give below 

the expressions for both of them in QIF. 

Let us denote the dipole-coupled nuclei as I and j. Also, 

let i be the non- resonant nucleus and j the resonant one. The 

transformation to QIF therefore, is carried out with respect to 
* 

of jth nucleus, i . e. , 


= expCi^* .t5 Se' exp C-iSf' . t:> . CII.55> 

Bchetero) o j n oj 


Using Eqn 11.36 for the e>qDression of X and Eqns II -4S and 

-V. > 

11.50, we write the following expression for SC : 

D^h'©' t © r o > 


X 


SX'h&t © r o > 


a>. .fa I . {-Cy3cos?->-sinf:» 

tj! X Xl 


xCC I * T^+I ® .*> cosK ’ t-C I ‘ * -I ® ~ ‘‘i si nK ’ t) 

Xj Xj Yj YJ 


+c I ^ T'‘+i ^ 7^5 +c yssi n? -cosr :> c i * T*-! * f 

XJ XJ Kj XJ 
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+a I -(Csi n^—/3cos?Z)CCI^.^+I 
y Yv. Yj 


^T*>cosK*t 

vj 


xj xj 


Vj Yj Yj Yj 


+a I .-{Cl‘.*- 

2 Zt ‘ 2 j 


I* .®:>+acos?ci‘ 

52^ Z) Zj 


■^Ssi n? CC I ‘ . "“-I ^T*:> cosf * t -CI ^ ~ * +I *• T “ 3 si nK ’ t)}- 

XJ XJ YJ VJ 


+bEI . -{Csinc-y3cosf:>C<^l‘ .^+I®"*:>cosK’t + Cl* ^-I®T'*>sinK*t) 

Ki YJ Yj Kj XJ 


t-C-»^8sin?+cose:>CI* 

vj YJ ' VJ YJ 

+ I . ■{-Ct^8cosf+sine':>(CI*T*+i^7'‘5cosK*t-Cl‘~*-I®~*:)sinK:’i> 

Yi ^ ^ XJ XJ VJ YJ 


+C I * T*-+I ^ +c ii^8si n? -cos ? > C I * T"*-! ^ T®:> }• 3 

xj xj ^ Xj Xj 


-•■ctl .^CI* *-I^ .^-J+acost:CI^ *+i^ .h 

Xt' Zj Zj Zj Zj 

+asi nf CC I ^ ^ cos K ’ t -C I ^ * +I * ■ ® :> si nK ’ t ))- 

Xj Xj YJ Yj 

+ I . ■{~cy3cosif+sin?-J(CI*~^+I®"*5cosK’ t-CI*7*-I®7*:>sinK* tj 

Zt ' ^ ^ XJ XJ Yj Yj 


+c I * T*+i +c -i<^as i nf -cos? > ci * T*-! * ~^y > 3 

XJ Xj Xj Xj 


+d [ I .-{Cl* ■*-! * T®:> +ecos?c i * ■*+! ^ T^> 

Yt‘ Zj Zj ^ Zj Zj 
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+asin?f(Cl‘”^-I*T‘‘>cosK t-CI^~*+l‘"®:)sinK t>} 

xj xj yj VI 


+I . •{Csin^-/3cos?>C<=I* .^+I®T*>cosK t + Cl .‘‘>sinK t) 

ssi yjyi xjxi 


-C I * ~ * -I ^ ® 5 +C VSsi n? +cosf :> C I ‘ T*+I }■ 3 ] 
vi Yl " yj yj * 


CII. 565 

In the homonuclear dipole case the transf orniation to QIF is 

» 

carried out with respect to SC of both ith and ith nuclei, i.e. , 

o 


Se . =expC“i^ . t5expCi^«f .tySC expC-iSC t5expC~i^ . t> 

IXhOTOO) Oi OjD OJ Ol 


CII. 575 


The expression for ^ 

0<howo > 


obtai ned f i nal 1 y is: 


St 

0{>»omo > 


w..[a {-cyScos^+sin?:>CCl‘.^+I^.*>cosK t-C I ‘ ® . *:>sinK t) 

iji** K* ^ ^ Xv Xu Yu yu 


+C I ‘ 7*+I * . +C Y3si nf -cos? 5 C I ‘ 7*-I ^ 

Kl XU % ► 


x-{-C yscos? +si n? :> CC I * . *-»-I ® '*> cosK * t -C I ‘ * -I ® si nK * t) 

xj xj yj Yj 




+c I ‘ 7'*+! * 7®:> +c ySsi nf -cos? > c I ‘ 7*-i ^ 7®:> ^ 

XJ XJ . V. XJ Xi 
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+a ^Csint-/Scosf5CCl‘7^+I^ *:>cosK t+CI ‘7^-I^7'‘5sinK t) 

y* ^ ^ Yi. Yl Xi Xt 


-c I ‘ ~ *-1 7 +c yssi nt; +cosf :> ci ‘7*+i * . *:> j- 

Yx YX Yi. Yx * 


x-{Csinf-y3cosf 2>(CI * .*+1* .^^cosK t+CI *7*-I^~*2>sirjK t) 

YJ vj xj xj ■' 


-C I * ~ * -I * ■ ® :5 +C -/3s i n? +cos?f > C I * 7* +I ^ 7®> )■ 

Yj Yj ^ ’ Yj yj ^ 


^ - 1-4 ^ 2-3 - .^ 1 - 4 ^ T 2 - 3 ^ 

-♦“a -{Cl . ~I . J -^-iBcos^^Cl -^*1 2) 

!K 5 Si asx ^ zi zl 


•<-£si nf CC I ‘ "^-I ^ “*3 cosK * t -Cl ^ * +I ^ ® 5 si nK ' t;)}- 

X\ XV Yi Yi 


..^ 1-4 - 2 - 3 . ...^ 1-4 « 2 - 3 . 

x^ci . -I . .^-^acos^ci . J 

2.1 SSj ^ Zj Zj 


+asi CC I ^ ""*-1 ^ . ‘‘i cosK * t -C I * ■ * +I ‘ ~ ^ 5 s i nK * 1 5> 

XJ Xj Yj Yj 

-•-bt -{-C -i^cos? ->-si n^D CC I ‘ ■^+1 ^ 7*:> cosK * t -C I ‘ ~ ^ -I ^ si nK ‘ t) 

XV Xv Vi Yi 


+C I ‘ . *-i-I 7®:> -^c yssi n^f -cos? 5 C I * . *-I * }• 

XV XV . % I 


x-{Csinf-/3cos?>C<^I* .^+I®7 ‘‘:>cosK t-*-Cl‘ ^-I^~'‘:>sinK t> 

Yj Yj Xj Kj 


-ci ‘7'‘-i*"®:>+cy3sin^+cos?>ci^7‘‘+i^7^:j)- 

YJ Yj Yj Yj ' 


+{Csine-y3cosf:>C<^I*“*+I®.‘‘>cosK i-^CI * . *-I ^"'‘■jtsinK t) 

^ Yi Yi Xi Ki ^ 
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-C I ‘ ^ +C ys'sinf +cosf :> C I ‘ T*'*-! 

Y\ Yi ► S 

Ki-C yscosf +si n? :> CC I ^"*+1 ^■‘‘jcosK ’ t-CI ‘ T ^ -I ^ T *:> si nK* t) 

^ ^ ^ XJ Kj Yj Yj 

+C I ‘'*+1 +C -/Ssi n? -cos?>C I ‘"■‘-I 

xj xj ^ ^ XJ Xj 

+c C -{-C -T^cos? +si nf 'J C< I * ^ cosK ‘ t -Cl ^ * -I ® " si nK ’ t) 

^ Xt Xt Yi Yt 


+C I * “*+I ^ . b +C -j^si nf -cos/T > C I * ■■‘-I ^ }■ 

XV Xv . s 


x-{Cl‘~‘‘-I^''b+8cosfCl‘ .*+1^ b 

Zj Zi ^ Zj Zj 


+8sin?C<^I ,*:> cosk’ t-CI^ .*+1* .®:->sinK*t)J- 

Xi xi Yi V.1 


+-{C I ‘ ■*-! ^ +£cosf C I ‘ "‘‘+1 

* Zv Kt ’ Zi 


a-3 

zl 




+asln?(CI ‘7^-1 ‘‘^cosK* t-CI *^“*+1 *~®5sinK* t))- 

XV Xv Yi Vi 


x-{-C-/3cosf +sinf3CCI* .^+1^ *:>cosk’ t-CI ‘ ^-I^"‘‘>si nK* t) 

xj X.l Vj vj 


+Cl‘7*+I*7®? +cy3sinf-cos?:Kl‘7‘‘-I*"^>}3 

xj xj ^ ^ xj Xj 


+dC-{Csin2^-y3cos?:)C<^l‘7^+I^.‘‘>cosK t+CI ‘“^-I®"‘‘5sinK t) 

^ Yv Yv Xi Xi 


-ci‘7*-i*. +cy3sin?+cos?:)Ci‘~‘‘+i*~®:)}- 

YV YV * , y^ y^ t 



Zj Zj Zj Zj 

+£sin?CCI^"®-I^T‘‘5cosK’ t-CI^T^+I^^T^'Jsinx’ t» • 

Xj Xj Vj yj 

^ ^acosf c I 

iSi Kt 25v SZv 

+asi nf (c I ‘ 7^-1 cosK ’ t -CI ^ ~ * fl ‘ ® 5 si nK * 

XV Xt Vv ¥t 

x-{Csi ne -y^cosei^CCI *"*+! ®7 *>cosK* t+CI t> 

yj yj xj Kj 

-Cl ‘7*-I*7®:>+cy3sinf +cosf>CI*~‘‘+I^7^'J>33 . 

Yj Yj ^ yj Yj ’ 

C 1 1 . 58> 

As pointed out earlier, in t.he QIF we can drop the terms from 

-V » ^ > 

9 ^ ^ and Sf , that oscillate at co or 2o^. TTiis leads to the 
ri tni 

following expressions for the terms in the Hamiltonian: 

= -CO Ct:>£kCl‘"^+I^~‘‘>+lCl‘~*-I®~*:>+mCl‘"^+I ®"*5 

rfi xx yy yy 

•^nCI -I >'^pCI -I -♦■I J] 

K X '^x X ^y y“* 

=-(^Ct>N , CII.593 

1 ^ 


where 
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k ■= sln^ C-»^cosif +sin^;>cosq> , 

Ik* Lj 

1 = “Cl/^Dsin^ sin^ CyScosf ■^sin?f>sl n(;^ ^ 

JL %mt 

m - “Clx;£Dcosd Ci'^Scosf ~‘sln^f2>cos<j6 » 

Imt 

CII.605 

n = Cl/^Dcosd cyScosif-siri^Dsi n0 , 

L« 

p = sin^ cos<p sin?cos0 , 

JL* JL* 

q = sin^ COS0 sljn^sin<^ * 

JL* £4 

and 


K X Y Y Y Y 

X X X X V Y 


CII. eOaZ) 
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Se’ = CO r a ■{Cl/2>cyi3cosf+sin?)’' 

0<homo> Ij X 


xi xi xi xj Y\ vv yj Yj 


xi XL X:) xj 


•+-cy3sinf-cos?><ci‘7*-^i^. ®:>ci‘7*-i^ .®:>+ci‘ . *-i*. 

XL XV. Xj XJI Xv Xt X) Xa 


+c-/3sin?-cosf:>^ci‘. * . “-I * .®:)y 

^ XL Xi Xj XJI 


+a^-{C 1 /-£:> C V3cosi!f -si n? > ^ 


x<ci‘. ^+1®, ‘‘:>ci‘"®+i®7‘‘>+fiV *-i®. .*-i®7*:>> + 

Yi Yi Yj YJ Xi Xi Xj Xj 


Yi Vi Yji Yj 


-c-/3sin?+cose-:><ci‘. *+i^ . ®>ci‘ .*-i*7^:>+ci ^ .'‘+i*7®:>> 

Yi Yi Yj Yj Yi Yi Yj Y J| 


+cy3sin^+cosfD*CI*"*+I®7®:>CI^ .*+1* .®>} 

N ' Yi Yi YJ Yj ^ 


+a -ICl . -I >CI . •-I . > 

JB* Ki ZJi SBJ Zj 


'^ScosJ^lCl -I 3CI . +I .i J)'+-Cl . *^1 . JCI . -I . >> 

Zi Zi ZJ ZJ Zi Zi Zj ZJ 


-*-4cos SCI . +I . JCI ,+I .J> 

Zc Zc Zjl zj 





+bCC 1 /S5 C 3cos*^' -si i) 


.^+1® .*5-ci*7*-i®. ‘‘:>ci *"*-1*7*5 + 

XL Xv YJ VJ VL YL XJ XJ 


Cl* *+I* *5CI* .*+1* .*5-CI*"*-I*. *5 Ci‘7*-I* .*:>> 

Yt Yt Xj Xj XL XL Yi Yj 

-CCI*"‘‘+I*"*5CI*"*-I*"*5+CI*"*-I*"*5CI*7‘‘+I‘"’"®» 

XL XL YJ Yj Vl Yl X.1 Xj 


-CY3sLrit-ac!stJ<CI -I 5Cl . -I . 5+CI -I 5 Cl . -I . 5> 

XL XL Yj yj Yl Yl Xi Xj 

+cy3sin2:+cos^5<CI*"*+I*"*5CI*7'‘+I*7*5+CI*~*+I*7^5CI *"*+I*"*5> 

XL XL Yj Yj Yl Yl Xj Xj 


+C3sin r-cos $5ccl . -I . 5CI . +I . 5+CI . +I . 5CI. . -I. . j>} 

’ ’ XL XL V j y j Yl V'l X j X j 


+cC-cy3cosC+sinf5sinC«:i*. *+I*7*5CI*7*-I^7‘*5+CI* . *-I*7‘‘5CI *7®+I^7*5 

^ XL XL Xj Xj yl Yl Yj Yj 


+Cl‘“*-I*7*5CI* .*+1* .*5+CI*7*+I*"*5CI*"*-I*7*5> 

XL XL XJ XJ YL YL YJ Yj 


+<ci‘7*+l*. *5CI*7*-I^ .*5+CI‘. *-I*7*5CI* .*+I*"*5> 

XL XL Zj Zj Zl Zl Xj Xj 


+acosi:cci *"*+!*. *5ci*7*+i* .*5+ci*. *+i*7®5Ci*7*+i* .*5> 

^ XI XV Kj !Zy Zv Zv Xj Xj 

+cy3sin^-cosf5<CI*"*-I*7*5CI*"*-I*"®5+CI*7*-I*7*5CI*"*-I*7*5> 

^ ^ Ki Xi Zj Zj Zi Zv Xj Xj 


+Scosecy 3 sinf-cos? 5 <CI* 7 *-I* 7 ® 5 CI*"*+I* 7 * 5 +CI*. *+I* 7 * 5 CI* 7 *-I^"® 5 >) 

^ ^ ^ XV Xv ZJ Zj Zv Zv Kj Xj 
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Xi Kv Yj Vj V\. Yv Xj Xj 


-<ci* ‘‘-I* . ®:)ci*"*-i^ .^5+ci*. ‘*-1*. 

Y^. Yv !ZIJ Zlj 5Cv aSv Yj Yj 


-^cosf<CI . -I . >CI . . i+CI . -^I DC I . -I . D> 

Yt Yi 5Ci ZJI m\. 2S\. Yj Yj 

YV Vv 25J (Zji 21 21 Y) Y 


+ecosfc/3sinf+cose><c'i‘ . ‘‘-•-I* . .^D+CI^T^+I^ . ®:)Cl‘ .*+1* .®:»>3 

^ Yv Yl 2j 2j Zt Zt Y). Y) 


CII.64D 

A few comments about the truncated Hamiltonians in QIF are 

in order. The process of truncation has projected out the 

ti me-i ndependent par ts from the vari ous Hand 1 toni ans that commute 

with the main term in the total Hamiltonian namely* the quadrupole 

Hami 1 toni an 9C . These truncated Hamiltonians are also known as the 
a 

secular part of the respective Hamiltonians. From Eqns 11.33 and 
11.61 it is clear that- save for the coefficients CD K in 11.33 and 

Oi 

A*^ in II.61Z) - the secular part of the torsional Hamiltonian 
has the same form as the efg inhomogeneity Hamiltonian. It is 
obvious also from the physics of the situation. Torsional 
oscillations of the units In the sample of which the quadrupolar 
nucleus Is a part, modulate the efg at the site of the nucleus. 
This results in a ti me -aver aged shift in the NQR frequency of the 



system. D K in Eqn 11.33 also is a description of distribution of 

Cl> 

N<3R: frequencies Cwith proper weightagesD around the centre 

frequency. Both can be taken into account by assigning 

appropriate values to the parameter A in our calculations. TTie 
secular part of the torsional oscillation Hamiltonian is also 
responsible for quadrupolar relaxation of the system C£,100>101] 


II. F. The Observables in the NQR Expert inent 

The signal that is observed in any NOR experiment is due to 
the voltage that is induced in the detector by the uiagnetization 

of the S 3 ^tem under study. In other words » the signal is 

proportional to the bulk magnetization of the HQR sample along the 

axis of the detector coil. We define therefore^ the laboratory 

frame observable at the radiation frequency a;- in the NOR 

experiment to be the magnetization operator given by : 

B = cos (Ia>t-+0 Z)I * CII.65r> 

i i y 

where tells us about the phase of the detected NQR signal. 

^ is an experimentally adjustable parameter and is related to the 
phase of the detector. A point to be noted here is that the phase 
<p of the pulses and also of the magnetization operator 
correspond, in NMR, to the directions in the x-y plane of the 
rotating frame, i . e. Zeeman Interaction Frame . As pointed out 
earlier, Qjp has no simple geometrical interpretation. An 
analogous assignment of directions corresponding to phases in NQR 


is therefore meaningless. 
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For an arbitrary orientation of the axis of the detector coil 
with respect to OPAS> Eqn 11.66 will take the form : 

B = cosCo>t-»"<^ >Csiri^ sin6 I •♦*cos6 I -^cos^ sin6 I !> . 

1 3 L L X JL Y L. L Z; 

CII.66:> 

Using the transformation Eqn 11.30 to transform 11.66 to the 
frame where is diagonal » we obtain : 

B ’= cosCa>t+-^ >[ 

1 ^3 XX XX 

+cySsine-cos?:>Cl‘~'‘-I^''^>>sin© sin^ 

X X Lt {.» 

+< -c yscose -si nf :> c I ‘ “ *+i ® ■ *> -ci * ~ ^ 

V Y Y Y 

+C ySsi nf +COS? 5 C I * ~ *+I ^ ~ >cos6» 

Y Y 

Z ZJ z z 

+£sinfCI*~^-I^"*j>sin© cos4> ] . CII.67> 

XX L L-* 

In QIF Eqn 11.67 will take the form : 

B *= cosCajt+9!> :>[ -{-cy!3cos?+sine:><Cl‘~^+I®"*:)cosK*t- 

i 3 X X 


cl -“I DsxnK t>-^CI > 

Y Y XX 
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+cy3sin2f-cosf >CI * sin^ 

7C Li Ii.4 

+-{-C yScos^ -si n? ><C I ‘ ~ * +I * “ *> cosK ' t +CI * ■ -I ® ■ *:) si nK ‘ t > 
-<;i‘"*-I^”®>+Cy3sin^+cosf>Cl‘~*+I*~®:>>cosd 

V y V Y lu 

* Z 5C 5S SC 

+esin2r«;i‘~^-I*~‘‘:>cosK' t-CI^“*-M*"®>sinK't>sln6> cos0 ] 

X ?c y y Li Li"* 


CII. 68> 


The part that is time dependent in Eqn 11.68 does not 

contribute to the signal we observe. Therefore* we truncate in 

the usual way, i.e. ignoring terms which oscillate rapidly Cor of 

order a> and higher!) we obtain : 
o 


B=B =k '^5+m CI^ 

13X X 3y y 3Y y 


^ ^^1-2 
-^n Cl 
a K 


K 


-.2-4.^. 1“3 .-.2--4n^ 

>*^p Cl -I Cl Z> 

^3 X X ^3 y y 



CII.695 
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where 

= -Cl/'^')sin8 sinj^ CVScosf +sin?;>cos^ 


~ -Cl^'SUsin^ sin^ cyScosf +sinf Z>sln^ 


m = -<!;i.-'‘S!>cosd cy3cosf-sin?>cos^ 

3 L* 3 


n = C 1 /''SO) c:as& C y3cosf “si n^>si 

3 L 3 


p = sin(9 cos<;^ sin^cos^ 

3 Li l-» 3 


q = sin^ cos^ sin2^sln0 

3 Li 3 


CII. 70> 


and 


N = k Cl* ^+1^ *:>+! Cl* Cl* *> 

3 3XX 3YV 3YY 

+n Cl -I I)+p Cl -I D+q Cl +I J 

3x X K X ^ay Y 

CII. 70a3 

Eqn II. 6Q Is the expression used for the calculation of the 
laboratory frame signal under various experimental situations as 
will be seen in the following. The in-phase signalCB, D will be 

i.r» 

obtained when ^ = O. Similarly when <p = the signal 

3 3 

CB > is out of phase. 

out 
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I 1.6. The Density Operator and the Observ^ations in an MQR 
exper i»ent 

The interaction frame is more than a mathematical trick. EXie 
to the usual design of a typical pulsed spectrometer^ 

experimental observations have a lot to do with the interaction 
frame. To show thls» it is first necessary to present the quantum 
mechanical description of a statistical as:;eii±>ly of spin systems 
and observables £115-1183 

The E>ensi tv O perator 

A spin system Is described by a density operator pCtZ>. If 
the system is in a pure state |^t2)>» the correspond! ng density 
operator is : 

pCt> = j^t3><v<t>j . CII.71:^ 

The system may also be and generally is» in a mixed state 
or an incoherent superposition of states. In that case it cannot 
be described by a single ket, but can be described by the 
density operator : 

pCt> ==r C Ct:>jn><mj » CII.72:) 

Where tjn>> is a complete basis of orthonormal states and is 

chosen according to the convenience that can result in tackling 

the problem at hand. pCtO is Hermitlant so that C = C 

^ nm mn 



The evolution of pCtiis governed by the Liouville-von Neumann 


equation [119,115,1163 : 


d^tD 

dt 


tCpCO,Pe 3 


CII. 73> 


The formal solution to Eqn 11.73 is ; 


pcty = uct5pco:>uct) 


CII. 74:) 


where the time - evolution operator UCt5 may generally be 
written : 

UCtD = TeiqsC-if SfefC t * !) dt * !) . CII. 755 

J o 

Using the Dyson time-ordering operator T tl£03; UCt5 is 
commonly called the propagator in magnetic resonance literature. 
The same propagator governs the evolution of pure states : 

=UCt5jv<05> CII. 765 

Given the density operator the expectation value b of an 
observable represented by Herinitian operator B, can be expressed 
as a trace ; 

b = TrCBp5 . CII. 775 

% 

The laboratory frame observable in an NQR experiment is 
typically the y-coinponent of the bulk nuclear magnetization in the 
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x-y plane of* the laboratory. The corresponding observable 

operator Is B in Eqn 11.65 with Thus the observed 

signal SCt>, is ; 

SCt5 = Tr tB U CtypCO:>V Cty^l . CII.78:) 

ou t L L 

Here U Ct-!) is the propagator in the laboratory frame. 

Eqns 11-37 and 11.43 imply that the propagator UCtZ) in the 
interaction frame is related to that in the laboratory frame by : 

veto = expC-iSr'oiKtD , CII.79:) 

L» o 

Substituting Eqn 11,70 into II.78» it can be shown that : 

SCt5 = TrtB UCt^pCO^UCO"^] coscot 

out 

-TrtB, UCt3pCO>U{:t>~‘3stncot , CII.80> 

L n 

where is obtained by putting = O in Eqn 11.60. 

In N<3R experiments the detected signal corresponds to the 
coefficients of costi'i and sina'^t in Eqn 11,80 and they are known as 
in- and out-of-phase components of the signal respectively. The 
signal components in an NQR experiment can either be in - or 
out -of -phase or both, 

% 

Experimentally* the signal SCtD is divided into two* which 
are mixed separately with two rf reference signals* one 
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proportional to coscot and the other to sincot. The tv>o results 
after mixing are then passed through separate low-pass audio 
filters. Sending the signal through low-pass filter is 

mathematically equivalent to the throwing away of high-frequency 
terms, i.e. , going from Eqn 11.68 to Eqn 11.69. This process of 
mixing and filtering has the effect of extracting signals that are 
proportional to the coefficients of cos<iit and sinoit. The two 
extracted signals are stored as the real and imaginary parts of a 
complex signal S^CtD : 

s ct:> = -TrCB^uct5pC0>Lk:t:>"‘] , cii.8i:> 

■T" “T 

with = B. •+'iB . Tlie experimental arrangement is illustrated 

in Figure II . 3. 

The above discussion illustrates the point that in NQR* 
just as in NMR» the experiments are performed in an interaction 
frame. However* there are important differences between am NMB 
and an NQP experiment. In NQR this interaction frame does not 
rotate with the magnetic field component of the external rf field 
perturbation* while in NMR the Zeeman interaction frame 

coincides with a frame rotating in the same sense as one of the 
two circularly polarized components of the rf magnetic field. In 
NQR there is no simple geon^trical relationship between the 
quadrupole interaction frame and any frame associated with the rf. 

It is pertinent to mention here that* unlike in NMR of spin 
1x2* the classical notion of precession of angular momentum 
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Figure 11*3: Schematic design of an NQR receiver that allows 

observations to be made of the signal in QIF. Signals 

(dB. near the rf receiver freqaencv a^re 

mixed i ndependentl y with rf references with a 

phase difference. TTii s mixing with and O rf 

reference is equivalent to multiplying the B and 

ou t 

B components of the si anal by cos6.^t and stna.vt* 

V n -..r « 

respectively. The complex audio signal after 

filtering corresponds to the in- and out -of -phase 
components of the signal that we observe. Passing 
through low-pass filters is mathematically equivalent 
to throwing away of the high-f requency terms in OIF 
in lowest order perturbation theory. 
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SIGNAL FROM PROBE 



IMAGINARY REAL PART 
PART 


(GOES TO DATA ACQUISITION) 


Figure II-3 
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completely breaks dovn in the case of NQR tl213. The full 
quantum mechanical solution of the problem Cl 21 3 in case of spin 
1 and 3/2 in quadrupole field does not give anything analogous to 
the classical Larmor precession. In therefore, it is, useful 

to talk only about resonant absorption of the nuclei from the 


applied external Tield 

that 

causes the 

transitions. In 

NMR of 

spin the dropping 

out 

of 

oscillating 

ter ms , when 

vi ewed 

classical ly» corresponds 

to 

the 

negi ect 

of 

the effect 

of the 


counter-rotating component of the oscillating rf magnetic field. 
Vhen not neglected, this component of the oscillating rf magnetic 
field gives rise to the Bloch-Siegert shift C1223. This, however, 
is not the case in NQR. We can not associate the oscillating 
terms in QIF to the dropping of one of the components of linearly 
oscillating rf magnetic field. As stated earlier, we can only say 
that this process corresponds to a first order perturbation 
treatment of the time -dependent quantum mechanical system. The 
selective excitation C+3/2 + 1/2 or - 3/2 *-* — 1/2I> experiment 

£123-1253 using circularly polarized rf field can be understood 
with the help of Wigner-Eckart theorem £126-1293 without invoking 
the misleading concept of precession of angular momentum in 
quadrupolar nuclei and its coupling with the circularly polarized 
rf of the correct sense. 


II. H. Pulse Sequence Propagators 
1 . General Remarks 


The effect of an rf pulse sequence on a spin system is given 
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by the propagator, defined niost generally in Eqn II. 73. Eqn 11.75 
can be viewed in several ways. The exponential operator can be 
expanded in a series ; 


UCtZ) = T 


[1 -ij* Sf Ct’>dt’-{1/^J Se Ct’Ddt*>*+. . . 3 . 


CII.8£> 


The iirae “-ordering operator Ttperforms the folioving function; 


TPKt 5«?Ct 3 

1 z 


sect :)9ect d . t 

> 

t 

12 1 



seit 'jm:t > , t 

< 

t 


Z ± i 2 


CII.835 


Use of Eqn 11.63 enables us to write Eqn II. 6S as 


uct> « 1 ~i f dt sect D 

1 ^ * 

•io 



dt se ct :>m:t . . 

2 12 


cii.e4> 


Eqn 11.64 is the Dyson series expression for a propagator 
[1603. If the Hamiltonian is piecewise constant, equal to 

...» dulling successive intervals of length t^, t^, t^, • • . » 

the propagator can be written as 


UCt> = expC*“iS^e" t 5 - . .expC-“iS^f t :>expC“i^ t > 

nn 22 41 


CII. 85> 



CII.86^ 


If ihe Hafidltonian varies continuously in time* one may 
imagine the total time interval t to be divided into many small 
subintervals during which the Hamiltonian is essentially constant. 
Then an expression like Eqn 11.85 holds to an increasingly good 
approximation as the number of subintervals increases. Thus, for 
a general time -dependent Hamiltonian : 

UCt> = lim expC~l^t:>t/^). . . 

OD 

CII. 87> 


The expression in Eqn 11.75 can be thought of as a shorthand 
f or Eqn 1 1 . 67 

The interaction frame Hamiltonian during an ideal pulse 
sequence is piecewise constant. Propagators for Hamiltonians that 
are not piecewise constant are considered in detail in Chapter IV 
and also in Chapter V, 

2- The propagator and the geometrical view of e vol utl on of the 

densi tv o perator i n a mul ti pie-pulse NOR experiment 

In NQR experiments, it is generally a good approximation to 

consider the individual quadrupolar nuclei to be isolated spins 

during a pulse, since the scalar couplings are very small compared 

to typical values of (tP. The Hamiltonian during a pulse is 

1 


therefore ; 
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^ ^ . CII. 88:> 

r f 

Meniion may be made here of the fact that when the duration 
for which the pulse is on is very long, as for example in 
zero-time resolution pulse, as explained in Chapter III, the 
internal Hamiltonians can no longer be neglected even during the 
pulse. This will be elaborated in detail in Chapter III. 

If the pulse is of duration t, the propagator for this 

duration is : 

UCt:> ■= exp [ -1 ^’^t] . CII.89> 

From Eqn II. S9 it is clear that in QIF the Hamiltonian of the 
system during the pulse is time-independent. When the pulse is 
off , i.e. , during the time interval between two pulses, the 
tl me-devel opment of the system is governed by 




int 


CII. 905 


With the corresponding propagator given by ; 


UCt5 = exp , t.] 

L rit - 


CII. 915 


The density operator at any time t^is CEqn 11.785 
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other by an orthogonal transforination. This orthogonal 

transformation Is a rotation in a Cai+1?*-1 dimensional vector 
space of colunm vectors whose components are the independent 
parameters defined above. In the particular case of spin 3/^ the 
vector space is 15-dl mensional . 

The most general state of a quantum mechanical spin s>rstem is 
best described by the well-known irreducible spherical tensor 
operators 112X^-12^2 . They provide physical insight into he 
process of creation and destruction of various order polarizations 
during the evolution of the system under various perturbations, 
both external and internal. The density operator for spin I, 
resolved into the components of irreducible spherical tensors, has 
the form 1117.130,131} 

2 1 k 

pCt5 = Cl .-- 21+1 DEI +r E T, CIDP, CtD] , CII.96D 
^ k-lci-l' 

where T CID is kth rank irreducible spherical tensor for spin I, 
kq 

and the polarizations P, CtD are in general time-dependent. Both 

kq 

T and P satisfy the Hermiticity constraint ; 
kq kq ' 


CT, )'*'=C-1D*’T, 

^ kq k-q 



c-id'^p 


k-q 


CII . 97D 


Clearly, the evolution of the system of nuclear spins I is a 
rotation of column vectors whose components are the polarizations 
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P|^^CtZ>; k=l,...2I; q=-)c , . . . , +Jc . Evolution under rf pulses Crank 1 

tensor, Eqn II.3aD conserves the rank k of P , but changes q; 

kq 

evolution under internal Hamiltonians Crank £ tensors, Eqns 
II.4a,II.10a and II.lta5 does not conserve the kth rank 
pol ar i zat i on . 

For spin 1-C£, k=l ; q=— 1,0,+1. The rotations referred to 

above are in 3-dimensions and the polarizations are all rank 1, 
i.e. the vector polarizations. The rotations of these vector 
polarizations can be visualized by using Bloch vector model C13£] 
and they represent the rotations we are familiar with in 
3-dl mensi ons . In usual magnetic resonance experiments only vector 
polarizations Crank 1 polarization!) are detectable. 

For spins I 5: 1, higher rank polarizations i.e. tensor 
polarizations, also contribute to the evolution of the system and 
a simple Bloch vector kind of geonustrical picutre is then not 
possible. For example in spin 1, rank S polarizations come into 
play and the evolution of the system is on a surface of 
8-dimensional sphere. 

For spin 3/CS nuclear ensemble, the tensor polarizations of 
rank i,£ and 3 contribute to the evolution of the system. 
Mention may be made here that in usual NQR experiments also for 
nuclei with spin I > 1 , i t is only the vector polarization Ci.e. 
rank 1 polarization!), that is detected. 

In spite of the physical insight the irreducible spherical 
tensors provide, we have not used this basis for our calculations 
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here. Instead, we have resolved the density operator for the 
mathematical convenience in single transition operator basis ; 

pcty = a 1 +E £ . CII.9©:* 

o . . t> p 

KJ p=X.V,2! 

For spin 3/S; 1* J 4. The single transition 

operators are an over -representation of the space spanned by spin 
I. Nevertheless, as will be seen in the following, it is quite 
useful for NQR spin— dynamics problems. 


II. I. The equilibrium dei\sity operator 

The pulse sequence propagator acts on the initial density 
operator, transforming it to a final density operator according to 
the Li ouvi 1 1 e- von Neumann equation 11.73. If the initial density 
operator describes a spin system at equilibrium, it is given by ; 

p =fs>cp C-ffe klD^TrCexpC-tff/'IcT)] . CII.QO> 

eq 

wiiere Sf is the Hami 1 toi\l an of the system with no rf fields 
present. Since the dominant term in PC is the quadrupole term and 
since aj^<<kT at temperatures above a few degrees K. it is a good 
approximation to write ; 

p = tl ykTJ^ 


CII. 100> 
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Here N is the total number of spin states, or the dimension of the 

system. This is the high teii 4 >erature approximation 

£3,4,17,133}. The unit operator part of p , commutes with all 

eq 

propagators and does not contribute to the observed signals. 
Therefore, it is usually dropped along with the constant of 
proportionality multiplying which only determines the absolute 

signal amplitude, leaving an initial reduced equilibriujn density 
operator of 

The equilibrium density operator being proportional to is 

a second rank: tensor CEqn II. a>. The first or the preparatory 
pulse prepares the system for observation in the x— y plane of the 
laboratory. The density operator during this time evolves in a 
5-d±mensional subspace of the full space. It is only after this 
pulse, when internal Hamiltonians come into play, that various 
order polarizations start getting mixed and for ail subsequent 
times the evolution is in the full space of spin I. 


II. J. The NOR Line Shapes 
1 . Gener al Remar k s 

The various external rf field perturbation schenuss affect the 
different internal static and dynamic interactions in a bulk 
sample differently. The overall magnetic resonance signal, 1 . e. 
the resonance line shape, depends ma.rkedly on these interactions 
and any modification of them results in a suitably modulated line 
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shape. It is Lhererore esservLial Lo evaluate very closely the 
relation between the internal interactions and the line shape. 

All contributions to a magnetic resonance line can be broadly 
classified as homogeneous or inhoinogeneous. The purest form of a 
homogeneous line is the one without separate contributions to 
different parts of the line, i.e. all the nuclei in that line 
make Identical contributions with their intrinsic linewldths the 
same as the composite. There are then no spectral shifts with 
respect to each other due to individual contributions. An 
i nhomogeneousl y broadened line, on the other hand is a sum of 
non -over 1 appl ng individual lines. The line shape in this case is 
determined by a distribution of shifts. Although there may be 
some life-time broadening connected with individual spectral 
packets* no coupling between them exists . Figure II. 4 

illustrates pictorially the ideas so far discussed [134“'1363. 

Th& following are some of the interactions in solids that 
will give rise to a homogeneous broadening of the magnetic 
resonance line shape : 

CaD Dipolar interaction between like spins 

Cb> Spin-lattice relaxation 

CcD Interaction with the radiation field 

Cd5 Motional ly narrowed fluctuations in the local field* 
etc . 

In case ther.e is a coupling between the spectral packets, the line 
is known as a heterogeneous line. 



Figure II. 4 : Schenvatic diagram showing Ca) homogeneous broadening 
and C b> i nhomogeneous br oadeni ng f 1 35 ,1363. 
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Some t.ypica.1 sources of inhomogeneous broadening include : 
CaD Hyper fine interactions 
C b!> Anisotropy broadening 

Cc2> Dipolar broadening between spins with different Larmor 
frequencies Cheteronuclear dipolar interaction^ 

Cdl) Inhomogenei ties in the electric field gradient Cefg^ in 
pure quadrupole resonance and 
Ce> The applied magnetic field in the case of NMR. 

The broadening due to dipolar interaction deserves a special 
note. Both homo- and heteronuciear dipole couplings give rise to 
a variation in the local fields at the site of the resonant 
nucleus. Thus, in principle, we should have an inhomogeneous! y 
broadened line in both the cases. However, in the case of 
homonuclear dipole interaction, the energy absorbed at one 

frequency of the dipolar broadened resonance line is quickly 
transferred to all the spins, whether or not they are in a local 
field exactly corresponding to the applied rf frequency 
[134,137,1383. This assumption is Justified because of the 
possibility of mutual spin-flips between neighboring nuclei 
brought about by the dipolar interactions. Thus in this case the 
line becomes homogeneous. In the heteronuciear dipole case, the 
mutual spin flip between neighboring nuclei becomes difficult as 
their Larmor frequencies are ‘too far apart. Therefore, the line 
r emai ns i nhomogeneous . 

In case of spin-lattice relaxation the spin-flip, and hence 
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the transfer of energy to the lattice, comes about because of the 
coupling of the spins with the appropriate lattice modes. The 
broadening is again homogeneous. Since the other interactions do 
not interest us in the present work* we shall not discuss them 
here. 

The spin 3^^ HQR Li ne Shape 

The static inhomogeneous contributions to the HQR line shape 
comes mainly from efg inhomogeneity and also to a lesser extent 
from heieronuclear dipole interactions. 

The homogeneous line-width in an HQR system has its major 
contribution from coupling of the quardupolar nuclei to the 
oscillating efg due to various modes of lattice vibrations. In 
our calculations this line- width will be described by the secular 
part of the torsional oscillation Hamiltonian, Eqn.II.61. Unlike 
in solid state NMR of protons, the homonuclear dipole Hamiltonian 
contributes to the homogeneous line- width of the line in a very 
Insignificant manner. Generally, this is because of the large 
size of the quadrupolar nuclei and their low magnetogyr i c ratios 
as compared to that of a proton. The calculations of Gutowsky 
and Woessner E1013 show that, the torsional oscillation is 10 to 

times more powerful in relaxing the quadrupolar nuclei than 
the dipole coupling. 


II. K. Stunmary 

This Chapter serves as the foundation for the ones that 
follow. After introducing all the Hamiltonians that are at work’ 
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in spin 3/^ NQR in the usual Cartesian basis of angular momentum 
operators, we have presented these Hamiltonians in a 

representation where it is most convenient to follow 
mathematically the time-development of the system. To begin with, 
the Hamiltonians are defined in their respective principal axis 
systemsCPASD . The quadrupole Hamiltonian is the main term in the 
overall Hamiltonian and, therefore, all the other Hamiltonians are 
then expressed in the Quadrupole Principal Axis SystemCCtf^ASi . To 
simplify later mathematical developments, all the Hamiltonians 


are now converted 

to 

a representation 

where 

the 

quadr upol e 

Kami 1 toni an 


is 

di agonal . The 

concept 

of 

quadr upol e 

interaction 

frame 

CQIF5 is then Introduced, 

and 

al 1 the 


Hamiltonians are transformed to QIF. In complete agreement with 
the lowest order perturbation theory, the oscillating terms are 
dropped here from all the Hamiltonian operator terms. The choice 
of the operator formalism for simpler description of the system is 
the generators of the SUC4> Lie algebra, known in magnetic 
resonance as single transition operators or fictitious spin 1/^ 
operators. This completes our task of defining the Hajiii 1 toni an of 
the system for applications to different problems in later 
Chapters. 

We also have described in this Chapter the observables in the 
experiment. The concept of the density operator and how it 
helps to follow the time -development of the NCa? system using 
Liouville-von Neumann equation has also been introduced. The 
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experimental procedure of detecting an NQR signal has been 
correlated with the calculation or the expectation value of the 
observable In NQR using the density operator. 

A detailed description of the propagators in multiple pulse 
experiments has also been included in this Chapter along with a 
geometrical picture of the evolution of density operator in the 
space spanned by the basis operators of the spin system. 

This Chapter concludes with a discussion of the nature of the 
internal Haml. 1 toni ans and how they modify the NQR line shapes. 

In Chapter III, we shall apply the formalism developed here 
to some simple problems, namely the single pulse response and the 
two pulse response. Also, the formalism will be applied to the 
zero-time resolution technique in non-axially S3^minetric HCS5 of 
spin 3/^ nuclei. 
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ONE- AND TWO- PULSE RESPONSES AND APPLICATION TO ZERO-TIME 

RESOLUTION CZTR) TECmiQUE 

In this Chapter the formalism developed in Chapter II is 
employed for some very well-knowi phenomena in NQR, namely^ the 
one-pulsi response known also as the free induction decay CFIDD 
and the two-pulse response, otherwise known as the spin echo 
phenomenon. The explicit solutions presented here have been 
obtained in the past using the standard ti me— dependent 

perturbation theory to solve the Schrddlnger equation [63 and also 
the density operator technique [1,139,1403 where the solution is 
obtained by manipulating explicit matrices. Our purpose is to 
demonstrate the fact that the diagonalized representation and 

the single transition operator forraallsra lead to a relative 
simplification of the typical NQR dynamics calculations. Also, we 
shall demonstrate theoretically, perhaps for the first time, the 
validity of zero-time resolution technique to non-axially 

symmetric NQR of spin 3/12 solids. 


III. A. A Special Property of the Quadrupole Hamiltonian for 
Spi n 3/12 , 

The quadrupole Hamiltonian in 3/i2 space satisfies 


SC 


a 


= K^l 


2 


€111.15 



vhere K 


£ 


(1 + f 


This was first observed by Pratt riil3 in a special . basis. 
Later it was noted by Ravishankar Ci£l] that this property is 
basis independent. Hence 

iSf 

exp C -iSe^tD = cos Kt - sin Kt . CIII.a:> 

so that the corresponding density operator at any time t, can also 
be written in a basis-independent laanner as 

pCt:5 = pCo> + C-4ir- sin £Kt + Sf )[pC'o5.Sl[f ] , 

^2 Ok Ok 

cm. 3> 

where pCoD is the density operator at time O, We shall be using 
this property in the following wherever required Csee the two 
pulse response. Section III.O. 


III. B. One-Pulse Response 

Let us consider the effect of a single rf pulse Cas shown 
in Figure on the equilibrium reduced density operator. We 

have made the simplifying assumption that our assen±)ly is a system 
of non-interacting nuclei. This problem may now be solved for a 
very general case of arbitrary flip angle and phase for the pulse. 
The solution in this form shall be used in forthcoming Sections 


and later Chapters. 



Figure III.l 


One-pulse response. 
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From Eqn Il.lOO, Lhe reduced equilibrium derxsity operator is 


eq 




KL 


cm . 4!) 


ITie one pulse propagator, from Eqn 11.89, is 


= exp f-iS? .t 1 , CIII.5> 

W rt W 

* 

vhere is defined in Eqn 11.59 . Using Eqn 11.74 : 

pCt > = UCt > UCt 5* 

w w a w 

= exp l-i9e\ t 3 k'l exp li^\t 1 . CIII.6D 

rf w ^ rf W 

From Eqn 11.89 * we have 

»(? i = -<x>° Nt = N 8 , CIII.7> 

rf w 1 W o 

where - & . Therefore Eqn III. 6 gives 

1 w o 

pc i D* — K exp C -i N 8 !> L expC iN8> . CIII. 8D 

^ w ^ a o 

i 

To evaluate Eqn III. 8, we use the Baker -Campbel 1 -Hausdorff 
CBCID formula [141-1511. It states that for two non-commutati ve 
operators A and B, i . e. if £A,B1 0,then 
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expC-iAtZ> B expCiAtD 


=B-Cit:)C A.BJ + 


Si 


[ A.[A.BJ1- 


Cit)' 

3i 


[ A.CA.CA.BS] J + 


cm. 93 


Using Eqn III. 9 in Eqn III. 8, we obtain : 


pet 3 = K CLcosCO i'3 
w o 


sin C8 r3>. 
r o 


cm. io> 


where 


M = -ikCl‘“^-I® *3+ilCl‘"*+I®"‘*3-HmCl‘ ^-I®~*3 

V Y XX XX 

-i nt I +I 3 -i pC I +1 3 +i qC I -I 3 , 

V Y Y Y ^ X X 

cm. 113 


and k-,l»m,n,p and q were defii-ied already in Eqn 11.60. Also, 
can be shown that 

IN. L 3 = M . ■ l-£a3 

[N. M 3 = . cm. 1853 

IM, L 3 = 


it 


N 


cm. isc3 
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where 


2 ,2. ,2. 2. 2 Z 2 

r = k •+ 1 + m + p. ■+ p -I ' q 


CIII. 12d^ 


Eqn III. 10 is the density operator in QIF after a pulse. In 
principle we can go ahead and calculate the expectation value of 
the signal using Eqns III. 10 and 11.69 in Eqn 11.77. In order to 
compa, e the result with that in the literature, we shall further 
simplify the situation. We shall assume the case of a single 
crystal where the rf is applied at resonance (IA=0Z> along the 
y-axls C<p-0-, <p^=0'3 or along the x-axis <i4> =0;^^ =rr/£!> of the 

i 

CS^AS. Therefore K = K , and 

Ca> detector along y-axis 

L = l‘-^-I^-* 

21 2 ; 

N = .l;t'sin 5' -yscos , 

Cj X * 

M = |iCstn ? -yScos . 

Also r reduces to r given by ; 

V 

r = 5 Csin ^ -VScos O . CIII.16> 

y B 

1 

Thetr 0 ^ Etcjn III. 10 in Lhis spGfCi3.1 C3.st5 r«&dicic<©s to : 

p Ct :> = K<Cl‘~*- I^'^^^cosC^ I^"*>sinCd 5> . 

'^y w Z Z Y X X V 


cm. i3> 

cm. i4> 

cm, is:) 


cm . 175 



no 


where 6 - 6 r 

V o V 

CbZ) Detector along the x-axis : 

L = . 

z z 

M = t i^~‘‘:)cyScos ^+sin if:> , 

Cji Y y 

H =- I; I^“*5C-/3cos f+sln . 

Cl 3 hC JC 

and r will reduce to given by : 

r = I; cyScos f+sin 
X i 3 


till . 17a:> 


cm. ie> 


cm. is:) 


cm. 80!) 


cm. 8i> 


Therefore 


p’ ct :>=K<ci‘ *-jcasc& :>+ I® ‘‘:)sinC©J>> » 

X w z z X y Y X 

cm. 885 

wher &&=&!'. cm. SSa5 

X o X 

* 

The corresponding magnetization operator required 
calculate the signal from Eqn 111.88 is Cfrom Eqn II. 695: 


Lo 



ill 


f+sin , CIII.83:> 

where the phase of the detector <p^ = rx.^. Therefore, 

<B > = -C§ 5Cy3cos2:+sirif>sinCd :> . {:ill.£4:> 

X Cl X 

In the laboratory frame, from Eqn 11.80; 

<B > = +c| >cy3cosf +sinf'lsinC8 :)sin Kt . Clll.as:? 

K C X 


Substituting for various parameters gives the value of the signal 
to be : 


<B > = )slnC^.C3V:>t). 


C3+}) > 


cm. a6D 


From Eqn III.S6, the laboratory frame signal for the axially 
synuiietrlc case Ct)=o> is : 

u> 

<B > = -^•/3sinCCy3-''8:>8 )sinCCc,i t) . CIII.a7> 

X >5=0 4 o Ol 

The results III. 25 and III. 26 compare with the reported 
expressions for the same £1,23. It can also easily be shown that 
there will not be any signal along y and z-axis of Cff'AS. Also, 
even along x--axis a signal is possible only when the phase of the 
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detector <p^~ The point to be noted is that phases of the 

applied rf and detector respectively differ by n/z. 

III.C. Two-Pulse Response t 

In this section we consider the effect of the following 
C Figure III. £> two pulse scheme on the density operator at 
equilibrium . If the density perator at time t > is given by 

pC tZ) , then 


pCt:>=OCt-r> ^PCt :> ‘qCt> ^pCt > V PCt >QCr^PCt ■:>C8:t-T>. 

W V w w 

ciii.ae:) 


ViThere 

PCt :>= expCi^ t > =e.xpCiN & ':> = expCiCl‘"*+I^“*:)e !> . 

W rfw o xxX 

cm. aQD 

CiCT> = expCiSK* T> = exp<iA.K:Cl‘"^-I^~*:>T) . CIII.30> 

QC t-TZ> =expCi A^KC I ^ • CIII.SiD 

In wi'iting Eqn 111.28, we assumed that both the rf pulses 
are applied along the x-axis of CiPAS of a single crystal sample, 
so that e = <P -n.'' 2 , the phases ‘ of both the pulses being <p=0. The 

Kw 

evolution in CXTF in between the two pulses and also after the 
second pulse is under the quadrupolar Hamiltonian. This is 



Figure III. 2: Twc-pulse or Spin Echo sequeiice. 
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Figure I1I-2 
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because we have assumed electric field gradient inhomogeneity in 

the crystal to see the spin echo effect. A is a factor that 

i 

incorporates the effect of efg inhomogeneity in Eqns III. 30 and 
III. 31 . 

Using the commutation relations of the single transition 
operators and Eqn. III.3, we get from Eqn III. £6 the expression 

for densitv operator at time t T-+^t as follows: 

w 


pCf!} = aCI*~^- 
1 as 


CI*' 

Z 4 V 


2 r-T*-2^T3- 

-I >+>' Cl +I 

Y IX X 




where 


cm . 32:) 


a = K Ccos^O ~sin*0 cosCA.Kr) , CIII,33a:) 

4 XXV 

ft - K [sin CAKrDSinO C A KC t-T)>-sin© -cosO Cl+cosC A.Krt:) 

' 1 i XV XX V 

X cosCA.KCt-T)>3 . CIII.33b> 


= K [sinC^ COS& Cl+cosCA. KT)DSinCA.KCt-T>2>- 


si nC A Kr3>si n© cosC A.KCt-r) 3 . CIII.33cD 

V X V 

The signal B will appear along the x-ax±s of QPAS. The 

K * 

expression also clearly shows that there can be two signals 
corresponding to two detector phases 0^ = O and ny'z. 



ilQ 


When 0=0, 

9 

^^Kd> =0 “ I . CIII.34:) 

Also» when 0 = tx/ 2 » 

3 

~ I <^'>^cos?+sin?:)r . CIII.35:> 

^3 * 

Eqns III. 34 and Il.c'i can be easiiy translated into the 
expression reported £1,23 by substituting for parameters . In 
this case also there will not be any signal along y and z axis of 
OP AS. 

III.D. TIv 5 Zero-Time Resolution Techivlque applied to NOR of 
non-axlally symmetric systeins containing spin 3x2 Nuclei 
1 - Gener al Remarks 

In solids the initial portion of a free induction decay CFIDD 
signal. that has important information, often becomes 

unobtainable, especially when strong spin-spin coupling among the 
resonant nuclei is present. This is due to the following; 

CiD The FID function FCt:>, has its temporal origin at 
approximately the centre of the rf pulse; the N<3R signal is 
not normally observable while the rf excitation pulse is on , 
so that the effective measurement time of the FCtI> is 
typically only one-half of the rf pulse length; 

Cil> more seriously, electronic causes also contribute , that is, 
the initial portion of the FID is masked due to rf transients 


arising from 
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C a? finile Q Ccjuaiily— facLor) of Ihe various resonant 
circuits excited by the rf pulse, and 

Cb> saturation of the high - gain signal amplifier circuits 
Con the rf receiver side? due to the feed through of rf 
power from the excitation pulse £93,1523. 

In solid state NMR spectroscopy, several techniques have been 
developed to overcome this problem £90,93,153-1563. Zero-time 
resolution C2rrR:2> technique £93,1563, so called because it permits 
one to resolve the signal near time t = o, has been found to be 
simpler and much superior to other methods available in the 
1 i terature. 

In hlQR of spin 3/^ with a non-zero dipole-dipole coupling 

among resonant nuclei also, fiiitte recovery time potentially can 

distort the FID. The possibility of overcoming this finite 

recovery time was discussed for the first time by Pratt £1113. 

He used the solid echo sequence of Powles and Strange £1533. 

More recently Singh and Ai'mstrong £943 have applied the ZTR 

technique to a spin 3/2 system in an axially symmetric electric 

field gradient Ci)=o5. But the scope of their work is limited 

because most of the compounds studied in NC^ spectroscopy possess 

non— axially symmetric electric field gradient {Ij}?'o3) £ 2, 157’, 1563 . 

The purpose of this section is to analyse theoretically the 

general validity of ZTR method for non-axialiy synuiietric spin 3/2 

% 

systems also. 

2. ZTR technl aue aopl i ed to NOR of non-axi al 1 v symmetric spin 
3/2 systems 

ZTR pulse is shown in the following figure ; 



iie 


Figure III. 3: The first stage of 

short initial n^z 
spin-locking pulse 
with respect to it. 


the 

ZTR 

pulse 

technique. The 

pulse 

is 

followed by a long 

that 

is 

phase 

shl f t ed by tr/'z 






lao 


As is apparent from the above figure, it is essentially a 
combination of the basic one— pulse experiment followed by or>e or 
more spin-locking pulses. The success of this method in 
recovering the shape of the free induction decay of the 
equilibrium magnetization following a single rf pulse has been 
predicted [93] to be dependent on the ability of the long rf pulse 
to prepare the spin system in the spin-locked state [159,160], 
that is, on the existence of a spin temperature in the QIF. The 
effect of the "ZTR pulse sequence on a spin 3.^2 quadrupolar system 
can now be evaluated. 

The reduced equilibrium density operator of the system under 
study is CEqn 11.1005: 

p ~ af = KL . cm . 45 

sq O 

Also, we assume a homonuclear dipole coupled nuclear 
assembly. Therefore, apart from the driving rf Hamiltonian the 
system will evolve under secular part of the dipole Hamiltonian 
during ZTR sequence. Let us assume that the pulses are applied at 
resonance along the y-axis of the laboratory frame. The rf 
Hamiltonian for the first pulse C<p=0^ is, from Eqn 11.59, 

= -a> Ct5 N ‘ , CIII.365 

rftl) i«> i 


where 



lai 


N 

1 


i X X 


IV Y 


+ p Cl 

1 X 


1-9 


X 


with 


CIII.37:> 


k - sin & sin <p C -y^cosf +sin?> , 

* £ IL» Udi 

m = (^3 cos & Csinf-yScosCO , 

A 2 Ljh 

p = sin ^ cos <p sinf 


cm . 38> 


Eqn III. 37 has been derived from Eqns 11.60 by putting 0 = 0. 

Siniilarly, the rf Hamiltonian for the second pulse C0 =tt.-''z> 
is gi ven agal n , f r om Eqn 1 1 . 59 , to be : 

dt' = -o Ct> N . cm. 395 

r f t zi H2> 2 


where 


N = k Cl 

2 i Y 


•z_i3-*5_ m + pCI^ ‘‘+1* 

Y 1 X X i V Y 


cm. 40D 


t 


Also, 


tN ,LJ = -iN 

i H 


cm. 415 



iBB 


tN .L3 = IN 
z 1 


cm. 421 ) 


[N .N 3 = ir*L 
12 1 

wtiere 




r 


2 

4 


, 2 , 2 

k + m 
1 1 



cm. 44> 


The effect of the first pulse of duration t » on the initial 

w 

density operator is : 


p Ct :> = expC-^ia^ t N >KL expC-lix^ t N :> 

W ^ 1<4) W X ^ 4C1) W 1 

cm. 45> 


N 

~ K{L cosCoN CtZ)t r — sinCo-v t r } 

* i<i> w 4 r 4a> w 4 ' 

1 

cm. 46> 


For a tT ./2 pulse 


c»3 t 
Kli w 


rr/'a * t her ef or e » 


p Ct y 
w 



cm . 47:> 


9€ 1- 

When Vi^e are using a powder sample, *n^z* pulse means a pulse that 
gives the maximum signal intensity and it corresponds to 0.66 n, 
rather than the single crystal value of 0.5 n. A good discussion 
and derivation of this is given in Petersen's Ph. D. thesis [1613. 
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If we assume the initial applied tt/'S pulse to be a 6~ 

function, then t = O. 

w 

Alter the excitation pulse, and in the absence of any further 
rf fields, the evolution of the density operator for the spin 
system is described by 

+t > ^ UCrp N UCr . CIII.4ey 

1. W X X Z dL 


where 


UCt > = exp C-i^’ ^ T ) . CIII.49:) 

i JXhonio) i 

In Eqn III. 49, ti me-i ndependent secular part 

of the homonuclear dipole Hamiltonian given in Eqn 11.5©. 

Now when the second long spin-locking rf pulse is applied, 
the operator governing the evolution of the system is : 

UCt > =expCiCSe' + ffiT' D t ) . CIII.S05 

Z rI(Zi 0 < homo) 2 

The second pulse is a soft pulse and its duration is long. 

* 

Therefore, the effect of SC , can no longer be neglected even 

DChomo) 

during the pulse being on, and hence its inclusion in Eqn III. 50. 

Also SC >>SC • therefore, it is useful to go to a second 

rft2) Bfhomo) 

interaction representation 11531 defined fay : 

U = ejqpCiSf’ t:> . CIII.5i:> 

IFO-f{Z» rf(2) 
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Label i ng 


•\* 

le“ 

I><homo> 


as 


that par t of 


X ^ which commutes 

0lTiOT>1O> 


w-ith only is effective in di'lving the spin system 

evolution In this interaction representation during the time r - 
Denoting the density operator in the interaction representation of 
the rf as p, one has 


pCr Z> = 

2 


expC 


■LX T 3pC02>expC 

D<hotnoi 2 


iX 


DvhomoJ 




cm. 525 

•V 

where p C05 = p Ct +t 5 . Cl 1 1.535 

i w 

At a time t after the spin-lockring second pulse has been 
switched off , 


pet +t + t + t 5 
w i 2 


K 


exp<-i^ ^ r)-{expC-i T^)expC-i x'^ tW 

r Dfhomo) Dthomo> 2 D<hoTno> 1 2 

1 


X expClSf ^ T )expCi S<f _ t )>expCi^*^^ .. t) 

0<homo> 1 D<hoTno> 2 D<.hoino> 


cm. 545 


For the first pulse along the y-axls of the laboratory frame, the 
magnetization operator is given by putting <p = n^z in Eqns 11.69 
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and 11.70. This can easily be seen to be nothing but N given in 
Eqn III. 40. Therefore, 

<N^Ct>> = Tr{^kt:> , CIII.5S:> 

where t in this case is given by 


t + T +T +T 
W 12 


cm, 56 :> 


If Eqn III.5S is to be identical with the FID then pCt +t +r D 

w 1 2 

must be diagonal, or else the off-diagonal elements must have no 
effects on the further evolution of <N^> after the spin-locking 
pulse, so that <N^Ct2>> evolves as though the matrix were diagonal 
at T = O. It appears that for latter seems to be the 

case, although this cannot be shown rigorously. We, therefore, 
present a semi -rigorous argument that the off-diagonal matrix 
elements will make a negligible contribution to the evolution of 
<N^Ct2>> provided that Po*' further justification we have 

to rely on experimental results. 

The useful analogy given by Vbllmers et al [933 in case of 
NMR can also be invoked here. If there exists a component of 
magnetization along the y-axis of laboratory frame that is out of 
phase by n/'z, it disappears in a time equal to several due to 


the effect of the Hamiltoqian X 


D<homo> 


The Hami 1 t oni an 


X has the same form with respect to X _ in the second 

achomo rtiZi 

^ ' 

interaction frame as fiP ^ has with respect to X in QIF, 

0fhOTOO> Q 



1S6 


except for a constant factor . Therefore, Its effect in the 
second interaction frame should be to cause the decay of any 
component ot magcnetization in a direction transverse to the 

-V * 

direction of a rate that is inverse to the constant 

♦ 

factor . Since , and X commute with L and N , 

Dihomo) Dfhomo) 2 

respectively, the magnet! zati on components* <L>and <N > would be 

2 

unaffected. After the magnetization transverse to the 

wv, 

direction of decayed, we may expect that when the 

spin-locking pulse is removed, the out-of -phase magnetization 
along the y-axis will be the only net magnetization, and will 
commence decaying again with the same time dependence as 
originally, but with the time origin at the end of the 
spin-locking pulse. 

To return to our derivation, the exponential operator in 
Eqn III. 48 can be expanded in a power series so that the density 
operator at time t , can be written as Cputting t = O in the 

1L V 

6-pulse limit:> : 


i T 


K 


pCT >= N - .N 3 

i. 2 1: D( homo > 2 


, 2 2 
i r 

1 


ai 


*■ D<homo> D(hoino 2 


cm . 57> 


N 


p<T^> 


r X Ct 5N * E X -Ct > N. L, 

,"^"^12 i j2 ,■*- 4 J4 k 


+ . . . cm . 583 


J = i 


_ 

The magnetization component corresponding to L is not directly 


detectable in NCS^ experiments. 



Eqn III. 58 follovifs from Eqn III. 57 when the commutators are 


evaluated and various terms collected. The first term in 
Eqn III. 58 is diagonal while all the rest of the terms are 
nondiagonal, becoming more complex as the order of the terms 
increases. If we have physically equivalent nuclei in the system, 
is independent of J and the first term in Eqn III. 58 reduces 
to y Ct ?N . For the case where all the nuclei are not 

Z 4 Z 

equivalent, it is advantageous to rewi'ite the first term as 


N 


r y. Ct >N. = y Ct D N 

. ^,iz 4 jz Z 4 S 




N 


r y . Ct D - y Ct :>N. ; 

^ ^iz 4 ^Z 1 JlZ 


.. N 




cm. 59:) 


cm . 60!) 


It is rigorously shown in Eqn 111.61 below that 
the dominant term in . does not affect <N > during the 

JDinomO 2 

spin-locking pulse. 


< N Ct +T +T >> = Tr CpCt +T +T D N > 

ZWiZ W4Z2 


- Tr<expC-i X , r DexpC-iiif ,, 
p ^ 0 < homo ) 2 0^howo> 1 3 


X expC i r >exp Ci 

^ 0<homo> 5L 0<nomo> 2 2 
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K -v * - ^ » _ 

- Tr<expC-i SC r 5N expC+i^l" ^ t >N 

^ ^ 0<homo> 12 D < homo) i 2 


< N C t ■*-T 2> > 

z w 1 


Ct :»Tr<N^> ; 

2 1 2 * 


cm. 6i:> 


in the above the first exponential has been cycled to the end of 

the product, where it conunutes with N and cancels with its 

2 

complex conjugate. Thus. regardless of r^, the magnitude of 
< N^> is the same at the end of the spin-locking pulse as at the 

•v# * 

beginning Cthose terms in SC . which produce spin-lattice 

rela.xation in the second interaction frame, and thus 

somewhat attenuate < N >, have been truncatedl>. 

2 

In view of the semi -rigorous nature of the arguments 
developed above the ultimate test is the experimental results. 
The successful use of the ZTR method in NQR applications is 
lindted only by the ability of the spin-spin interactions to 
achieve a‘ state of quasi -equi li brl urn during the time of the 

spin-locking pulse. For a strongly inhomogeneous system, where 
neighboring spins are isolated from one another, the dipolar 
coupling cannot bring about this situation in the second 
interaction frame. The response of the spin system following the 
long pulse can then not be expected to reproduce faithfully the 
true FID. If there is inhomogeneous broadening but still 
effective spin-spin interaction either throughout the entire 



CHAPTER IV 


SPIN 3/2 NQR IN POWDERS * MULTIPLE-PULSE SYMMETRY AND 
APPLICATION OF AVERAGE HAMILTONIAN THEORY CAHT3 

In this Chapter we present a theoretical study of the 
response of a spin 3/^ powder sample to several 

multiple-pulse sequences which are frequently exploited in solid - 
state NMR. Also proposed are some suitably adapted multiple-pulse 
sequences with their possible influence on internal interactions 
of a powder sample. The guideline for designing our multiple- 
pulse sequences has been their symmetry property with respect to 
the spin 3/"S NQR powder sample. A critical anal 3 rsis of the 
cyclicity criterion of these sequences in the present context has 
been made. It is interesting to point out that the synmietry of 
some of the multiple-pulse sequences proposed by us is 

independent of the crystal orientation parameters. Also some of 
them preserve their cyclicity irrespective of the flip angle 
chosen for the pulses in the sequence. The Average Hamiltonian 
Theory CAHTD within the framework: of Magnus expansion has been 
used to study the effect of some t 3 q:^ical Ce.g. PAPS and WAHUHA3 
cyclic multiple-pulse sequences on the system's response. Our 

Ir 

t-hieor eti. ca.1 findings for PAPS have been demonstrated by a Cl 
experiment on KCIO^ powder. 
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IV. A. Some useful Properties of the Operators L, M and N 

Before we enter the main theme of this Chapter, it is 
pertinent here to point out the mathematical properties which 
characterize the operators L, M and N defined in Eqns II.3£a. 

II. 60a and III. 11, respectively. The subscripts and 3 on M 

and N distinguish the different phases C<p^ involved in k:,l,m,n»p 
and q in Eqn 11.60. The properties defined here, therefore^ are 
true for these subscripted operators also. 

It has been observed by us earlier C Chapter HID that the 
operators L, M and N in spin 3/3 NC^ are analogues of the 
corresponding angular momentum operators I^, I^ and I^, 
respectively, in spin l/S NMR. L. M and N have the additional 
properties of carrying information about the crystal orientation 
paranusters, the asymnvetry parameter and the phase of the 
externally applied rf field. We list here their most important 
properties : 

CaD The commutation relations C reproduced from Eqns 

III. 13a-dD are : 

[ N, L 3 = M , CIV. laD 

t N, M 3 = r^L . CIV. IbD 


f M. L 3 


N 


CIV. IcD 



13a 


where 


r* = fc* + 1 * + + n* + + q* 


CIV. IdU 


Also 


i 1 


CIV. aa:> 




i 24 
=-4 ' 


CIV. ab> 




1 2 . 

4 ^ 


CIV. ac> 


Eqns IV. aa-c lead to the following properties of their 
respective e.xponential operators ; 


expCiLtD = lcosCt/''a 5 +aiLstnCt/^ 2 > 


CIV. 3a5 


expCiMt:> = lcoshCtrxa:>+^i^^sinhCtr/a:) . CIV. 3 b 3 


.aiN 


expCiNtO = 1 cosC -p sinCtr/"l£^ . CIV. 3 c J 


Eqns IV. Sa-c lead to the follov.ring form of the evolution of 
the density operator : 

pCty = expC “iLt^ pCo]>expClLt;> 


f * -1 - 

pCoj>-*--|isint+ 4 sin Ctr/a:)L /[ pCo!) , L] 


CIV. 4 a:? 
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pCiJ = expC-iMt^pCo^expCiMt? 




si nhC tr/3:> 

= pCo:)+-{CiXr>sinhC tr) + MV [pCo^.M] , 


LnhCtr/3:>_' 

^M. 

CrxeD^ 


CIV. 4b> 


pet!) = expC-lNt3pCo>expCiNt:> 


sinCtr/^D 

= p<lo2> +-{Ci.'^r;)sinC tr5 + N 

Cr/^D^ 


} [pCo:>,N] . 

CIV. 4cD 


We will use ail the abovementioned properties of L,M and N in the 
Sections that follow and also in later Chapters. 


IV. B. Cyclic, Periodic and Aperiodic Perturbations and Average 
Hamiltonian Theory CAHTl 

The general trend in much of the recent magnetic resonance 
i-esearch is to prepare a system, or to let it evolve, under an 
ext er nal 1 y cont r oi 1 ed Hand. 1 ton! an or under a Hand. 1 1 onl an wi th 
certain synnnetry or transformation properties Csee Figure IV. 15. 
It is necessary on many occasions to implement a specific desired 
Hamiltonian, perhaps different from the natural unperturbed 
Hamiltonian, of the system. In other words, whereas the system 
might naturally evolve under its propagator UCt5,sve apply a 
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Figure IV. 1 


Natural and externally modified evolution of 
state of a system Cl 623 


Ir 


the 
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per t/urb<it-i on so that the system evolves under a diriferent 
propagator UC ti , due to a Hamiltonian ai'riving perhaps at a 

different state at time t. The perturbation needed may be a 
coherent one, such as a sequence of rf pulses, a mechanical 
rotation, or perhaps an incoherent one such as heating the sample 
or applying pressure to the system. In order to exploit the 
symmetry or transformation properties to its maximum, we shall be 
concerned here about the design of cyclic and periodic sequences 
of coherent pulse perturbation schemes. We, therefore, define 
below the terms cyclic, periodic and aperiodic perturbations, 
i • Cyclic Perturbations 

A cyclic perturbation has the property of returning the 
system to its original state after the cycle is complete. In 
other words, if the state of the system at the initial time t is 
specified by pCtZ>, and t is the cycle-time of the cyclic 

C 

perturbation scheme, then 


pCt + t:> = pCt> 
c 


CIV. Si 


A cyclic perturbation can be depicted as shown in Figure 
IV. a. Eqn IV. S implies that, when U is cyclic, 

pCt+t i = UCt , ti pCti UCt , ti”‘ = pCti . CIV.6i 

c c , c 


Eqn IV. 6 makes it evident that, for a cyclic sequence of rf field 


perturbations , 
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Figur e IV. £ 


Evolution of a system under 
Perturbation 


% 



cycl ic 
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t 

c 


UC , to = T exp [ -i 


•V * * > 

Sf ,ct )dt 1 

tf ■* 


t 


1 


CIV. 7^ 


£' Per i odlc Per t urbati ons 

The requirement that a sequence be periodic is : 

96 CtO = fit’ Ct+t 5 , CIV. 

rf rt p 

where t^ is the periodicity of the Hamiltonian. In other words, 

if a perturbation keeps recurring at fixed time intervals of t , 

p 

it is periodic. Pictorially it can be depicted as in Figure 
IV. 3Ca>. 

3. Aper i odi c Perturbations 

The perturbations that do not have the aboveraentioned 

property of recurring at fixed tiaie intervals Cthe characteristic 

time t is absent!) are aperiodic. This is shown pictorially in 
p 

Figure IV. 3Cb>. 

If a cyclic perturbation sequence is repeated more than once, 

each cycle constitutes a period and the perturbation scheme also 

becomes periodic, with is noted that there may 

be sub-time-i nter vals t , with nt = t = t , where n is a 

h h p c 

non-zero positive integer, when the perturbation is periodic Cbut 
need not be cyclic^. Also, it is very clear from the above 
discussions that periodicity does not imply cyclicity and vice 
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Figui'e IV. 3 ; CaZ> A periodic perturbation. CbD An 

perturbation [i0£3 


aper i odi c 
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versa. When an rf pulse perturbation scheme is both periodic and 
cyclic , then 


se ^ CtD = ^ Ct+Nt :> 

ri rf c 


CIV. 9a> 


U , CNt ,t:> = 1 

rf c 


CIV. 9b:5 


where N is a non-zero positive integer. It is to be noted that 
Eqns IV. 9a and IV. 9b do not impose any constraint on subperiods 
and Cor> subcycles. 

When the perturbation on a system is cyclic it is enough to 
know the behavior of the system over one cycle. For subsequent 
cycles the same behavior is repeated so that [16,633 


p CNt > = [ p Ct :> 3 ” 

c c 


CIV. 10> 



The theory that best accounts for the design of specific 
Hamiltonians Cand therefore propagators2> under coherent 
perturbation is Average Haiti! 1 toni an Theory CAHTD. Suppose the 
Hamiltonian is time -dependent Cowing to, for example, an applied 
per tur batiorD , as depicted in Figures IV. 3 and IV. 4, A question 
that arises is : can we find a time-dependent HaBiiltonian 
which can induce evolution of the system through UCt2> to the same 
state as XCt^ would induce at time t The answer is yes ♦and 
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that 8C is the Magnus HaJiU. 1 tonl an given by 1 16-19. 57,58, 163, 1643 : 


PC 






CIV. 115 


where 


♦ 


ae 


foi 


1 

t 


t 

-io 


civ.ia:> 


is the average Hamiltonian, The next two terms in this series 
are : 


kj‘ ft 


flIXt *5 . dt"df 


CIV. 13 > 


and 


it J J J S^fC t’ 5 . [SifCt'O.I^It* * ‘ 5 ] 3 


o o o 


+dect' “5 . [ PKfO 


,£*«:t*>3 3 J 


dt* * ‘df'dt ‘ 


CIV. 145 


If Eqn IV. 11 converges rapidly’, so that PC are small for k O, 
then the average Hamiltonian provides a good description of 

the system and other are correction terms. The problem of 
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convergence of Magnus expansion has been discussed in the 
literature ESI >64# 165»1663 . By applying a perturbation we havie» 
therefor e» taken a system which would have evolved under some 
Hamiltonian St and caused it to evolve effectively under as in 

Figure IV. 1, The trick is to realize the desired by 

appropriate perturbations that can be implemented expert mental ly> 
such as employing PAPS and WAHUHA sequences to be discussed in the 
next Section. Note that* in general » St and depend on t * but 

if Stity is periodic and the evolution is cyclic* then the same St 
and hold at all integral multiples of t^. If StCty is not 

periodic the approach is still useful* since one can effectively 
impose an St which produces a particular desired final state pCt> 
at time t CFigure IV. i>. It is interesting to point out here that 
the average Hamiltonian theory for cyclic perturbation schemes is 
related to Berry’s geometrical phase [leT-lSSJ. In the Section 
that follows* we shall adapt the ideas discussed here for the 
particular case of spin 3/^ NC^ to design cyclic rf pulse 
per t ur bati on schemes . 


IV. C« Cyclic Multiple-Pulse Sequences f or Spin NQR of Powders 

^ • The Phase A1 ter naied Pulse Sequence CPAPSD 

It is a two-pulse cycle in which the phases of the pulses 
differ by n. It may be represented as 




CIV. i55 
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li is seen from above that the cycle time t = 4 t. We shall 
examine below its cyclicity. In the 6-pulse limit we obtain, for 
the cycle pulse propagator, 

U = expC -1 expC -1 , Cl V. 1 6:) 

where N is given by Eqn 11.60. The subscripts i and z on N 
represent the phases ^+rr and 0, respectively, of the pulses in the 
C 2 A:le. Since cosC0+fr> = -cos0 and sinC0+rT}> =-sln0, Eqn IV. 16 
reduces to : 

U = expC-ii9N > expC+i©N y = 1 . CIV. 171) 

Therefore PAPS is cyclic for any value of 0 in IV. 15. From the 
above discussion it is clear that all the following cycles 
represent experimentally feasible PAPS sequences for spin 3^ NCS^: 


CID Cn^y ^,—tT-Cey -ar-C&y -rl- , civ. 18a> 

0 JT o Jn 

Cii> Cfr/a:)^.-|T-C6>^-£r-C8>^-Tj-^ . CIV. 18b> 

. CIV. 18c3 


& civy cn^ay 


fr-c&y -ar-cey , -r\ 

^ TTz-z an/z Jn 


CIV. 18d:> 
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It is also clear from the above discussion that the 
raul ti pi e~pul se sequences in Eqns IV. IS remain cyclic irrespective 
of the orientation of the individual crystallites in a powder 
sample. This conclusion for spin 3^ NC35 is similar to that of 
Osofcin [S3, £5, 31] in the context of spin 1 NC^. It, therefore, 
seen\s to be true for any spin I. Also S and , the flip angle 

of the cycle pulses and phase of the preparatory piii. se, 
respectively, can have any value. Usually one takes & = rT.x£ and 
<p' = O. 

2. Tlie Waugh -Huber -Haeber 1 en CWAHLMAD pulse Sequence 

It is a four -pulse cycle with cycle time t = 6t and is 

C 

represented as ; 








-r-cey . -ar-cey. ^ -r-ce:). -r 


n 


CIV. IQJ 


In the 6-pulse limit, the pulse propagator for this sequence is ; 

U = expC-i©N :)expC-i©N :)expC-i©N DexpC-iSN D , 

CIV. £0> 

where, as before, N is given by Eqn 11.60. 'line subscripts i.z,3 

— 4 

ai'jd 4 on N represent the phases and of the 

pulses, respectively, in the pulse sequence. Again, since cosCfj^H-jx^ 
= -cos 4> and slnC^+rt:) = -sin we have 
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1 . 

CIV. £i:> 


Ilierefore WAHUHA is cyclic for any value of 0 in IV. IQ. One of 
the examples of VAHUHA for spin 3/2 NC^ is 


Crr/2D 


4’ 


-{ 


T-C8D -t-c 8:> -2T-ca:) , -t-c8d 

TT nr/a anr/a 


o 4n • 


CIV. 225 


Aiiother example is : 


CTr./25 


0 


{ T-Ce5 , -T-C85 , -2 t-C85 , -r-COy , -t I 

SfT./* rr./* 5fT/* rr/* In 


CIV. 235 


The Car r -Fur cel 1 CCP5 Pulse Sequence 

This pulse sequence consists of two pulses in the cycle, both 
with a flip angle of rr and having the same phase. It may be 
represented as : 

The cycle time t is 4t. 


The cycle propagator is 
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U = e>cpC-i nN> expC-i7rN:> CIV. £51) 

= expC -£i ftNy 

iSi S 

= Icos Cfrr>+ ~ — sinCrrr^ . CIV. S6> 

Eqn IV. 26 is obtained using IV. 3c. Iv is clear from Eqn IV. 26 
that Carr-Purceil sequence is, in general, not cyclic for a spin 
3/^2 powder sample. However, it will be cyclic for all those 
crystallites for which r = 1. In the Carr-Purcell-Meiboom-Gill 
CCPMG> modification of the CP sequence, the phase difference 
between the preparatory pulse and the cycle pulses is rz/^. As 
Eqn IV. 26 is independent of CPMG pulse sequence is also not 
cyt-lic, in general, for spin 3/^ N<3$. 

'4^- The Ost rof f — Vaugh COWD Pul se Sequence 

This multiple pulse sequence consists of a train of 
pulses. The first pulse and the rest of the pulses in the train 
have a phase difference of This sequence may be represented 

3lS 

Cn^zy -T— rT-CTr/C2> -tT . CIV. 27‘i 

o n/z Jn 

III NMR, four of the pulses in the sequence constitute a cycle, 

i 

with cycle time t^ = 8 t. We shall examine here the case of spin 
3^ NQR. The cycle pulse propagator for IV. 27 is 
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U = exp N)exp C~-^ N)exp C ~ N>xp C N) 

CIV. aey 

where N is given by substi t-ut-lng 0 = rrxfi in Egn 11.60. Therel'cre, 

U = >xpC -£jTi ND 


= 1 cos Cnry + -p— sin Crrr:? . CIV. £9> 

where the last step has been obtained using Eqn IV. 3c. Therefore, 
from Eqn I V. S8 it is clear that OW is cyclic Cwith cycle time t 

c 

= 8t> only for those crystallites in a powder spin 3/^ sample that 
have their crystal orientation and asymmetry parameter such as to 
make r =1. However, as will be seen in Chapter VI, a different 
approach based only on the periodicity criterion for the pulse 
perturbation can still be made for handling the problem using AHT. 

In the following Section, we examine the effect of some of 
the cyclic multiple-pulse sequences on the internal Hamiltonians 
of spin 3/''S powder sample. Our approach is based on AHT for 
cyclic multiple-pulse sequences within Magnus expansion framework. 

IV. D. Response of Spin 3X2 Powders to Cyclic Multiple-Pulse 

Sequences 

1. general Remarks 


When a spin 3/''2 NOR powder sample is irradiated by a cyclic 
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multiple-pulse sequence » the behavior of the system is described 
by AHT, provided we observe the system only at specified 
observation windows, 1. e. , only after a cycle Is complete. By 
making a Judicious choice of cyclic multiple-pulse sequences we 
can introduce time dependence in the 1 ine— broaden! ng internal 
interactions in such a way that their average Hamiltonian commutes 
with the initial density operator pCcO prepared from the 
equilibrium density operator p by a preparatory pulse fl70],i.e. . 

C pCo:>. ^3=0 . CIV. 30!) 

pCoD thus becoQjes a constant of the motion during the cyclic 
interaction with the rf, and the line-broadening internal 
interactions under consideration do not contribute to the 
evolution of the system. In effect, we achieve llne-i'iarrowlng. 
In this Section we calculate the average Hamiltonian using Magnus 
expansion for the line-broadening Internal interactions , namely, 
efg 1 nhomogenel ty and heteronuclear dipole coupling. The 

commutators of these average Hamiltonians with pCo5 are then 
evaluated. Homonuclear dipole coupling is very weak and is 
discarded for this study. Also, conclusions regarding torsional 
oscillations will be the same ^.s that for efg inhomogeneity. 

2. The Effect of Cyclic Mul tipi e-Pul se Sequences on the efg 
I nhomogenel t v and Torsional Qsci 1 1 ati on Haml 1 toni an 
The major line-broadening interactions in quadrupoiar solids 
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are efg inhomogeneity and the effect of torsional oscillations. 
The Hamiltonians corresponding to both these interactions have a 
common structure and may be represented as CEqns 11.33 and II. 611): 


A, = A KL 

V Q i 


CIV. 31> 


where A^ is a parameter that takes into account the shift of 
resonanc frequency from the centre frequency due to efg 
inhomogeneity or due to the secular part of the torsional 
oscillation Hamiltonian. 

Ca5 The Phase Alternated Pulse Sequence CPAPSD 

We shall now consider the pulse given in IV. iSa and its 
effect on the Haaviltonian given in Eqn IV. 31. The pulse cycle in 
IV. 18a is : 


{ 


T -C &D -St -C -T 

TC O 




CIV. 335 


For the average Hamiltonian, therefore, we look at the 
transformations of L during PAPS. We present the result in Table 
IV. i. 

Therefore, using Eqn IV. IS, we obtain the average Hamiltonian 
for efg iahomogenei ty and torsional Hamiltonaian during the 

cycle IV. 3S as. 


A,«f 
L a. 


Co5 


■Si 


.Or 


= A^K-j^LcosC'^ 3+ 


iM 


sinC 


.Or 


3 


CIV. 335 
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Table IV, 1 : Transformations of L during PAP'S** 


Interval 

Time 

r f 

rf 

■-k*- 

C 

6 

1 


0< t<T 

0 

1 

L 

P'ulse 

1 

T 

-N 

expt i ©N> 

Undefined 

6 

z 


T<t<3T 

0 

1 

Lcost ©r 5 +--st nC ©r 5 
r 

Pulse 

s 

3t 

+-N 

expt -i©N5 

Undef i ned 

6 

3 


3t< t<4T 

0 

1 

L 


If the preparatory pulse has flip-angle from Eqn III. 10 

we can write the initial density operator as : 


pCo> = K IL cost© 


|l cost sin t©^r>^ 


till . i05 


The commutator is easily evaluated. After some algebra, 
the result is 

r pt o'l , A. = A, K^[ L . M] >-Jsi nt ©^r :> +si n(t r ) j- . 

CIV. 345 

1 

Eqn IV. 34 will go to zero only when ©r= ©^r -n. Therefore, 

those crystallites in a powder that have r = 1 will not see the 


**« is obtained from Eqn II. 60a by putting 0=0.- 
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effect of efg irshomogenet ty and torsional oscl nations* when 

For all other orientations of the crystallites efg and 

torsional oscillations will contribute to the i i ne'^broadeni ng. 

For a single cr 3 rstal with r =i * PAPS with 0 ==0-Tt will show a 

o 

1 i ne-nar rowi ng. 

1 The. lifAuah "^HUber --H Aeber 1 en C WAHUH A> Pulse Sequence 
We will consider the following 0301:1 e CEqn IV.^t:0 : 

lr-C0':) -r-ce:) ^ -ar-ce-:) ^ . civ . sb :) 

t Tt rx/'z svcf^z o J n 


We look at the transf orniations of L during WAHUHA. In Table 
I V. £ below we present the results . 

From Table IV. £, the average Hamiltonian for the efg 
inhomogeneity and torsional oscillation Hamiltonian can be 
calculated. W'e present the final results below : 


V. a 3 1 


“ £i 

{i -f cosC ^r y -^coshC ^r J cost Or 5 si nC 0r 5 


si nhC 6 fi-:>cosC 6 fr:> 


CIV. 36':) 


When Or = , Eqn IV. 36 reduces to 




CIV. 37) 
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Table IV. g : Trans for ftiat ions of L during WAHUHA** 


I nter val 

Time 

r f 


*%«■ 

C 

6 

1 

0< t<T 

0 

1 

L 

Pul se 1 

T 

-N 

expCi^N:) 

Undefined 

6 

z 

T< t<2T 

O 

1 

— i Q 

LcosC Or :) + -^i nC rr y 
r 

Pul se IS 

Sr 

+M 

ejqpC -i OhD 

Undef i ned 

c5 

3 

£t< tS4T 

0 

1 

■{LcoshC Or y “^-^i nhC Or y 





xcosC Or > +^— s i nC Or > 
r 

Pul se 3 

4t 

-M 

expC i ONZ) 

Undefined 

6 

* 

4r< t <5 t 

0 

1 

LcosC Or y +--si nC Or > 
r 

Pul se 4 

5r 

N 

expC — i OhD 

Undef i ned 

<5 

5 

5t< t <6t 

0 

1 

L 

It 

can be shown 

that Eqn 

IV. 36 does 

not conunute with the 


initial density operator CEqn III.IOI? for 6^0. Therefore, the 
interesting conclusion may be made that WAHUHA will not be 
effective in removing broadening due to efg inhomogenei ty and 
torsional oscillations from the NOR line shape of a spin 3/2 
powder sample. 


’‘n and M are obtained from Eqs II. 60a and III. 11, respectively, by 


putting <p = O. 
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3 • EUCSSii. q£. paps on Heteronuclear Di pol e Int-eract-ion i n. a 

Spin 3x3 Si ngl e Crystal Sample 

The calculation of the average Hamiltonian corresponding to 
heteronuclear tit pole coupling in a powder sample becomes very 
complicated. We> therefore, present here the result of our 
calculations for a single crystal samp’ e [1711. We consider the 
rf field as being applied along the x-axis of the Cg^AS and we 
assume the single-crystal to consist of physically equivalent 
quadrupolar nuclei. We give here only the final result : 

The average Hamiltonian is 


I>< hetero) 


^ - 1-4 - r^-l--4 -2-3^ 

= G Cl iCI -I )cos(~:0 

IX X X K ^ Sr 


-Cl +I :>sinC~';i>rcos('“:0 

y V £^J ^ 


*c^c I ‘ - Vi I ‘->-1 cos <®|> 

.►C I ‘ - Vi » -V si .i(%)}<=os(§|3 


+C Cl 
e» 


1-4 

z 


+I 




z 


CIV. 335 
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v-fhere 


S " " F/V3sinC-cos!r3; C^ = -F^; 

= F^Cy^lnf+cosf>;C^ = F^; and cosC . 

CIV. 39> 

CFor the definition of F’s, see L’qn II.63> 

A1 so t 

[pCoi). * CIV. 405 

0<h©tro> 

when - & r = n , CIV. 4il> 

o 

Therefore, we can conclude that, when condition IV. 41 is 
satisfied, heteronuclear dipole interaction does not contribute 
to the NQR line shape in a single crystal containing spin 3/^ 

nuclei- We infer that, for particular orientations of 

crystallites in a powder sainple, the heteronuclear dipole 
coupling is removed even in powder samples. In the next Section 
we present our experimental result of the response of spin 3/^ 

powder sample during PAPS, which essentially deiiionstrates the 

theory developed above. 


IV. E* Experliaenial Details 


In this Section 


we present our experimental result on 
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application of PAPS to the Cl NQR* of a powder sample of KCiO . 

3 

The experiment was performed using a Bruker CXP-90 NMR 
Spectrometer at the Central Leather Research Istltute Madras. 
India. Of course, the dc magnetic field had been switched off 
during the experiment. A two— pulse program, in vdilch the second 
pulse with alternating phases would occur between 1 and 4095 
counts during a loop, was written, and data were acquired at the 
end of each cycle. Tbe first pulse was set to be a *n' pulse. 
The peak transmitter power employed was :5i400-500 Watts 
corresponding to a ‘rr* flip angle of £6 psecs. A standard wide- 
line NMR probe CIO mm diameter? tuned to oscillate at a frequency 
of £0-3£ MHz was employed. The phase of the pulses in the cycle 
could be controlled in steps by using the pulse-list feature 

of the system's Z-i7c pulser. The spectral width chosen was 1/''4 t, 
where £t is the separation between the & pulses in the PAPS cycle. 
The ambient temperature was (i\ai ntai ned to be £1±1°C. 

KCIO is an ionic crvstai and has a monoclinic structure 

3 

with £ molecules per unit cell at room temperature t£.17^]. There 
is unit positive charge on the positive ion and unit 

negative charge on the C10~ ions, and the direct influence of 
neighboring ions on the field gradient can be appreciable [1731. 
We present the experimental parameters and the spectrum obtained 
on the next page. The prolongation of the signal for a time as 
long as 400 msecs is an indication of the removal of efg 
inhomogeneity effect, torsional oscillations and heteronuclear 
dipole coupling. 



Figure IV. 4* 


Time domain response of *=*€1 signal in KCIO^ powder 
C<o^ = 28.100 MHzD to PAPS ; experimental parameters 
p^secs. & = iTy A * 0.021. 


T *=s 
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IV. F. Sumnvary 

In this Chapter, we have set forth our results regarding 
useful mathematical properties of the operators L,M and N. It 
turns out that these three operators describe completely the 
dynamics of an assembly of non-interacting spin 3^'S nuclei. We 
have critically examined the cyclicity of the well-known multiple- 
pulse sequences aiid how they can be adapted for spin 3/''2 NC^ 
spectroscopy of powder samples. The Average Hamiltonian Theory 
has been used to calculate the effect of PAPS and WAHUHA on the 
major internal Hamiltonians of spin 3/^ NC2? system. Also, 

explicitly brought out in this Chapter is the manner in which the 
pulse sequences influencing the Internal Hamiltonians depend on 
crystal orientation and asymmetry parameter. Finally, we have 
presented our expert metal results on the use of PAPS in Cl 
resonance of KCIO^ powder sample. Both theory and experiment 
enable us to conclude that the efg inhomogeneity and torsional 
oscillation Hamiltonian can effectively be removed from the system 
for those crystallites having © r = ©r = tr, during a PAPS 

" O 

experiment. Furthermore, our PAPS calculations for a single 

crystal show that heteronuclear dipole coupling is also removed 

completely for flip angles - & = n the same should hold true 

in a powder at least for some orientations of individual 

% 

cr.Ystalli tes with respect to the rf coil axis. In case of WAHUHA, 
line-broadening features due to quadrupole coupling and torsional 
effects, are not suppressed. However, a WAHUHA experiment was not 


performed by us. 
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NQR AND SPIN THERMODYNAMICS IN SPIN 3/2 SYSEMS EXCITED 
PERIODICALLY IN TIME : USE OF THE FLOQUET ThCOREM 

In Lhe rest of the s-hesis we present our studies of the 
behavior of spin S/'S quadrupolar nuclear ensemble when subjected 
to the mul t i pi e— pul se scheme referred to earlier as spin— locking 
pulse. This combined system of spin S/'S nuclear spins and the 
spin-locking multiple-pulse sequence constitutes a quantum 
mechanical system characterized by a Hamiltonian with temporal 
periodicity. This is an example of a class of quantum mechanical 
systems the analysis of which is best done using Floquet Theorem. 

The use of Floquet theorem enables us to model the natural 
peri odic-ti aie -dependent Hamiltonian by a tinje -independent Floquet 
Hami 1 toni an . 

To understand the equilibrium properties of a many-body 
system with an energy-conserving Hami 1 toixi an , one often takes 
recourse to the powerful tools provided by statistical mechanics. 
However, a quantum mechanical system of which the Hajjii 1 toni an is 
time-periodic is energy non-conservative and is, as such^ 
Intractable by the methods of * statistical mechanics. It is here 
that the model Floquet Hamiltonian C which is time-independent and 

hence energy— conservlngZ) comes into its own in rendering the 

* 

methods of statistical mechanics more amenable. We 


shal 1 
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elaborate upon the abovementioned ideas in this Chapter, Keeping 
in 0 vlnd the spin 3/£ quadrupolar system with time-periodic 
Hamiltonian. To start with, however, we present a brief review of 
the applications ol Floquet theorem to quantum mechanical systems 
in general. The earlier applications of Floquet theorem in 
0 ^gnetic resonance are also pointed out in this review section. 

V. A. Applications of Floquet Theorem to Periodic Quantum 

Mecliuanlcal Systems 
i • General Remarks 

Mathematically, Floquet theorem helps us to write down the 
solutions of a system of linear differential equations with 
periodic coefficients £60,1743. The solutions of the differential 
equations of a quantized system interacting with a sinusoidal 
monochrojnatic field can be written as a Fourier series. As far 
back as the years 1883-1900, Floquet £593 studied the solutions of 
such linear differential equations with periodic coefficients, and 
Poincare £1753 considered the construction of the solutions. 
It was next used by Bloch to analyse the Schrbdinger equation for 
a periodic potential. The statement of the Bloch theorem £1763 
is in fact the Floquet theorem as applied to the spatial 
periodicity of a crystal lattice. The Floquet theorem was later 
used by Autler and Townes £1773 to obtain the wave function for 
the two-level system in terms of infinite continued fractions. 
However, the application of Floquet theorem to quantum systeins 
began to grow only after the mid— 1960s £61,178—1803. 
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In a well-known paper Shirley [613 reformulated the problem 
ol the interaction of a quantum system with a strong oscillating 
field by relating the solution of the Schrddinger equation to a 
time-independent infinite Floquet matrix. While this is literally 
a semi class! cal theory which nowhere involves explicit field 
quantization, Shirley has shown that his Floquet states can be 
interpreted physically as quantum field states. Indeed, his 
Floquet quasi energy diagram is identical to the dressed atom 
picture derived by Cohen-Tannoudji and Haroche C1813. Shirley’s 
work and its subsequent application to various problen\s by Young 
et al C18S3 was ail for f inite-basl s -set systems. A noteworthy 
development after this was its generalization to a system 
characterized by infinite Hilbert space C633. Barone et al [6£] 
hav'e applied the Floquet theory within the framework of Shirley's 
method to the interaction of stationary and moving, neutral 
two— leva! atomic or molecular systems with monochromatic linearly 
and circularly polarized electromagnetic radiation. Recently, 
Ho et al £183-1863 have e.xtended the single-mode Floquet formalism 
of Shirley to a generalized many-mode FI oquet theory yielding a 
practical non— per turbatl ve technique for the semici assi cal 
treatment of the interaction of a quantum system with several 
monochromatic oscillating fields. Ho et al have illustrated the 
theory by a detailed study of the population dynamics of a 
three-level S 3 ^tem driven by two monochromatic radiation fields. 
Also, Mar icq £1943 has applied it to solve the problem of the 
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iow-energy scattering by a vibrating molecule. 

In fact the list ol applications of Floquet theorem Is quite 
extensive. A comprehensive survey of different Floquet approaches 
for two-level systems published before 1976 is given by Dion and 
Hirschfelder E1873 which also includes a description of the nine 
Floquet theorems. Some of the earlier applications '"f Floquet 
theory to excited atomic and molecular states in strong fields can 
be found in the review by Bayfield [1883. Recent developments in 
semiclassical Floquet theories for intense-field multi photon 
processes have been reviewed by Chu [1863. 

Appl icati ons to Magnetic Resonance 

The earliest known application of Floquet theorem in magnetic 
resonance is by Bucci et al [189-1913. A strong, continuous rf 
field can cause a raul tipi e-quantum resonance, in continuous wave 
Ccwl> experiments, to shift in frequency [198-1933 in analogy to 
Bloch -Si egert effects [1823. Bucci et al [1893 using Floquet 
theorem confirmed these effects quantitatively and formulated a 
new theoretical approach which circumvents the rotating frameC in 
spin NMR3) by using the second quantization formalism [1903. 

This work led to mul tipi e-quantum e.xperlments with two different 
rf fields which involve the absorption of p photons from one wave 
and the simultaneous emission of p±l photons to the other wave 
[1913. Mar icq has used the Floquet solution for the evolution 
operator of a spin 1/2 NMR [64-663 and also the single resonant 
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frequency irradiation of a spin 1 NQR system in a multiple-pulse 
experiment [33], The periodicity referred to in these experiments 
is of the multiple-pulse sequences in the interaction frame. The 
effective Hamiltonian thus obtained has been used to • examine 
the equilibrium properties of the system [33,64-66,1051. Recently 
[1961 Mar icq has also used the method of second quantization [1971 
to study the relaxation and equilibrium of a spin system coupled 
to a radiation field. The idea of Shirley's Floquet states have 
been successfully used in the study of multi -photon Nl® [198,1991. 
A semi class! cal approach involving the development of the method 
of transformations to rotating coordinate frames, particularly 
suitable for magnetic resonance problems, has been recently 
reviewed by Series [2001. 

We now present in the following the Floquet Theorem and its 
application to study the dynamics and equilibrium properties of a 
spin 3x12 solid when irradiated with a periodic multiple- pulse 
sequence. We have adapted the method of Mar icq C 64-661 for 
this purpose. In the following Section we present Floquet theorem 
and a derivation of Floquet Hamiltonian to study the experiments 
of the kind mentioned above. . 


V. B. Tiiwe-Perlodic Quantum ‘Mechanical Systems and the Floquet 
Theorem 

The time-dependent Schrddinger equation is 


d UCt> 
dt 


-isecouco , 


CV. 
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assumption, = SKO, it follows that VCO must be a 

solution to azhrbdinger equation with the initial condition 

VCt=0>= UCt^3. Next we set tfCt!) = UCOUCtO. By operating with 

the left on both sides of Schddinger equation for 

UCt!) CEqn. '/.ID it follo'ws that WCtD is also a solution with the 

initial condl ^ ion WCt=OD= UCtD. From the uniqueness of the 

solution we must then have WCtD = VCtD or UCt+t D = UCtDUCt D. 

P p 

From the definition of the operator PCtD we have 

PCt-*-t D -- UCt+t Dexp [iCt+t D^l . CV. 4D 

P P ^ P 

Using the previous result we obtain 

pct^-t > = uct:>uct :)exp [ict+t . cv-s:^ 

p p p 

Because UCt > = expC ~i » the above two time-inc;Jependent 

factors cancel , leaving UC tI>expCii*ft!> which is Just PC t> . Thus 

PCt-^t > = PCi> completing the proof. 

P 

It is to be noted that the order of the factors on the right 
hand side of Eqn V. 3 is important. If these factors are 
reversed^ the resulting PCtD will not be periodic unless UCt D 

p 

commutes with UCtD for ail t» which need not be true. 

It may be pointed out here that in a sense 9€ is an 
effective ti me-lndependent Hamiltonian for the system. More 
precisely, if the time-scale for secular changes is large 
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compared to the period t^, the time-averaged state vector^ averaged 
over a period t^ evolves as if ^ were the Hamiltonian. 
Essentially we have modelled the natural Hamiltonian of the 
system which is ti me— dep>endent by an effective time— independent 
Hamiltonian, namely, the Floquet Hamiltonian. 

We now substitute Eqn V. 3 into the Schrodinger equation 
CEqn V. 1>. This gives 

= -iS«^t>F<t5 +iPCt>^ , CV. 6> 


with PCO> = 1. Eqn V. 6 is a system of linear differential 

equations. A perturbation scheme is obtained by invoking the 
following power series expansions. 

00 

PCt> = V PCO, CV.7> 

n=i 

00 

se = ^ , cv. 65 

n:=l 

and v/iih the substitution >• It is to be noted that Eqn 

V. 7 is a time-dependent power series in X» whereas Eqn V. 8 is a 
time -independent power series in X. The factor X has been 

introduced to keep track of the various orders and is set to 1 in 
the end. Next we substitute Eqns V. 7 and V, 8 into Eqn V, 6. After 
some algebra we obtain the main result : 
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P CtO 

n 


t 


-i 


f cfD 

Jo ^ 


n-l 

k=l 


f 


CV. 9> 


and m 
n 


t 

£ I l^t •.)P^_^Ct*> 

P Jo ^ 


n- 1 

k=i 


CV. 10> 

The are obtained with the additional requirement that 

each P Ct5 be periodic with period t . It is to be noted that 
n p 

the two zeroth^-order terms are P = 1 and X O. 

o o 

In magnetic resonance applications, all the above 

calculations are assumed to have been done in the interaction 
frame of the main term in the total Hamiltonian of the system. 

Mar icq C643 has shown the equivalence of Magnus expansion 
terms CEqn IV. 11> and Floquet Hamiltonian terms given by Eqn 
V. 10, provided the latter is evaluated at times t = t^. The 
correspondence between the two series is 

^ ^ , cv.ii:) 

n<F) M 

where ^ is nth order Floquet Hamiltonian and the 

nCFJ " 

Cn-l!)th order Magnus Hamiltonian. 

The nth-order term in Eqn V. 10 is easily obtained and the 

expansion of PCt) is evaluated concurrently with This can 



171 


potentially be used to follow the evolution of the system within 
each cycle. This^ in essence > removes the restriction of 
observation of the system only at stroboscopic points in the 
observation windows. Mar icq [643 has given a proof of the 
fact that the series for expC'-i^Sft D converges in mathematical 


sense. Thus> the power series ^ ^ must al 


so converge. 


Often in many experimental situations the energy Level is not 
sharp but spread around a central value because of the various 
broadening mechanisms which we have earlier discussed. Under such 
conditions » Mar icq has concluded that in order that the Floquet 
Hand. Itonian series converges rapidly one must irradiate the 
system always above resonance. In fact» it has been shown 
by Mar icq that the decay of the quasi -stationary magnetization at 
long times is the consequence of the failure of series Floquet 
Hamiltonians Cor Magnus Hamil tonians> to converge. We shall 
adapt Maricq*s arguments for our particular S 3 rstem under study. 

In the next Section we will show how the model Floquet 
Hamiltonian can be used to study the equilibrium propeties of a 
periodically perturbed quadrupolar spin 3/^ system. 


Equilibrium Properties from Time Averages in Spin 3X2 Solids 
The purpose of this Section is to derive the expression for 
the equilibrium magnetization of a system of spin 3X2 nuclear 
ensemble which is subjected to a time-periodic external 
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perturbation. The corresponding Hamiltonian will be 

explicitly ti me— dependent. This precludes the direct application 

of ensemble averaging procedure to arrive at the eguilibrium 

properties of the system. One approach to surmount this 

problem is to approximate the time— periodic systeni by another 

that evolves under a time— independent effective Hamil tonian» 

namely^ the Fioquet Hamiltonian described in the last Stection. 

The disadvantage of this procedure is that one is then 

confronted by the question of the uniqueness of 96 C Fioquet 

Hamiltonian is not unique; in fact each of its eigenvalues is 

only defined within an additive constant of nn/t » for a non-zero 

P 

positive integer , In order to arrive at an equilibrium value 
that is independent of the choice of 96 ^ an alternative route 
is taken to calculate the equilibrium magnetization from the 
time average of MCtD. This procedure is based on the following 
ar gument . 

A macroscopic collection of quadrupoiar spin 3/E nuclei in 
a solid sample will obey the mixing property C20S-S04I. The 
mixing property guarantees that if a nonequii i hr i um population of 
spin states is introduced locally* so that only a small portion 
of the sample is affected* then after a time the nonequilibrium 
state will spread uniformly throughout the sample via spin 
diffusion . A dynamical system that obeys this mixing property is 
ergodic [803,3043. In the case of a conservative Hamiltonian 
governing the spin system, this guarantees the equality between 
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the time average of the magnetization C expectation value taken 
along the rf detector coil a.xis> 


M 


lim 

T-*oo 



T 

TrBpCOdt, 


CV. ia:> 


and the ensemble average : 


M 


•q 

M. 


L 


TrCQ SIOTr<Sf(B> 
Tr 


CV. 135 


where Q is a traceless Her mi ti an operator and is defined to be the 
initial reduced density operator C immediately after the 
preparatory pulse> and is the same as pCo]) defined in Eqn.III»10 
Ct^i^ OD , Also» B is the magnetization operator already defined in 
Eqn 11.69. We explore in the following the extension of the 
time -aver age approach for calculating equilibrium properties in 
spin 3X2 NQR experiment » to a periodic time-dependent Hamiltonian. 

Calculations of the magnetization using the Floquet form 
for the evolution operator leads to an equilibrium value of : 


M 


lim 

T-#oo 


T 

BJ Tr |bPC t5 expC -i Set'j pC o5 «xpC i P 1 1 > Jdt 


CV. 14> 


Here 


UC t5=PCt5e>cpC-iS^5 is the propagator for the time-periodic 
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Schrttdinger equation CEqnV.l>.By Interchanging the order of p'^'ctD 
and the remainder of the expression in the trace, vw? can 

consider the magnetization to arise from the product of ts^ 
terms the quantity ltCt5 = expC-i^>pCO>expCi^D . in which the 

spin system appears to evolve under an energy— conservative 
i t. onx ©.n X A.n.ci t.i m<0***“pfS'riociic obs^&r VclIdI 

p'^ct:>BPct5 . cv. IS) 


The limit, as T — ► oo, of the integral in Eqn V. 14 is non-zero 
only if the integrand has a component that is constant. This is 
because, as the time variation is periodic, the integral will be 
sweeping two equal and opposite areas such as those of a 
sinusoidal function. Hiis results in the value of the integral 
becomi ng zer o . 

Whether the integral in Eqn V. 14 will have a constant con^o- 

nent^ depends on the Fourier spectra of the operators P^Ct)BPCt) 

and .Sttt). Because PCt) has a periodicity of t , we can expand 

P 

P*^Ct>BPCt) in a Fourier series, obtaining : 


p”^tt)BPCt> = y expCikanrt/t^) , 


tv. 16) 


whereby the spectrum is discrete. In order that a particular 
frequency component ® jc * contributes to the time average, the 
operator i<itt) must include a component that oscillates with 
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frequency -k:8n^tp. Application of the spectral theorem to 
provides the following expansion for ^t> : 


^t:> 


exp^-i t J\d\ j \>< X 


p pC 05 expC i t JXdX j XX X 



CV. 175 


where S€ is assumed to have a continuous spectrum of eigenvalues 
X. Now» since 


{|XXX|}'^ “ 


jxxxp 


CV. 185 


therefore. 


exp l^-i t JXdX |X><X j J = ejqsC— itX5-^J^dX jXXX 


CV. 195 


Antd Eqn 


V. 17 


JtCt5 


reduces to 

= expC - i tX5 j XX X I ^pC 05 -jj^X } XX X j ^xpC 1 1 X5 

= JJdXdX exp^-iCX-X 5 t^|XX X jpC05 jX XX | . 


CV. 805 
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From Eqn V. 80 it is apparent that the range of frequencies 
is dictated by the differences between the eigenvalues of and 
the line shape is determined by the matrix elements of pCOD or 
equivalently of Q in the basis of the eigenstates of X Cthe 
line shape is the Fourier transform of the autocorrelation 
function of MCt!) [SOS]. Now the counter-rotating component from 
Eqn V. SO is given by 

etxp l^-iCX-X Dtj- = exp |^-ikSnt/^t^| » CV. Sl> 

i . e. K = kSrr/t 

P 

or X - = X.* . CV. SaD 

p 

Combining the term in Eqn V. 16 that is proportional to 

expCikSfrt^t D with the counter -rotating component given by Eqn 
P 

V. SS, we get from Eqn V. 20 the equilibrium magnetization as 
= e’ rdX<X-kSTT/'t^}Bj^|X><X|pC05 |X-kSTr/t^> . 


CV. S33 

1 

The prime over summation indicates that the sum over k includes 
only those terms for which X-kStrXt is an eigenvalue of This 

is an alternative way of saying that we are taking only the 
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diagonal elements to find out the trace which gives us the 

magnetization vector. The terms with k O represent the 

contribution from transitions between states with different 
energies. The trace in Eqn V. 23 will be non— zero only when 

P 

The* uso ot Ecjn s/. 20 requires knowledge of Lhe eigenst-a.'tes of 
This is readily obtained for systems in which a small groups 
of spins interact strongly among themselves but weakly with the 
rest of the sample. In that case the matrix operators that 
describe the interactions within the group are of a manageable 
size and can be diagonalized. However > in more general 

cir cumstances » where strong interaction^ between many spins are 
involved the above expression is not very convenient. 

Eqn V. 20 for the equilibrium magnetization can be transformed 
into a more appealing form under the assumption that only the k =0 
term in the Fourier series contributes signi f icaintly to the sum» 
i.e. » only the constant components of P CtJBPCtD and ^tD are 
significant. This assumption is valid under the following 

circumstance: IX --X I < &nyt for all eignevaiue pairs of If 

* i j « p 

this is the case» then the spectrum of J«Kt> lies entirely within 

the frequencies -fT/"t and , whereby X~k2n/^t is an eigenvalue 

p p P 

of ^ only for k = O. Eqn V. 20 reduces, under these circumstances, 
to * 


= J dX<\l|ptt>BPCt>J^^jX><\jpCO-l jX> 


CV. 245 
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v^here is the zero Fourier coefficient of the Fourier expansion 
of p‘**Ct:)BPCt> , i.e. 

T 

" T I P'^ct^BPct:»dt . cv.a5> 

O 

Each of the operators ■|p”*’c t>BPCt> and pCO> can be written 
as the sum of two terms : one operator that commutes with Se and is 
di agona.! in t^hc jX)> basis i and a second bKab is contpl > 01.01 y ofiT 
diagonal. Clearly, the diagonal part of pCO> is the equilibrium 
density operator. We consider for a moment a conservative system, 
i.e. , one with PCt> =1. If the energy is the only constant of 
motion, two important properties hold true : First, any operator 
that commutes with the Hamiltonian must be expressible as a 
function of it, and second, the system is ergodic £SO3>2043. In 
the case of a time-periodic Hamiltonian, plays the role of the 
constant of motion, although it does so only at intervals of t . 

p 

Nevertheless, we use the first property to write the diagonal 
portion of the initial density operator as 

diag-JpCO^J- = = fC^J = E fCX:>f^><Xj . CV.S65 

If there were a second constant of motion, p should be expanded 
in a basis of simultaneous eigenfunctions of the two constants of 
motion, and therefore would be a function of both the 
operators. The second property enables us to use the laws of 
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statistical mechanics to determine the function fCX>. This is 
achieved by maximizing the entropy S, of the system 

S = - , CV.STD 

gubject^ to the constraints that the density operator is 
nor!Balized» i«e. § and subject to the energy constraint 

<E^> = TrCS^p^^3. Using the methods of undetermined multipliers 
CS063 we obtain the result as 

fCXD - t expC-(? X5S 1-R X , CV. 28) 

■=- ” ' Trt? 

where Z is the partition function. The high- temperature 

approximation is invoked to expand the exponential » and the 

equilibrium inverse temperature ft ^ is. written in terms of the 

eq 

initial inverse temperature ^ft^ from energy conservation £1331 
Eqn V- £8 is substituted into Eqn V. 26 to find p ♦ which in turn* 

eq 

is used to replace the diagonal part of pCO> in Eqn V. 25 to give 

cv. a9> 

This is easily rearranged into the desired form 



leo 


M « -«<^Tr^CTr<p'^COBPCt:» « 

•q ' -av CV. 30> 

Trflf * 

and in the application of this equation in the next Chapter to 
spin-locking experiment in N<lR of spin 3^ it is normalized by 
dividing by the initial magnetization M. given by 


M, = - /?.Tr{QM} 

X- V 


CV. 315 


Theref or e , 


t 

CTrfleCO (Tr<P Ct5BPCt5> Sf) 
CTr^5CTr<2B5 


CV. 325 


It is useful at this pointy to indicate some relevant 
features concerning the above derivation. Eqn V. 32 is quite 
similar to the ensemble average given in Eqn. V. 13; simply 

the magnetization operator B is replaced by i3> , the 

av 

average taken over one period in the numerator. This similarity 
underlies the relevance of the effective Hamiltonian, which is 
responsible for the evolution of the system over one period, in 
determining the equilibrium properties of the spin system. On 
the other hand, the appearance of the periodic operator PCt^^ 
signifies that the behavior qf the spins during the entire 
period must be accounted for^ in order to predict the true 

•q V 


beha vi or of M 
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The solution for given by Eqn V. 3a is valid only under 

the condition that the specturm of ^ lies completely within the 
frequency interval C -n^t^ , n/t^5 . This is the same condition that 
applies to the truncation of the series expansion for Sf in Eqn 
V. 10. Therefore, in the case in which we wish to use the 
pier t. urbaLi on solutions for the Fioc^uet form of the propagators the 
condition tlct allows the series to be truncated also allows us 
to use the result of Eqn V, 32 to calculate the equilibrium 
magnetization. Of course, if an internal interaction of the S 3 rstem 
causes the line to spread beyond ± r/t , then Eqn V. 32 applies 

p 

only approximately. 


V.D* The Ki-ylov-Bogollubov-Mitropolsky CKBIO Method for CMbtalning 
Average Haiail ionlan: Cojnparison with the Magnus Expansion 
and Floquet Theorem Methods 

The substitution of a time-dependent Hamiltonian by a 
corresponding time-independent effective or average Hamiltonian in 
solving different physical problems is of much current interest 
not only in magnetic resonance studies £29533,63,66,208-2103, 
but also in coherent optics £1863 and solid state physics £2113. 
A careful evaluation of the various techniques employed for 
obtai ni ng aver age Hami 1 ton! an , there! ore becouies per tx nent . 

Choice of a particular technique closely depends on many factors. 
The average Hamiltonians obtained are in general approximations 
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to the corresponding natural time-dependent Hamiltonian and are 
usually represented as series expansion. The convergence of 
this series solution , both in the mathematical and practical 
sense Cfrom the viewpoint of experiments^ , needs to be studied 
thoroughly. Also the mathematical conditions imposed on the 

existence of average Hamiltonians in various techniques need to 
be kept in mind so that» in a given exper t mental situations the 
technique of choice describes the true experimental 

behavi or . 

We confine ourselves here to three averaging procedures 
employed so far in explaining and designing magnetic resonance 
experiments. Two of these techniques > namely* Magnus expansion 
Csee Chapter I\0 and Fioquet theorem Csee Section B of this 
Chapter > based average Hami 1 toni ans have already been described. 
We now give a brief account of the KBM method and then examine 
the relative merits of the three techniques in the context of 
pulsed magnetic resonance experiments. Our description of this 
technique is based on the account given by Mehring CiBI which 
is a simplar version of the original derivation given hyBuishvili 
and Manabde [693. 

The KBM averaging method is connected with the existence of a 
certain change of variables^ whereby it becomes possible to 
eliminate the explicit time-depfendence with any desired degree of 
accuracy relative to a small paratmeter » i.e. » it permits the 
exact solution of the problem to be approximated by the solution 
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of a corresponding conservative problem much simpler than the 
original one. It is to be noted here that the KBM method is 
the most general way of setting up an average Hamiltonian 
because it does not depend upon the nature of the 
time-dependence of the system's true Hamiltonian. It is 

therefore useful for problems v^ere stochastic averaging or 
aperiodic perturbations are involved. Its use for the particular 
case of time-periodic systems may now be illustrated. 

Since the excitation of the system is periodic both the 
Hamiltonian and density operator can be conveniently e^q^ressed 
by the Fourier series 


9Kt5 


n 5=00 



H expCico tD » 

n n 


CV, 33> 


and 


nm=:00 


p(It> s: \ P Ct]>expCia> t> ; 

lu m 


nrisr-O 


v/her e 


CV. 34> 


r ^ 

H I ^ dtSKt])expC— t5 ; o> ■ n. CV. 353 

n t I ^ ^ 

p J p 

^ o 

Also p Ct> corresponds to thd slow motion evolution. Inserting 

m 

Eqns.V. 33 and V. 34 into Li ouvi 1 1 e-von Neumann equation CEqn II.73:> 


1 eads to 
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p Ct>+ico p Ct5+i 
w mm 






= O . 


CV. 36> 


The zeroth component p^Ct>^ known as the slowly moving part of 
pCt>, can be easily derived from CV. 365 and is given by 

~ p/t> i H^p^Ct:>+iJ H_j^p/t5 = O . 

CV. 37:> 

In a zeroth order approximation the sum over ki^O can be neglected 
and we obtain 

p Ct5 = expC-itH Dp Ct=OD , CV. 38D 

o o o 


where 


H = i r P dtSKtD , CV.39D 

o t J 
p o 

according to Eqn V. 35. Eqn V. 39 is the same result as was 
obtained earlier for the zeroth order term using Magnus expansion 
CEqn IV. 12D and also using Floquet theorem CEqns V. 9 and V. lOD. 
The deviation, however, occur’s in the higher order corrections. 
To obtain the higher -order terms it is convenient to Laplace 
transform Eqns V. 36 and V. 37, whereby we obtain 
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kVo 


C V. 401) 


for m ^ O, and 


Sp CSD+IH p CS>+iy H p CSD 

o o -k^lc 


pCci'J 


CV. 41 > 


where 


p CS> 

W 


r 


dip CtZ) expC-si> 


CV. 42> 


and where the initial condition p Ct=0>=0: p Ct=0>=pC05. We are 

Tn?*o o ~ 

essentially looking for an expression of the form 


cs+ii?ip cs> =pCo> 

o 

where 


CV. 43> 


H 




which has the solution 


C V. 44:t 


p Ct> = expC-itH>pC0:> , CV. 45:> 

o 


for the time evolution of the ‘secular’ part of the density 

i 

operator. The higher order terms are obtained by solving Eqn V. 40 
for p CS) and inserting it, after redefinition as pCS>, into 
Eqn V. 41. For the secular motion part C jSj<<co^> ^ 
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replaced by and we gel the first two terms of the series 

in Eqn, V, 44 as 


H » H 


CV. 46> 


-J 

1C To 




-k k 


C V. 47> 


Expressing H in the notation of Eqn IV. 13, we get H 


r;«J 




H 


<t> 


o ^ L 


sect :>-se' 


x:s<o> 


CV. 47a> 

Comparing Eqns IV. 13 and V. 47a, it is evident that both 
the average Hamiltonians are equal if = H =0. If, however, 

O 

H p* O there Is a ‘‘non-secular * contribution to which 

o 

leads to improper results in the Magnus treatment. This will 
be illustrated below with one example : 

The Bloch-Siegert Shift has been calculated using all the 
above techniques by different workers. The results of these 
calculations give the shift as follows 
Ca> the Magnus e>q>ansion based AHT C893: 
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* 


Cb:> the KBM method based AHT C80] : 


z 

CO 

<*> «= <0 — i + 

o 4w 

o 


CV. 49> 


and Cc> the Floquet theorem based AHT [61} : 


CO 



35b** 

38co® 


CV. 505 


where b in Eqn V. 50 is a parameter proportional to co^, the rf 

amplitude. 

It is clear from the above, that Magnus ejqaansion based AHT 
predicts erroneous value for B1 och-Si eger t Shift. The KBM and 
Floquet theorem results match upto second order at least in 
the form of the expression. 

Now we are in a position to compare qualitatively the three 
techniques for obtaining average Hamiltonian. Of all the three 
the Magnus expansion method seems to be most-restrictive. It 
imposes a special structure Ci.e. of cyclicity!) on the scheme of 
perturbation. Floquet theorem method depends only on the 
periodicity criterion to arrive at the average Hamiltonian. KBM 
is more general and, using it, we can set up an average 
Hamiltonian for any arbitrary time dependence. Magnus expansion 
method allows the observation of the system only at stroboscopic 
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points. On the other hand, both KBM and Floquet allow the 
observation o.f the system at all times. The convergence of Magnus 
expansion and Floquet theorem results has been discussed by 
Mar icq and that of KBM can be found in the book by Krylov and 

Bogoliubov [883, 


E* Si.iiiwiary 

In bhis Chapter we have reviewed the application or Floquet 
theorem to systems characterized by time-periodic Hamiltonians. 
Floquet theorem based AHT has been discussed and has been used 
to derive equilibrium properties in a spin 3/2 Finally we 

have made a comparative study or the various techniques ror 
setting up the average Hamiltonian. KBM appears to be the most 
general method to set up the average Hamiltonian. However, we 
have llirdted ourselves in this Chapter and the rollowing one to 
Floquet theorem-based methods only; the abovementioned generality 
of' the KBM approach, although realized towards the end or our 
work, could not be gone into rurther. 

In the next Chapter we shall conrine ourselves to Floquet 
theorem based AHT to explain the results or pulsed spin-locking 
experiments in NQR of spin 3^2 solids. 



CHAPTER VI 


PULSED SPIN-LOCKING IN THE NOR OF SPIN 3/2 SOLIDS : 
THEORY A^O EXPERIMENTS 


In ihis Chapt-er we present our theoreiica . calculations of* 
the response or spin 3/^ solids to spin—lockingr mul tiple-'pulse 
sequences, A computer program has been written to evaluate and 
plot graphically the final expression for equilibrium 
magnetization during the abovementioned multi pie —pulse sequence . 
In this program^ a Monte Carlo based method called ''systematic 
sampling* has been used to take into account the polycrystalline 
C powder > nature of the samples used in our experiments. A 
complete listing of the Fortran program appears in Appendix A. 

The experiments performed by us on powder samples of NaCiO^ 
and HgCl^ C Ci resonance^) at room temperature are then present-ed 
and interpreted in the light of our theory. 


VI. A. The Equilibrium Magnetization in the NQR of S^in 3X2 Solids 
in the pulsed Spin-Locking Experiment 
In pulsed spin-locking experiments the time domain signal is 
sampled at the echo-tops after every pulse in the sequence. 
Therefore, the effect of static inhomogeneous interactions. 
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namely. the time-independent terms in the spin Hamiltonian 
involving single resonant spins, such as efg inhomogeneities and 
heteronuclear dipolar coupling is removed completely. 

^^tions that determine the behavior of' the 
time-domain signal are therefore the ones that contribute to the 
width of the frequency domain signal in a homogeneous fashion, 
namely, the torsional oscillations of units containing the 
qua.dlr upol a,r nuclei in t^he solid samplo and aJ.so t^o a iessier 
extent the homonuclear dipolar couplings. The contribution of 
the homonuclear dipolar couplings is in general, negligibly small 
in spin solids. Nevertheless, we will include it in our 

calculations. The total Hamiltonian of the combined system of 
spin 3x2 solid and spin-locking multiple-pulse sequence in QIF is, 
therefore written as : 

ct:)N~AKL+^' , , cvi , i > 

1 DChomo), 

where N, L and Pt are defined in Eqns Il.SOa, II. 32a and 

0<homo>, 

11.64 respectively. For this particular experiment, the rf 
Hamiltonian in the <5— pulse limit Cthe first term in VI. 1.* will 
take the following form for the first pulse ; 

Se' = -B «5Ct^N ; > CVI. 2:? 

rtu) o i 

and for all other pulses Ci.e.the pulses in the sequence^ , it is : 
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se 


rf{z> 


00 - 

<5jt~TCak+13lN , 
k » o ^ J ^ 


&hl 00 

~ J2 "^expC i nnih'J , 

-00 


CVI.Si 


CVI.4D 


where ar is the inter pulse duration or period of the pulses in 
the sequence. The subscripts 1 and a in Eqns VI . 8 and VI. 3 
distinguish the phases of the first pulse with respect to the rest 
of the pulses. A is a parameter that takes into account both 
resonance offset and the average effect of the torsional 
oscillations over the timescale of the experiment. It is to be 
noted that Eqn VI . 4 is Fourier series representation of the 
spin-locking pulse sequence. 

For the choice of the flip angle © < 1 radian the QJF is 

appropriate for applying the perturbation solution for Sf and PCt>, 
The recursion relations given in Eqns V. 9 and V. 10 can be used to 
calculate the Floquet Hamiltonians upto third order. It can also 
be calculated using the Baker-Campbell -Hausdorf f theorem £141—1513 
to the evolution operator over one pcrioct : 


UC2T0=exFJ r-iC-AKL+Se’ ^ ^lexpr +1-SN lexp [-iC -AL+.^^ . ^ . 

1 D < homo > J i- 4-1 *- D<homorJ 


cvi. 5:> 

A straightforward but tedious algebra leads us to the 


**For example, the evaluation of the integral over the 
discontinueties of the 6-function has to be done painstakingly 
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Also, the periodic operators up to second order are given as : 


P Ct5 - 1 

o 


CVI. 105 


P Cf5 » - ^ N r — — expCin)Tt*/^T5 , 
i 2n * n 

-00 
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ril^o 
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cvi. i2:> 


In our calculations of equilibrium magnetiration we also 
need to have the average values of periodic operators. 

Integrating the Eqns VI. 10. VI . i . and VI. 12 o’mr the period from 
O to 2r and dividing by 2 t, we obtain the following expressions 
for the average values of P Ct>,PCt> and P Ct5 ; 


P Ct> = t 

o 


cvi. i3:> 


PCL> = O 

i 


CVI. t4> 


P CtD 

Z 


p 

la 





-akl+h' ^ 

D C howo > 



CV. i5> 


UaWt the quasi-*- equllibri um magnetization is given by Eqn V. 32» 
which is reproduced here: 



C Tr »KP C Tr t3 BPC t’J 

CTr Si^CtrCfS'J 


CV. 32> 


Taking into consideration the perturbative expansions for 36 and 
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PCl>ln Eqn V.32 and arranging lernia according to their order 
keeping the overall order only opto three, we have : 

<P^Ct>+pJ'ct>-t-p\-t>+. . BCP Ct>+P CO+P Ct>+. » 


C VI. 165 


Now 



therefore Eqn VI. 16 reduces to 


<P^Ct>BPCt>> =B+-fBP+P B-KP’^’ft^BP Ct55 V 

a.v * 2 1 i a.v J - 


CVI. 17> 


Cupto third order5. 

Substituting ^^^’=0 CEqn VI. 75 and ■CP'^C t5BP<:i5> from Eqn VI. 16 

into Eqn VI. 16, we obtain, upto third order : 


M 


eq 


M. 


Tr 


i 


QCl?f ^l-Tr 

* ^ J 


BCSC +X >+/bP +P B+CP’Ct5BP Ct> 

i 9 Z Z 1 1 < 




TrCSf + Sf > 
i a 


TrCQB5 


CVI. 195 


substituting the values of P^Ct5 and P^ from Eqns VI. 11 

shown straightaway that 


ind VI . 15 . it is 
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Ct>BP col = 

* * L i i Jav 

It [ N..B3.SJSJ 

/* — ^ f 

*ia ■^''tE ■ cvi.ao> 


Substituting Eqn VI . 20 into Eqn VI. 19, we obtain the final result 


as 


H 


«q 


M, 


Tr {«* >1. |i Tr{[ [ S.,B].SJ*} 

>3-q^.} 1 . 




0<homo>' 


TrC^ + X TrCC^y 

19 


CVI.21> 


The first term in the square brackets is the result that 
corresponds to the one obtained for the NMR of spin l/<2 nuclei by 
Pedfield [56], who applied spin thernK>dynami cs to the truncated 
Hamiltonian in the Zeeman interaction frame tthe rotating frajne in 
NMRD . The second term is a correction that arises from 

considering the entire dynamical history of the spins and not 
merely the stroboscopic observation pjoints. In other words, the 
oscillatory behavior of the magnetization vector about its average 
trajectory contributes to the equilibrium value. It is also to be 
noted here that the analogy of spin 3^2 N(^R with spin 1<'^ NMR 
[64,65] and spin 1 NQR [335 under multiple-pulse spin-locking 
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experiment is quite striking once the Hamiltonians are expressed 
in diagonalized representation of Se^ and STO formalism. 

In the next Section of this Chapter shall describe our 
procedure for computing the response of a powder sample of spin 
3X2 solid when the spin-locking multiple-pulse sequence is 
appl i ed , 


VI eB* Colwputation of Powder Average Response 

The expression in Eqn VI . 21 is orientation dependent, 

due to the fact that the quadrupoiar tensor C interaction matrix> 

is i ntri nsi cal 1 y anisotropic. In a finely crushed powder specimen 

quadrupolar spin sites occur with all possible orientations of 

their QPAS with respect to the axis of the rf detector coil. Ihe 

''powder pattern^ computer simulation of the line shape is 

well-known in other branches of magnetic resonance such as 

EPR C 212, 2133, and these techniques may be adapted for our purpose. 

Now the orientation dependence of M is in terms of & 

and S , which are Euler angles between <FAS major axes and rf 

coil axis CFigure II. I>. The overall powder response of aspin 3x2 

quadrupolar spin ensemble <M^^xM^>, is obtained by integrating 

* 

Csumndng and averaging3> the expression C which is the single 

crystal response for any arbitrary and in the whole space 

C& , 4 , >.The ranges of & and are given as 0< n and O < e^< 

lU JL# ^ 
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an Cln Figure II. I ^nd in the following derivations, it may seem 
at first sight that we are rotating the detector coil axis at 
various 0^, for a fixed crystal orientation. Ihis is however, 
equivalent to the physical reality where the -powder’ 

corresponds to a collection of single crystals that are at all 
possible orientations, which fully validates the procedure 
d«?scribed b©low>. 

In one finds the usage of various powder averaging 

techniques by different authors [1,3,161,814,0153. For example 
Petersen C1613 prescribes a procedure where averaging is done 
only in 9 space Ca^cially symnvetric cases onlyD. In the present 

lw> 

work, we take the more general case of non~ax±ally symmetric NQR 
systems. Here the averaging is done in complete C9 ^4^ ^ space. 

JLt Zw 

Considering the fact that the number density of orientations 
with respect to each 9^ is proportional to sin 9^ 12173 and in the 
(p direction all orientations are efquailly probable, one can write 

ILi 

the pov^der response as : 



L 



2TT TT 

I f (M \M, ) sin9 d9 d<p 
J J eq t. L 


o o 


tVI. 22> 


The problem at this stage amounts to estimating this 
integral. Various numerical techniques including Monte Carlo ISIQI 
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procedures are in vogue for this purpose. Recently in our 
laborator/ Cdl9,ag0] a much faster Monte Carlo procedure based on 
systema.ttc sampling Caat.aaS] has been developed, which may be 
briefly described as follows: 

In essence, our procedure considers the 'weight factor* for 
and as 'probability density functions’ Caail, Accordingly, 
Lhe normalized density function of Q can be written as 

L» 

sia & 

L 

PCd :i = ; 0 < (9 < rr , CIV.23> 

lu* cii L* 

and t/h<e normal ization facLor can be correot.f?cL labier by 
multiplying the above response by S. Similarly the density 
function of may be understood to be of a ^uniform^ type 

Generally in a Monte Carlo procedure one generates pseudo 
uniform random numbers corresponding to the integration range* 
and estimates the average of the integrand for the values thus 
generated CS163. The only difference in systematic sampling^ is 
that these uniform random numbers are generated systeraaticail y 
such that 


p =- , where i - 1,2* ...N , CVI,24!> 

i N 

It is easy to see that these deterministic also follow 

a uniform distribution and the estimates of averages will not 
change very much from those of earlier procedures. 



109 


Finally, Tor the range of Q thi» p »€? 

y t, ® thus generated are 

converted to follow the density function of © . PC^ > CEqn VI.a3>, 

L L> 

in such a way that *s are obtained by solving. 



PC© >d© 

L L 


o 



sin© 

2 


d© 

L 


CVI. 25> 


Similarly, the ’ s corresponding to the powder are obtained by 
simply scaling the uniform R^’s for the range 0 ^ i n. Thus Eqn 

VI . 2S is approximated as : 


& 


LK 




LN 


<^i> - Z 


L Is L, L 
i i 


where fCQ,<p2> = M /M, 

L eq-' i 

The above equations have been programmed for our 
calculations. The complete program listing CFortran I\0 appears 
in Appendix A. 

In the following few pages we give the theoretical plots for 

M yH versus vaj'ious experimental parameters. The specific 

•q i. 

parameters for vthich the computations were performed are specified 
in the Figure captions. The Important point to be noted is that 
in the Figures VI . 1 Ca>, CbD. Cc>, Ce> and Cgi) we see that a 
finite value for the equilibrium magnetization persists even when 
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the contribution of the dipolar Hamiltonian to X is zero. The 

irit 

interesting point is that spin-locking can be attained even in the 
presence of the torsional oscillation Hamiltonian alone. This has 
relevance to the experiments that we have performed in this thesis 
Csee Section VI . C. > The internal Hamiltonian of spin 3x3 solid 
has a very negligible contribution from the homonuclear dipolar 
interaction, as pointed out in Chapter II. In the case of offset 
irradiation, the behavior of the equilibrium magnetization can be 
explained by similar curves. Another point to be noted in these 
the difference in the & dependence for the case where 
dipolar Hamiltonian is dominant as compared with the case when 
torsional Hamiltonian is dominating Crefer to Figure V 1.1 oncK.<ti^ - 
This difference in the O dependence, for example, could 
distinguish the mechanism that is at work in spin 3x3 solid for 
the establ ishment of .spin-locked Cor quasi- equilibriunD state and 
its subsequent decay. However, we have not pursued this point 
further by means of appropriate experiments. Of course, to 
ascertain unambi guosl y, the particular interaction that provides a 
mechanism for attainment of spin-locked state, we have to do other 
kinds of measurements. For example, the combined investigations 
of resonance line-widths and relaxation time Cdi polar or 
qua.drupolar> ought to lead to more definite conclusions. However, 
it is to be noted that the asymmetry parameter tj, does not change 
the nature of the behavior of the with respect to the 

different parameters, as for example comparison of Figures 
VI-4.Cdl> and C«5 reveal. 
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Figui t? VI , 1 


Q>jasi stationary magnetization vs v'arious 
e>q[:jeri menial parameter s ; 

^ ^ iradians:!; F’aranieters ■. rj = 0.0, H** 

Li 

=0.0, A=0. 01 

(. b!) XM, vs A ; Parameters : t) = 0.0, H 

»=0. 0> 1.0 radian 

Cc5 M /H, vs T ; Paranteiers r-n = 0-0. H 

«q t ^ L. 

« 0 . 0 , A=0.001 

Cd> M /'hi vs S ; Paranieiers 177 ~0. £8> H 

©q t L 

5 Gauss, A~0, 0 


Cey H vs O I Paran^eiers : » = O.S 8 > H 

#q i 1- 

=0.0, A=0.01 

Cf> M vs & ; Parameters : >7 = 0.0, H 

eq C I- 

=0. 50. A=0. 0 

Ca5 M /M vs ^ ; Parameters .-77 = 0.0, H 

” •q i 

=0.0, A*0.01 1^^ = 45° 


M 

is local dipolar field. 

the phase of the pulses in the sequence s^ith respect to th 


preparatory pulse 
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(a) 0 = I- 0 radians 

(b) 0 = O*8 radians 
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VI. C. Experiwental Details and Results 

As stated in Chapter IV. an t k- - . 

iv, aii the experiments were performed 

at the Central Leather Research Institute Madras, India. Ilie 
instrumer^t used was a Bruker CXP-90 NMR Spectrometer with its do 
magnetic field switched off during the experiment. The 

sp"! n ""'1 ock X jri0 wd.s in ^<&n 6 i*P 2 Ll 

• CVI.36> 

A two-pulse program, in which the second pulse would occur between 

1 and 4005 counts during a loop was written and data were 

acquired at the end of each cycle, namely, at the ’echo tops*. 

The first pulse was set to be a *ny3' pulse in all the 

exper 1 nw&nts. The peak transmitter power employed was ss 

400-300 watts corresponding to ’irXS’ flip angle of 0 #.isecsC**Cl 

resonance in MaClO and HgCl 5. A standard wide-line MMR probe 

s t 

CIO mm diameter !> tuned to oscillate in the frequency range of 
20-33 MHz was employed. The phase of the pulses in the period 
could be controlled in steps by using the pulse-list 

freature of the system* s Z-17c pulser. The spectral width in this 
case also was i^Ar, where 2 t is the separation between the 0 
pulses in spin-locking sequence. The ambient temperature was 
maintained to be 24±i‘*C Cfor NaClO^ experiments? and 21±1°C Cfor 
HgCl ^ exper i nwnts? . 
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Th® samples chosen for the CW nr wm? i 

vne K/f( or MW experiments were finely 

crushed powder samples of NaClO^ and HgCL 

Sodium chlor.l. mor.ocrysl.ls h.v. . N.a .rr,r.gm.r.l of th. 
ior.s, dislortmd to .ccomod.l. th. CIO^ ,nior,s io such . way that 
the overall symmetry remains cubic £1733. The cubic unit cell 
contains four molecules. The four chlorine atoms, on one hand, 
and the four sodium atoms, on the other occupy equivalent 
lattice s. tes. Since oxygen is essentially nonmagnetic we negj.act 
its presence in the discussion of the properties of the spin 
system of sodium cholrate. The molecular bonds, joining chlorine 
atoms and the associated sodium ions, are parallel to the body 
diagonals of the unit cell. The efg at the location of sodium 
and chlorine nuclei show axial symmetry Cr)=05, resulting from 
Cl.l.O direction Cparallel to the molecular bonds> being a three 
fold symmetry axis. Sodium and the two isotopes *’ci and *^C1 of 
chlorine all have spin 3/^. As a consequence, only one 
quadrupolar transition frequency is associated with each of the 
chlorine and sodium isotopes. As pointed out earlier we have 


irradiated our sample only at **C1 C&> = 26. 


at a4± 1 O. 


The local field due to internal interactions in case of 
NaClO^ can be decomposed into the following components: 

Assuming direct dipolar couplings between nuclei on an 
undistorted rigid lattice, the dominant contribution C7.050 

comes 


from “Na-**Na interactions, ”Na-”ci interaction 


**Na-*’’Cl interactions contribute 0.742 G* . 
other interactions provide almost negligible contribution, the 


sontrlbute 3. 20S G 


Ml 


sum being - 11.10 G 



ail 


The HaCl, single crystal, has a teiramolecular orthorhombic 
unit cell tl?2K In this case also the local field has a 
negligible contribution from the homonuclear dipolar interaction. 
On the other hand, the angles in our theoretical computations were 

used for the dipolar interaction between two chlorine atoms in 

HgCl. 

z 

1 • e£ phase shift 

We generated *n^4* phase shift by simultaneously opening 
two quadrature phase channel outputs of the pulser— transmitter 
system of the CXP~©0 machine and scaling the corresponding pulse 
duration by a factor of l/'-jfe. 

2- The Ef ficiencv of Spin~-locklnq in the Present Expert ments 

In general > in mul tipi e-pulse experiments an effective field 
can be defined over a period Cconsidering a periodic perturbation 
or pulse trainD'. A fraction of the magnetization of the spin 
system gets locked and remains parallel to this effective field 
direction while the other fraction makes an angle with this 
effective field and evolves. The former component gives rise to a 
signal at the transmitter frequency corresponding to zero 
frequency in the interaction frame and is known as the pedestal. 
The latter gives rise to a signal with a non-zero frequency in the 
interaction frame- a frequency which in general reflects scaled 
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interactions that survive in ♦ k- t i 

urvive in the interaction franve under the 

multiple-pulse excitation. 

In the case of fresh. DriwHAr->w 4 , 

powaered NaClO^ sample* observations on 

th. oscilloscop. Lr.ees of th. FT of the qu.dr.lure delected 

time-domein sign.1, showed . sigoel to pedesl.l Intensity r.lio 

of approximately 2.-1. In other words* only 33>4 of the 

magnet i zation remains locked in the effective field direction 

under our experimental conditions, while - 67X of the 

magnetization evolves making an angle with respect to the 

effective field direction. It is thought that the above 

inefficiency in spin— locking arises essentially from poor rf 

inhomogeneity and/'or inadequate rf field strength which might be 

partially compensated for either by a suitable change in the 

preparation pulse flip angle, rf phase alternation during the 

spin-locking pulse train or both tl6] 

^ ■ The Short- and Long- ti me behavi or of the quasi egut 1 i br 1 urn 

magneti zation 

As far as the short-time behavior is concerened it is clear 
from all the spectra presented in Figure VI . 5 that the spin-locked 
st.ate is attained in spin 3/^ NC^ solids also. In the following 
we shall examine the initial oscillatory behavior as well as the 
subsequent long-time decay of the magnetization. 

As explained by Maricq C33.65.661 in case of spin 1^2 NMR and 
spin 1 NQf?. in case of spin 3/3 NQR also we observe an oscillatory 
behavior of the magnetization initially Csee Figure VI. 35. During 

of the signal the spins evolve effectively 


this initial portion 
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Fig’M'tf VI *3 ; Fr«9<iuency domain signal of the ”ci Ca> = 39.9333 
MHz'} r«rSonancie in NaClO^ powder ai room t/emperat^ure. 
CaD Fourier transform of the one pulse FID 
Parameters : A « Of *= 40^ 

Cb> Fourier transform of the ti'me domain signaJL 
fol lowing a spin-^iocking multiple-pulse sequence COW 
or MW> » showing clear evidence of line narrowing 
CParameters :A-0, r ^lOOpsec. & = 40^5 



Cl NQR in NaCI03 
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VI • 3 ; Th€i iaitiail poriion of the response to spin -locking 
mull tipi e*-pulse sequence COW or hPfD of ^C1 
29. 0232 >Wz> in NaClO^ powder. Parameters A - 
--0. 01 I 0 rr/'4 > r « 40 psecs. The signal has been 

convoluted with an exponential function to reveal 
the initial oscillations clearly. 



2l6 



Figure VI-3 



Figure VX*4 : Illustration of CaZ) the predicted evolution under 


effective Hamiltonian vs Cb-> the 
time -development of the magnetization 
time-periodic system t66]. 


obser ved 
for a 





Figure VI.4(b) 
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Figuro VI . 5 ; Tim® -domain response to OW or MV? raul tipi e-pul se 

‘.«H.juonce of the ®'’ci signal in NaCIO Cco 

9 

atJ. (123:) and HgCl^ ~ 28.072 MHz^ 

< a) NaClO l A «=“S. 0 KHz ; t = 100 psec ; & = 40°. 

II 

( b.) NaClO ;A ■+5. 0 KHz ; t = 100 psec ; B = 90°. 
t c) HgCl^jA »-10. 0 KHz ; T = 40 psec ; © = 45°. 

Cd> H<3C1^:A «-13, 070 KHz ; t = 40 psec ; d = 55°. 

Ce> HgCl A =■*■4 KHz ; t = 40 psec ; & = 45°. 

Cf> HgCl :A = 0. 0 KHz ; t = 40 psec ; ^ = 90°. 

M 

Cg.) HgCl j A » 0. 0 KHz ; t = 40 psec ; ^ = 45°. 

Ch’J HgCl A =-0.2 KHz ; t =41.955 psec ; O = 29.10°. 

CiD HgCl ‘A = 4.0 KHz ; t = 40 psec ; ^ = 63. 60°. =45°. 

CJ> HgCl A =-14 KHz ; T = 40 psec ; ^ = 90 . 

Ck> HgCl^jA =-^2.0 KHz ; r = 41.955 psec ; © = 29.10°. 

Cn HgCl^lA =-10.0 KHz ; T = 40 psec ; & = 90i?. 
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under the action of . This time period features a wide 

spectrum of oscillations of the system magnetization. This 

occurs because the precession rate of each spin is determined by 
the local magnetic field and electric field gradient surroun g 
it, in addition to the deviation from resonance. These 

oscillations arise from the component of the magnetization vector 
initially perpendicular to the effective rf field axis. There is. 
in addition a constant contribution to the magnetization that 
arises from the parallel component of the initial magnetization. 
After a time of approximately T^. the transverse component decays 
to zero and quasi stationary state is reached. The time required 
to. arive at this quasi stationary state is determined by strength 
of the dominant internal interaction Ctorsional oscillations in 
case of NOR of spin 3^ solids:) and the flip angle employed for 
the pulses in the spin-locking multiple-pulse sequence. It 

depends only indirectly on the interval between pulses, in the 
sence that if x is too large, does not represent the 

effective Hamiltonian even for short times. 

A prior knowledge about the long-time behavior of the system 
of an ensemble of nuclei is important even when we restrict the 
problem to torsional oscillation Hamiltonian only. A very small 
imprecision in the initial conditions, i - e. the initial values of 
the parameters. can accumulate a substantial error in the 
subsequent evolution of the magnetization. A study of the 

long-time behavior from the point of view of the stability of the 
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equilibrium magnetization also belcomes important. 

When the NOP response of spin 3^ solids to spin locking 
multiple-pulse sequence is sampled only at observation v 
once every period, the system Is supposed to behave as shown in 
Figure VI.4Ca3. In actuality, however. Figure VI.4<;b> represents 
the true behavior of the s;rstem; that is to say, the 

quasi equilibrium does not last long and decays over a 1 g 
time scale to an equilibrium state 165.661. In our experiments 
also, we have observed this behavior of the magnetisation in the 
long-time regime. We attribute this to a case of looked 

quadrupolar relaxation, because of which, as seen in the Figures 
VI. seal and Cbl . the magnetisation does not remain in the 
quasiequilibrium state under the spin-locked condition for long. 
But the important point to be noted is that the magnetisation will 
have a definite value for periods much longer than ST^ CT^ for 
NaClO,. for example, is ~ 45 msecs at room temperature but we see 
signal to be persisting for as long as 300 msec. Figure Vl.SCaO 
and fbl. Although, we have not carried out a theoretical 
investigation for 3 > 1 radian for the flip angles of the pulses 

t Vis* OW or MW experiment s we have also 
in the sequence for the uw or nir t- 

presented our experimental results for this case an 
conclusions are very similar to what has been said about 3 < 1 


radian case. 
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VI. D. Summary 

In this Chapter we have derived theoreticaly the expression 
for the equilibrium magnetization in the case of spin solids 
when irradiated with OW or MW multiple-pulse sequence. Both the 
theoretical curves and experimental spectra give a strong evidence 
of attainment of a spin-locked state in spin 3X3 NQR case also, as 
in the case of spin 1X3 NMR and spin 1 NQR but the different is 
that in case of spin 3X3 NOR it is a locked quadrupolar relaxation 
and its mechanism is dicatated by the torsional osciallations of 
the units containing the guadrupolar nuclei under consideration. 
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SUMMARY OF FINDINGS AND FURTHER PROSPECTS 


Although as pointed out right at the beginning of this th -sis 
spins 1/^, 1 and 3^ have distinct and special fea'tures of energy 
level splittings and transitions, it is noteworthy that the 
dynamics and equilibrium magnetization properties under different 
schemes of pulsed rf p>erturbati on are very analogous. Indeed, the 
generaliasation seems possible that, in an ensemble of nuclei with 
any spin I, as long as a single frequency radiation connects only 
a pair of levels ^degenerate or non-degenerate> , the behavior will 
be very similar to what has been demonstrated in the past in spin 
l/^'S NMR, spin 1 NQR and, in the present thesis, for the case 
of spin 3^ NC8?:. Theoretically we can discuss this similarity by 
choosing the appropriate representation and operator formalism to 
describe the evolution of the system. 

In this thesis this fact has been demonstrated for spin 
NOR by going to diagonalized representation and choosing the 

single transition operator C JaTOD formalism for the description of 
the Hamiltonians. An important result is the development of the 


special operators L, , and in Chapter II. We have set forth 
the commutation and other properties of L,M^ and ; I = 1 , 2 , 



X 
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respect! vely» used in spin Ixa NMR terminology. Once w® express 
the Hamiltonian of system in diagonalized representation and STO 
formalism the dynamics and equilibrium properties of spin 3/^ 
solids show striking similarity to what has been shown in the past 
in case of spin NMR and spin 1 

The difference between the observed response now relies 
predomi nantl y on only the magnitude of the various Hamiltonians in 
the interaction representation that governs the response of the 
system. In case of spin NQR» for example the broadening of 

the line, the spin-locking and the quasi -equilibrium state have 
been shown to be dictated by the torsional Hamiltonian^ in 
contrast to the preponderance of dipolar interaction in case of 
spin 1/^3 NMR. 

We shall now outline certian prospects for future 
i Dvesti gati ons : 

Ci5 The theory for mul ti pie— pul se perturbation of spin 3X3 hiQ^ 
system has been worked out only for the lindting case of 6-puises. 
The correction due to finite pulse width needs to be worked out 
CiiD The theory for spin-locking has been worked out for the flip 
angle of the pulses in the .sequence & <1 radian only. Extending 
the theory to investigate the results for &>1 radian would also 
be fruitful. The spin-locking experiment in case of & >1 radian 
has also been presented in this thesis* so we expect that the 
theory should give us a result that will be similar to what has 
been observed in the case ^<1 radian. 

Ciii> Going to submultiple regions of periods in the 
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per iodic -perturbation schemes to probe the behaviour of the syst&m 
deeper seems to be another interesting possibility. 

CivD Use of small magnetic fields to perturb the NQR energy levels 
and separate the transitions from each other, and then performance 
of a multiple-pulse experiments are also required. Tliis will help 
us understand how the various modes of relaxation in spin 3/^ NQR 
are affected during multiple-pulse sequences. 

Cv!> A study of the temper at ure-dependance of the response of 
multiple-pulse sequence to separate out information about various 
relaxation mechanisms and how they are affected during a 
mul iiple-pul se sequence perturbation will be an interesting study, 
and finally 

Cvi“J) Generalization to higher spin systems seems to be an 
interesting prospect for future. 
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APPENDIX A 


C 

C Program SQMGNT.( Scaled Quasiequilibriiun MaGNeTization. ) 

C This program calculates the quasiequilibrium magnetization 

C for a spin 3/2 NQR powder specimen when irradiated with an 

C Ostroff -Waugh (OW) Multiple-pulse sequence. The quasiequilibrium 

C magnetization is scaled down by dividing it by initial 

C magnetization so that its value is always between 0.0 and 1.0. 

C The internal interactions that has been taken into account In 

C the theory are the homonuclear dipolar coupling between two 

C nuclei and the distribution in electric field gradient (efg) 

C in the vicinity of a resonant nucleus. 

C The variable names are almost directly in correspondence wxth 

C the symbols used in expressions of Chapter Six in the Thesis. ^ 

C NOTE: A digit (1-9) in coluam 9 means continuation line. Also ; 

C00000000111111111122222222223333333333444444444455555555556666666666777 
C23456789012345678901 234567890123466789012345678901234567890123456789012 


C 

C 

C 

C 

C 

C 


PROGRAM SQMGNT 

COMPLEX YP1,YP2,YN1,YN2,YP3,RMBAR,RNBAR,RBBAR 
COMPLEX DCl,DC2,DC3,DC4,DC5,DC6,r>C7 

COMPLEX HD, DKLP, HI, H3,Q, RATIO, CON3 , / i nn ^ 

DIMENSION XXX(IOO) , YYY(IOO) ,SINTL{100) ,COSTL(100) ,SINFL{100) , 

1 COSFL(IOO) TITLE(7) 

COMMON/BLKl /AX , AY , AZ , B , C , D , ETA , TAU , DELWQ , THETAO , R , THETA , HL , RRKK 
COMMON/BLK2/XIC , XIS , XIPl . XINl , XIP2 , XIN2 
(X)MM0N/BLK3/XP1(4,4) ,XP2(4,4) ,XN1(4,4) ,XN2(4,4) , 

1 XN3(4,4),ZN1(4,4),ZN2(4,4),ZP2(4,4),AIDEN(4,4) 
(X)MKK)N/BLK4/YP1(4,4) ,YP2{4,4) ,YN1(4,4) ,YN2(4,4) ,YP3(4,4) 

C0MM0N/BLK5 /DC 1 ( 1 6 , 1 6 ) , DC 2 ( 1 6 , 1 6 ) , DC 3 ( 1 6 , 1 6 ) , 

1 DC4(16,16),DC5(16,16),DC6(16,16),DC7(16,16) 

CC»filON/BLK7/HD(16, 16) ,DKLP{16, 16) ,H1(16, 16) ,H3{16, 16) ,Q(16, 16) 
C0MM0N/BLK8 /RK , RL , RM , RN , RP , RQ 
C0MM0N/BLK8 1 /RIK , R1 L , RIM , RIN , RIP , RIQ 

<X)MMON/BLK82/R2K,R2L,R2M,R2N,R2P,R2Q . r. . 

COMMON/BLK9/RBBAR(16, 16 ) , RLBAR{ 16, 16) , RMBAR{ 16 , 16) ,RNBAR(16, 16) 
COMMON/BLK10/HDSMSQ,HLSQ,RSQ8 ^ 

OPEN(UNIT=20 , DEVICE= ’ DSK’ , FILE= ’ S^GNT. DAT .ACCESS- SEQIN ) 

OPEN UNIT=24 DEVICEr^DSK’ FILE=’ SQMGNT. OUT' , ACCESS =’SEQOUT' ) 
OPEN{aNIT=6,DEVICE='DSK* , FILE=' SQMGNT . PLT’ . ACCESS =’SEQOUT' ) 


WQ is resonance frequency of the quadrupolar nucleus (in MHz). 
(1+DELTA)*WQ is offset frequency(in MHz). . 

THETA is flip angle of pulses in the sequence (in Radians). 
THETAO is flip angle of the preparatory pulse (in Radians). 
2*TAD is time duration between two pulses in the 
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sequence { in microsecs ) . 

ETA is asymmetry parameter (0.0 < ETA < 1.0), 

THETAD and FID are angles that relate QPAS of a crystallite 
to the Dipolar Principal Axis System (DPAS) (in Radians). 

HL is local dipolar field (in Gauss), 

FIl is phase of the first i.e., preparatory pulse(in Radians). 
FI2 is phase of pulses in the sequence with respect to that of 
the preparatory pulse i.e., first pulse (in Radians). 

FIS is phase of the detector with respect to first pulse. 

IVARI indicates variation parameter: 

IVARI=1 is for THETA variation; 

IVARI=2 is for TAU variation; 

IVARI =3 is for DELTA variation, 

VINI=initial value of varying parameter. 

VFIN=final value of varyii'vg parameter, 

IVN=n\imber of variation steps, 

NT=number of THETA(LAB)S’ ,NF=number of FI (LAB)S’ , used in powder 
averaging. Typically NT=25,NF=25. 


READ(20,100)(TITLE(I),I=1,7) 
100 FORMAT (7A4) 

READ(20,*)WQ 

READ( 20 , * ) THETAO , ETA 

READ ( 20 ,*) HL , THETAD , FID 

READ(20,*)FI1,FI2,F13 

READ ( 20 ,*) THETA , TAU , DELTA 

READ( 20,*) IVARI 

READ ( 20 , *) VINI , VEIN , IVN 

READ(20,*)NT.NF 

CLOSE(DNIT=20) 


WRITE (24, 110) TITLE 
WRITE (6, 110) TITLE 

110 F0RMAT( IX, ’OUTPUT FROM SQMGNT’ /IX, 7A4) 

WRITE (24,1 20 )WQ 
WRITE(6,120)WQ 

120 FORMAT( IX, ’RESONANCE FREQUENCY^’ ,E13 . 6 , ’MHZ ’ ) 

WRITE (24, 130) THETAO, ETA 
WRI TE ( 6 , 1 30 ) THETAO , ETA 

1 30 FORMAT ( IX , ’ FLIP ANGLE OF PREPARATORY PULSE= ’ . F9 . 6 , ’ RADIANS ’ 

1 , 2X , ’ ASYMMETRY PARAMETER= ’ , F5 . 3 ) 

WRI TE ( 24 , 140 ) HL , THETAD , FID 
WRI TE ( 6 , 1 40 ) HL , THETAD , F ID 

1 40 FORMAT ( IX , ’ LOCAL DIPOLAR FIELD= ’ , F6 . 3 , ’ GAUSS ’ /IX , 

1 ’DIPOLAR THETA=’ ,F9. 6, ’RADIANS ’/IX, 

2’DIPOLAR FI=’ ,F9. 6, ’RADIANS’ ) 

WRITE (24,1 50 )FI1 
WRITE(6,150)FI1 

150 FORMAT( IX, ’PHASE OF THE FIRST i . e ., PREPARATORY PULSE=’,F9.6 
1 , ’RADIANS’) 
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WBITE{24,160)FI2 

WRITE(6,160)FI2 

160 FORMAT( IX, ’PHASE OF PULSES IN THE SEQUENCE W.R.T. FIRST PULSE =’ 
1 ,F9. 6, ’RADIANS’ ) 

WRITE(24,170)FI3 
WRITE (6, 170 )FI 3 

170 FORMAT( IX, ’PHASE OF THE DETECTOR W.R.T. FIRST PULSE= ’ , F9 . 6 , 

1 ’ RADIANS ’ ) 

GO TO(180,220,250),IVARI 

180 WRITE (24, 190) TAD 

WRITE(6,190)TAU 

190 FORMATCIX, ’TIME DURATION BETWEEN TWO PULSES IN THE SEQUENCE=’ , 

1 E13. 6, ’MICROSECS’ ) 

WRITE (24, 200) DELTA 
WRITE (6, 200) DELTA 

200 FORMAT (IX, ’OFFSET PARAMETER =’ ,E1 3. 6) 

WRITE( 24 , 210 ) VINI , VFIN 
WRI TE ( 6 , 2 1 0 ) VINI , VFIN 

210 FORMAT( IX, ’INITIAL AND FINAL VALUES OF THETA= ’ , F9 . 6 , IX, ’ , ’ , IX, 

1 F9.6, ’ (RADIANS)’ ) 

GO TO 270 

220 WRITE (24, 200 )DELTA 

WRITE (6, 200) DELTA 
WRITE (24, 230) THETA 
WRITE (6, 230) THETA 

230 FORMATdX, ’FLIP ANGLE OF THE SEQUENCE PDLSES= ’, F9 . 6 ,’ RADIANS ’ ) 
WRITE( 24, 240) VINI, VFIN 
WRI TE ( 6 , 240 ) VINI , VFIN 

240 FORMATdX, ’INITIAL AND FINAL VALUES OF TAU=’ ,E13.6,2X,E13.6, 

1 ’MICROSEC’) 

GO TO 270 

260 WRITE(24,230)THETA 
WRITE (6, 230) THETA 
WRITE (24,1 90 )TAU 
WRITE(6,190)TAU 
WRITE ( 24 , 260 ) VINI , VFIN 
WRITE (6, 260) VINI, VFIN 

260 FORMATdX, ’INITIAL AND FINAL VALUES OF OFFSET PARAMETER =’ , 

1 E13.6,2X,E13.6) 

270 CONTINUE 

WRITE ( 24, 280) IVN 
WRITE(6,280)IVN 

280 FORMAT (IX , ’ NO . OF STEPS IN VARIATION= ’ , 12 ) 

VDEL=(VFIN-VINI)/IVN 
WRITE (24, 290 )VDEL 
WRITE(6,290)VDEL 

290 FORMATdX, ’INCREMENT IN VARIATIONS =’,E13.6) 

WRITE (24, 300 )NT,NF 
WRITE(6,300)NT,NF 

300 FORMATdX, ’NO. OF THETA (LAB) VALUES=’ , I4/1X, ’NO. OF FI (LAB) 
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1VALUES=M4) 

C 

HL=HL>«'2621.4614E-06 

HLSQ=HL*HL 

R3=SQRT(3.0) 

PI=4.0*ATAN(1.0) 

WQ=WQ*2.0*PI 

XI=ATAN(ETA/R3) 

XIC=COS(XI) 

XIS=SIN(XI) 

XIP1=R3*XIC+XIS 

XIN1=R3=J=XIC-XIS 

XIP2=R3*XIS+XIC 

XIN2=R3*XIS-XIC 

RRKK=WQ*SQRT ( 1 . O+ETA^ETA/3 . 0 ) 

RRKK6=RRKK/6 . 0 

DELWQ=DELTA*RRKK6 

CALL INPMAT 

STD=SIN(THETAD) 

ST2D=SIN ( 2 . 0*THETAD ) 

STDD=STD*STD 

CTD=COS{THETAI>) 

CTDD=CTD*CTD 

R6=1.0/SQRT{6.0) 

R5=SQRT(1.6) 

SFD=SIN(FID) 

CFD=COS(FID) 

SF2D=SIN{2.0*FID) 
CF2D=COS(2.0*FID) 
E1=R5*STDD*CF2D 
E2 =R6* (3.0 *CTDD- 1.0) 

AX=E1-E2 

AY=-E1-E2 

AZ=2.0*E2 

B=R5’»:STDD*SF2D 

C=R5*ST2D*CFD 

D=R6’^ST2D*SFD 

CALL HDMAT 


Powder averaging preliminaries: 

A Monte Carlo procedure called ’Systematic Sampling’ is used. 
Reference: l)S.C.Sivasubramanian,Ph.D. Thesis , Indian Institute 
of Technology Kanpur, India (1989) . 

2)P.Raghunathan and S.C.Sivasubramanian,Proc. Indian Acad. Sci. 
(Chem. Sci. ), 96(1986)565. 


TWOPI=PI+PI 

NTNF=NT*NF 

CPHAS1=C0S(FI1) 

SPHAS1=SIN(FI1) 
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CPHAS2=COS{FI2) 

SPHAS2=SIN(FI2) 

CPHAS3=COS(FI3) 

SPHAS3=SIN{FI3) 

DO 310 1=1, NT 

COSTL(I)=((1,0-2.0*I)/NT)+1.0 
310 SINTL(r)=SQRT{{1.0-COSTL{I)**2)) 
DO 320 1=1, NF 
RND= ( ( I -0 . 6 ) /NF ) *TWOPI 
SINFL(I)=SIN(RND) 

320 COSFL(I)=COS(RND) 

XXX(1)=VINI 
YYY(1)=0.0 
DO 330 1=2, IVN 
XXX ( I ) =XXX ( I - 1 > + VDEL 
330 YYY(I)=0.0 


Powder averaging loop starts : 

TL and FL are the angles that relate QPAS to the detection axis, 
namely, the X-axis of laboratory frame. (Crystal Orientation 
parameters ) . 

MEQMI is called NT*NF times for each powder averaging. 


DO 370 1=1, NT 
DO 360 J=1,NF 

E1=-0.6*SINTL(I)*SINFL(J)*XIP1 

RK=E1»CPHAS1 

R1K=E1*CPHAS2 

R2K=E1*CPHAS3 

RL=E1*SPHAS1 

R1L=E1*SPHAS2 

R2L=E1*SPHAS3 

E2=0 . 5*COSTL{ I )*( -XINl ) 

RM=E2*CPHAS1 

RlM=E2=t^CPHAS2 

R2M=E2*CPHAS3 

RN=-E2*SPHAS1 

R1N=-E2=^SPHAS2 

R2N=-E2*SPHAS3 

E3=SINTL(I)*COSFL(J)*XIS 

RP=E3*CPHAS1 

R1P=E3*CPHAS2 

R2P=E3*CPHAS3 

RQ=E3*SPHAS1 

R1Q=E3*SPHAS2 

R2Q=E3*SPHAS3 

R=SQRT( (IiK*RK+RL*RL+RM*RM-ARN*RN+EP*RP+RQ*RQ) ) 

RSQ8=R*R*8.0 

CALL HDPROJ 

CALL LMNBAR 
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C ! ! ! ! ! JGENERATE Q! i I i ! ! 

CONI =SIN ( THETA0*R) /R 
CON2=COS(THETAO=t‘R) 

CON3=CON1*(0.0,1.0) 

DO 340 ID=1,16 
DO 340 JD=1,16 

340 Q( ID, JD) =RRKK*(CON2=kRLBAR( ID, JD) -CON3*RMBAR( ID , JD) ) 


Parameter variation loop starts: 

Each variation (increment by VDEL) requires one powder 
averaging. However , we have arranged this loop in such a way that 
powder averaging is done simultaneously for all parameter 
(THETA, TAD or DELTA) variations. 


DO 350 11=1, IVN 
TYPE*, II, I, J 

IF( IVARI .EQ. 1 ) THETA=XXX( II ) 

IF( I VARI . EQ . 2 ) TAD=XXX( I I ) 

IF( IVARI .EQ. 3 )DELTA=XXX( II ) 

IF( IVARI . EQ. 3 )DELWQ=DELTA*RRKK6 
CALL MEQMI (RATIO) 

YYY( I I ) =YYY( I I ) +REAL ( RATIO ) 

350 CONTINUE 
360 CONTINUE 
370 CONTINUE 

WRITE (24, 380) 

380 FORMAT (IX, ’VARIATION PARAMETER VS. EQ. MAGNETIZATION’) 

DO 390 11=1, IVN 
YYY( II ) =YYY( II ) /NTNF 
WRITE( 24 , * )XXX( II ) , YYY( II ) 

390 CONTINUE 


CALL PLOT(XXX,YYY,IVN) 
CLOSE {DNIT=6) 

CLOSE (0NIT=24) 

STOP 

END 


SUBROUTINE INPMAT 

Generates input matrices which are Single Transition Operators 
COMPLEX YP1,YP2,YN1,YN2,YP3 

COMMON/BLK3/XPl(4,4),XP2(4,4),XHl(4,4),XN2(4,4), 

1 XN3(4,4),ZN1(4,4),ZN2(4,4),ZP2(4,4),AIDEN(4,4) 
COMMON/BLK4/YPl(4,4),YP2(4,4),YNl(4,4),yN2(4,4),YP3(4,4) 

DO 100 1=1,4 

DO 100 J=l,4 , „ „ 

XPl ( I , J ) =0 . 0 ; XP2 (I,J)=0.0;XN1(I,J)=0.0; XN2 ( I , J ) =0 . 0 
XN3 ( I , J ) =0 . 0 ; ZP2 ( I , J ) =0 . 0 ; ZNl ( I , J ) =0 . 0 ; ZN2 ( I , J ) =0 . 0 
YPl (I,J)=(0.0,0.0); YP2 (I,J)=(0.0,0.0);YNl(I,j)=(0.0,0.0) 
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YN2(I, J)=(0 

CONTINUE 

XP1(1,2)=0. 

XP2(1,4)=0. 

XN1(1,2)=0. 

XN2{1,4)=0. 

XN3(1,3)=0. 

yPl(l,2)=(0 

YP1(3,4)=(0 

YP2(1,4)=(0 

YP2(2,3)=(0 

yNl(l,2)=(0 

YN1(3,4)=(0 

YN2{1,4)={0 

YN2(2. 3)=(0 

YP3{1,3)={0 

YP3(2,4)=(0 

ZN2{1,1)=0. 

ZP2(1,1)=0. 

ZN1{1,1)=0. 

AIDEN(1,1)= 

RETURN 

END 


.0,0.0); YP3 (I,J)={0.0,0.0); AIDEN( I , J ) =0 . 0 


6;XP1(2,1)=0.5: 

5;XP2(4,1)=0.5 

5;XN1(2,1)=0.5 

5;XN2(4,1)=0.5 

5;XN3(3,1)=0.5 

.0,-0.5);YPl(2 

.0,-0.6);YPl(4 

.0,-0.5);YP2{4 

.0,-0.5);YP2(3 

.0.-0.5);YNl(2 

.0,0.5);YN1(4, 

.0,-0.5);YN2(4 

.0,0.5);YN2(3, 

.0,-0.5);YP3(3 

.0,-0.5);YP3(4 

6;ZN2{4,4)=-0. 

5;ZP2(4,4)=-0. 

5;ZNl{2,2)=-0. 

:1.0;AIDEN(2,2) 


;XP1(3,4)=0.5;XP1(4,3)=0.5 

;XP2(2,3)=0.5;XP2{3,2)=0.5 

;XNl(3,4)=-0.5;XNl(4,3)=-0.5 

;XN2{2,3)=-0.5;XN2(3,2)=-0.5 

;XN3(2,4)=-0.5;XN3(4,2)=-0.6 

, 1 )=( 0 . 0 , 0 . 6 ) 

,3)=(0.0,0.5) 

,1)={0.0.0.5) 

,2)=(0.0,0.5) 

,1)=(0.0,0.5) 

3)=(0.0,-0.5) 

1)=(0.0,0.5) 

2)=(0.0,-0.6) 
l,l) = {0.0,0.5) 

:,2) = (0.0,0.5) 

5;ZN2(2,2)=-0.5;ZN2(3,3)=0.5 
5;ZP2(2,2)=0.5;ZP2(3,3)=-0.5 
5; ZN1{3, 3)=-0. 5; ZN1(4,4)=0. 5 

i = 1 . 0 ; A I DEN ( 3 , 3 ) = 1 . 0 ; A I DEN ( 4 , 4 ) = 1 . 0 


SUBROUTINE HDMAT . 

Generates the homonuclear Dipolar Hamiltonian. 

CO^ON^LKl /S!Ay?Az!B , C , D , ETA , TAU , DELWa , THETAO , R , THETA . HL . ERKK 
COMMON/BLK2/IIC , XIS , XIPl , XINl , XIP2 . XIN2 

COMMON/BLK3/XP1 (4 , 4 ) ,XP2(4 , 4 ) ,XN1 (4 , 4 ) ,1^2(4 , 4 ) , 

1 XN3(4 4).ZN1(4,4),ZN2(4,4),ZP2(4,4),AIDEN(4,4) 

COMMON/BLK4 /YPl ( 4 , 4 ) ,' YP2 ( 4 , 4 ) , YNl ( 4 , 4 ) , YN2 ( 4 , 4 ) , YP3 ( 4 , 4 ) 

COMMON/BLK5/DC1 (16,16),DC2{16,16), D^ { 1 ^ J > ' 

1 DC4(16, 16) ,DC5(16, 16 ) ,DC6(16, 16) ,DC7( 16, 16) T»c/ie ic\ 

COMMON/BLK6/Dl(16, 16),D2(16, 16) ,D3{16, 16) ,D4{16, 16) ,D5(16, 16) , 

CO^ON^LK7/HD( 16 , 16 ) , DKLP( 16 , 16 ) , HI ( 16 , 16 ) , H3 ( 16 , 16 ) , Q( 16 , 16 ) 
COMMON/BLKIO/HDSMSQ, HLSQ , RSQ8 
COMMON/BLKl 1 /HDD (16,16) 

DO 100 1=1,16 
DO 100 J=l,16 
HD(I,J)=(0. 0,0.0) 

CALL DPSMRC(XP1,XP1,YN1,YN1,DC2,1.0) 

CALL DRPRRR{XP2.XP2,D1) 

CALL DPSMRR(XP2,XN2,XN2,XP2,D4,1.0) 

CALL DRPRRR{XN2,XN2,D2) 

E1=0.5*XIP1*XIP1 

E2=XIN2*XIN2 
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110 


120 


130 


140 


DO 110 1=1,16 

DO 110 J=l,16 ,,, 

HD(I,J)=HD(I,J) + (E1*DC2(I,J) + (D1(I,J)+D4(I,J)*XIN2+E2*D2(I, J) ) 

1)*AX 

CONTINUE 

CALL DPSMRC(XN1,XN1,YP1,YP1,DC1,1.0) 

CALL DRPRCC(YN2,YN2,DC2) 

CALL DPSMCC(YP2,YN2,YN2,YP2,DC3,1.0) 

CALL DRPRCC { YP2 , YP2 , DC4 ) 

ElrO.S’t'XINl^XINl 
E2=XIP2*XIP2 
DO 120 1=1,16 
IXD 1 2 0 J **’ 1 16 

HD(I, J)=HD(I,J) + (E1*I>C1(I. J)+I>C2(I,J)-XIP2*DC3(I,J)+E2*DC4(I, J) 

DM 

CONTINUE 

CALL DRPRRR(ZN2,ZN2,D1) 

CALL DPSMRR{ ZN2 , ZP2 , ZP2 , ZN2 , D2 , 1 . 0 ) 

CALL DRPRRR(ZP2,ZP2,D3) 

CALL DPSMRC{XN3 , XN3 , YP3 , YP3 , DCl ,1.0) 


E2=2.0*XIC 

E3=E2*E2 

E4=2.0*XIS*XIS 


DO 130 1=1,16 

ro ( I , J ) =HD (I,J) + ((D1(I,J) +E2>^D2 ( I , J ) +E3*D3 ( I , J ) ) +E4*DC1 ( I , J ) ) *AZ 


CONTINUE 

CALL DPSMCR(YN1,XN1,XP1,YP1,DC5,-1.0) 
CALL DPSMCR{YP1,XP1,XN1,YN1,DC6,-1.0) 
CALL MTADCC { DC5 , DC6 , DCl , 1 . 0 ) 

CALL DPSMCR( YN2 , XP2 , XP2 , YN2 , DC2 , 1 . 0 ) 
CALL DPSMCR(YN2,XN2,XN2,YN2,DC3,1.0) 
CALL DPSMCR{YP2,XP2,XP2,YP2,DC4,1.0) 
CALL DPSMCR(YP2,XN2,XN2,YP2,DC5,1.0) 


E1=0. 5*XIP1*XIN1 
E3=XIP2*XIN2 
DO 140 1=1,16 
DO 140 J=l,16 

HD(I, J)=HD(I,J)+(E1*DC1(I,J)-DC2(I,J) 
1DC4{I, J)+E3*DC5(I, J))*B 

CONTINUE 

CALL DPSMRR(XP1,XN3,XN3,XP1,D3,1.0) 
CALL DPSMCC{ YNl , YP3 , YP3 , YNl , DC4 , 1 . 0 ) 
CALL MTADRC(D3,DC4,DC1,1.0) 

CALL DPSMRR(XP2,ZN2,ZN2,XP2,D1,1.0) 
CALL DPSMRR(XP2,ZP2,ZP2,XP2,D2,1.0) 
CALL DPSMRR(XN2,ZN2,ZN2,XN2,D3,1.0) 
CALL DPSMER(XN2,ZP2,ZP2,XN2,D4,1.0) 


-XIN2*DC3 ( I , J ) +XIP2* 


E1=-XIP1*XIS 

E3=E2*XIN2 
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150 


160 


170 


180 


C 

C 

C 

C 


100 


C 

C 

C 

C 

C 


DO 150 1=1,16 

DO 150 J=l,16 T^^ 

HD(I,J)=HD(I, J)+(E1*DC1(I, J)+(D1{I, J)+E2*D2(I,J)+XIN2*D3(I, J)+ 

1E3*D4{I, J)))*C 

CONTINUE 

CALL DPSMCR(YP1,XN3,XN1,YP3,DC6,-1.0) 

CALL DPSMCR(YP3,XN1,XN3,YP1,DC7,-1.0) 

CALL MTADCC ( DC6 , DC7 , DCl , 1 . 0 ) 

CALL DPSMCR{ YN2 , ZN2 , ZN2 , YN2 , DC2 , 1 . 0 ) 

CALL DPSMCR ( YN2 , ZP2 , ZP2 , YN2 , DC3 , 1 . 0 ) 

CALL DPSMCR{ YP2 , ZP2 , ZP2 , YP2 , DC4 , 1 . 0 ) 

CALL DPSMCR( YP2 , ZN2 , ZN2 , YP2 , DC5 ,1.0) 

E1=-XIN1*XIS 
E3=E2*XIP2 
DO 160 1=1,16 
DO 160 J~1 IS 

HD ( I , J ) =HD( I , J ) + (E1*DC1 ( I , J ) -DC2 ( I , J ) -E2*DC3 ( I , J ) +E3*DC4 ( I , J ) 

1 +XIP2^DC5{I, J))*D 
CONTINUE 

CALL MTMLCC(HD,HD,DC1) 

HDSMSQ=0 . 0 
DO 170 1=1,16 

HDSMSQ=HDSMSQ+REAL(DC1 (1,1)) 

DO 180 1=1,16 
DO 180 J=l,16 
HDD(I,J)=HD(I,J) 

RETURN 

END 


SUBROUTINE HDPROJ 

This subroutine projects the dipolar operator on the quadrupolar 
interaction frame (QIF) . Details are given in Chapter Six. 


COMPLEX HD,DKLP,H1,H3,Q,HDD 

COMMON/BLK7/HD(16, 16) ,DKLP(16,16) ,H1 ( 16 , 16 ) ,H3 ( 16 , 16 ) ,Q( 16 , 16 ) 
COMMON/BLKIO/HDSMSQ , HLSQ , RSQ8 
COMMON/BLKl 1 /HDD (16,16) 

WD=SQRT ( ( HLSQ*RSQ8 ) /HDSMSQ ) 


DO 100 1=1,16 

DO 100 J=l,16 

HD(I,J)=HDD(I,J)*WD 

RETURN 

END 


SUBROUTINE LMNBAR . ^ 

Grenerates matrices LBAR,tffiAR,NBAR and BBAR. vt 4 

BEAR is the operator that corresponds to the observable in an 
NOR experiment i.e.,the magnetization ^^^ponsible for si^al^ 


COMPLEX YPl , YP2 , YNl , YN2 , YP3 , RMBAR , RNBAR , RBBAR 



AlO 


C 

C 


COMPLEX XC1(4,4),XC2(4,4),XC3(4.4) _ 

COMMON/BLKl /AX , AY , AZ , B , C , D , ETA , TAU , DELWQ , THETAO , R , THETA , HL , RRKK 
COMMON/BLK3/XPl{4,4),XP2(4,4),XNl{4,4),XN2(4,4), 

1 XN3(4,4),ZN1(4,4),ZN2(4,4),ZP2(4,4),AIDEN(4,4) 
COMMON/BLK4/YPl{4,4),YP2(4,4),yNl(4.4),YN2{4,4),YP3{4,4) 
COMMON/BLK6/Dl(16,16) ,D2( 16,16 ),D3( 16,16 ) ,D4( 16,16 ) ,D5( 16,16 ) , 

1D6{16,16) 

COMMON/BLK8/RK , RL , RM , RN , RP , RQ 
COMMON/BLK81/R1K , RIL , RIM , RIN , RIP, RIQ 
COMMON/BLK82/R2K,R2L,R2M,R2N,R2P,R2Q 

COMMON/BLK9/RBBAR{16, 16),RLBAR{16,16),RMBAR(16, 16) ,RNBAR{16, 16) 
CALL DPSMRR(ZN1, AIDER, AIDER, ZN1,RLBAR, 1.0) 

DO 100 1=1,4 

^l(I^jt=(0%,1.0)*(-RK*YNl(I,J)+RL*XPl(I,J)-»-RM*XNl(I,J) 

1 -RN^f^yPl { I , J ) -RP»YP3 { I , J ) +RQ*XN3 ( I , J ) ) 
XC2(I,J)=R1K*XP1(I,J)+R1L*YR1(I,J)+R1M*YP1(I,J)+ 

1 R1N*XN1(I, J)+R1P*XN3(I, J)+R1QW3(I, J) 

XC3(I,J)=R2K*XP1(I, J)+R2L*YR1(I,J)+R2M*YP1(I,J)+ 

1 R2N*XN1 ( I , J ) +R2P=»;XN3 { I , J ) +R2Q*YP3 ( I , J ) 

100 CORTIRUE ^ 

CALL DPSMCR (XCl , AIDER , AIDER , XCl , RMBAR ,1.0) 

CALL DPSMCR(XC2, AIDER, AIDER, XC2,RNBAR, 1.0) 

CALL DPSMCR(XC3 , AIDER, AIDER, XC3 ,RBBAR, 1.0) 

RETDRR 

END 

SUBROUTINE MEQMI (RATIO) 

This calculates the ratio (eq. magnetization/ini. magnetization^^^ 

COMPLEX DCl , DC2 , DC3 , DC4 , DC5 , DC6 , DC7 
COMPLEX RMBAR, RNBAR,RBBAR 
COMPLEX HD,DKLP,H1,H3,Q 

™ . R. THETA 

COMMON/BLK5/DC1 ( 1 6 , 1 6 ) , DC2 ( 1 6 . 1 6 ) , DC3 ( 16 , 1 6 ) , 

1 DC4fl6 16 ) , DC5 ( 16 , 16 ) , DC6 ( 16 , 16 ) , DC7 ( 16 , 16 ) 

COMMON/BLK7/HD{ 16,16), DKLP( 1 6 , 16 ) , HI ( J > m ft ^ i R > ^PNRARM 6 ^ 16 ) 
COMMON/BLK9/RBBAR( 16 , 16 ) ,RLBAR( 16,16) ,RMBAR( 16 , 16 ) , RNBARC 16,16) 

CALL H1H3 

CALL MTADCC(H1,H3,DC1,1.0) 

CALL MTMLCC(Q,DC1,DC2) 

TR1=(0.0,0.0) 

DO 100 1=1,16 
100 TE1=TR1+DC2(I,I) 

TR2=(0. 0,0.0) 

CALL MTMLCC (RBBAR , DCl , DC2 ) 

DO 110 1=1,16 
110 TR2=TR2+DC2(I,I) 

TR3=(0. 0,0.0) 
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All 


TRNR=(0. 0,0.0) 

TRDR=(0. 0,0.0) 

CALL MTMLCC(I>C1,DC1,DC2) 

CALL MTMLCC(Q,RBBAR,DC3) 

DO 120 1=1,16 
TR3=TR3+DC2(I,I) 

1 20 TRDR=TRDR+DC3 (1,1) 

TRDR=TRDR=^TR3 
TR3=(0. 0,0.0) 

CALL CMTRCC ( RNBAR , DKLP , DC2 ) 

CALL CMTRCC (DC2,RBBAR.DC2) 

CALL MTMLCC(DC2,H1,DC3) 

DO 130 1=1,16 
130 TR3=TR3+DC3(I,I) 

TR4-(0. 0,0.0) 

CALL CMTRCC(RNBAR,RBBAR,DC3) 

CALL CMTRCC (DCS, RNBAR, DC3) 

CALL MTMLCC(DC3,H1,DC1) 

DO 140 1=1,16 

TRNR='ITil*(TR2+(THETA*TA0*TR3/12.0) + (TEIETA*THETA*TR4/24.0) ) 

RATIO=TRNR/TRDR 

RETURN 

END 


SUBROUTINE PLOT(X,Y,N) 

This subroutines creates a line printer plot in the xiie 
corresponding to UNIT=6.NAG subroutine 601AGF is usedill. 
Execution comand should be .ex sqmgnt,/sea pub:nag.rel. 


DIMENSION X(N),Y(N),ISORT(100) 
CALL G01AGF(X , Y , N , ISORT ,110,67,1) 
RETURN 
END 


C 


C 


SUBROUTINE H1H3 

This calculates first and 


third order Floquet Hamiltonian. 


COMPLEX DCl , DC2 , DC3 , DC4 , DC5 , DC6 , DC7 , RMBAR , RNBAR , RBBAR 

a! !b?C , D . eta . TAO , DELWQ . THETAO , E . THETA , HL . KRKK 

COMMON/BLK5/DC1 ( 1 6 , 1 6 ) . DC2 ( 1 6 , 1 6 ) , DCS ( 1 6 , 1 6 ) , 

1 DC4('16 16),DC5(16,16),DC6(16,16),DC7(16,16) 

COMMON/BLK7/HD ( 1 6 , 16) ,DKLP(16, 16), HI (16, 16) ,H3(16, 16),Q(16,16) 
cSn^LK 9:^BAR( 16 , 16 ) , RLBAR( 16,16), RMBAR( 16,16) .RNBAR( 16,16) 
! ! n ! ! 5 COMPUTE HI J I ! ! ! 1 !• i 


CON1=THETA*0 . 5/TAO 

CALL MTADCR ( HD , RLBAR , DKLP , -DELWQ ) 

CALL MTADCC(DKLP, RNBAR, HI, CONI) 

I J I 1 t I I COMPUTE H3 !!!!!!!! 
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CALL CMTRCC ( RNBAR , DKLP , DC 1 ) 

CONI =THETA/1 2.0 

CON2=CONl*TAU 

CONl=CONl*THETA 

CALL CMTRCC(DC1 , RNBAR, DC2 ) 

CALL CMTRCC ( DCl , DKLP , DC3 ) 

DO 100 1=1,16 
DO 100 J=l,16 

100 H3(I, J)=C0N1*DC2(I, J)+CON2*DC3{I, J) 

RETURN 
END 


SUBROUTINE CMTRCC (A1 , B1 , C) 

This subroutine is written for evaluating the commutators of 
operators (complex Hermitian matrices.) 

COMPLEX Al(16,16) ,B1( 16,16 ) ,C( 16,16 ), Cl (16,16 ) ,C2( 16, 16) 

CALL MTMLCC(A1,B1,C1) 

CALL MTMLCC(B1,A1,C2) 

CALL MTADCC(C1,C2.C, -1.0) 

RETURN 

END 


SUBROUTINE DPSMCC(A1 , B1 , Cl , D1 ,D,CON) 

This and the following three subroutines are written for the 
routine evaluation of one particular kind of operation. 
involves first taking direct product of A1 and Bl.then of Cl 
and D1 and finally summing the two results. The routine names can 
be decoded as follows: DP— >direct product; SM-->sum;C >complex 
R — >real. 

COMPLEX Al(4,4) ,B1(4,4) ,C1(4,4) ,D1(4,4) ,DC1(16,16) ,DC2(16,16) , 

1 D(16,16) 

CALL DRPRCC ( A1 , B1 , DCl ) 

CALL DRPRCC ( Cl, D1,DC2) 

CALL MTADCC(DCl,DC2,D,CON) 

RETURN 

END 


SUBROUTINE DPSMRR(A1 ,B1 ,C1 ,Dl^D,CON) ^ 

DIMENSION Al(4,4) ,B1(4,4) ,C1(4,4) ,D1(4,4) ,D2(16, 16) ,D3(16, 16) , 

1 D(16,16) 

CALL DRPRRR(A1,B1,D2) 

CALL DRPRRR(C1,D1,D3) 

CALL MTADRR(D2,D3,D,CON) 

return 

END _ 

SUBROUTINE DPSMCR( A1 , B1 , Cl , D1 , D , CON ) 
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C 

C 


COMPLEX Al(4,4) ,D1(4,4) ,DC1(16, 16) , DC2( 16 , 16 ) ,D( 16, 16 ) 
DIMENSION B1(4,4),C1(4,4) 

CALL DRPRCR(A1,B1,DC1) 

CALL DRPRRC(C1,D1,DC2) 

CALL MTADCC ( DCl , DC2 , D , CON ) 

RETURN 

END 


SUBROUTINE DPSMRC ( A1 , B1 , Cl , D1 , D , CON ) 


COMPLEX Cl(4,4) ,D1(4,4) ,DC2(16,16) ,D(16,16) 
DIMENSION A1 ( 4 , 4 ) , B1 ( 4 , 4 ) , D2 { 16 , 16 ) 

CALL BPPRRR(A1,B1,D2) 

CALL DRPRCC{C1,D1,DC2) 

CALL MTADRC(D2,DC2,D,CON) 

RETURN 

END 


SUBROUTINE DRPRRR(A1,B1,C1) written 

This subroutine and the following many subroutines are written 
for matrirmanipulations-The routine names can be deeded as 
follows; MT->matrlx; AD-> addition; ML-> multiplication, 
DRPR->direct product; R->real; C->complex. 


DIMENSION Al(4,4) ,B1(4,4) ,C1(16, 16) 
IR=0 

DO 120 1=1,4 
IC=0 

DO 110 J=l,4 
DO 100 11=1,4 
DO 100 Jl=l,4 

100 CKIR+Il, IC+J1)=A1(I, J)*B1(I1, Jl) 
110 IC=4*J 
120 IR=4=t^I 
RETURN 
END 

SUBROUTINE DRPRRC(A1 ,B1,C1) ^ 

iii?oi^A?(4.4,,Bl(4.4,.Cl(16.16, 

IR=0 

DO 120 1=1,4 
IC=0 

DO 110 J=l,4 
DO 100 11=1,4 

DO 100 Jl=l,4 Ti ^ 

100 C1(IR+I1» IC+J1)=A1(I, J)*B1(I1, Jl) 

110 IC=4*J 



120 


IR=4*I 

RETURN 

END 


SUBROUTINE DRPRCC ( A1 , B1 , Cl ) 


COMPLEX Al(4,4) ,B1(4,4) ,C1(16,16) 
IR=0 

DO 120 1=1,4 
IC=0 

DO 110 J=l,4 
DO 100 11=1,4 
DO 100 Jl=l,4 

100 C1{IR+I1,IC+J1)=A1 ri, J)*B1(I1,J1) 

110 IC=4*J 
120 IR=4*I 
RETURN 
END 

SUBROUTINE DRPRCR(A1.B1,C1) 

COMPLEX Al, Cl 

DIMENSION Al(4,4) ,B1(4,4) ,C1(16,16) 
IR=0 

DO 120 1=1,4 
IC=0 

DO 110 J=l,4 
DO 100 11=1,4 
DO 100 Jl=l,4 

100 C1(IR+I1. IC+J1)=A1(I, J)*B1(I1, Jl) 

110 IC=4*J 
120 IR=4*I 


RETURN 

END 


c 

SUBROUTINE MTADCC(A1,B1,C1,C0N) 


c 

100 

COMPLEX A1(16,16) ,B1C16,16) ,C1(16,16) 

DO 100 1=1,16 

DO 100 J=l,16 

Cl ( I , J ) =A1 ( I , J ) +©1 ( I , J ) *CON 

RETURN 

END ___ _ 


c 

SUBROUTINE MTADRE(A1,B1,C1,C0N) 


c 

100 

DIMENSION A1(16,16),B1(16,16) ,C1(16,16) 

DO 100 1=1,16 

DO 100 J=l,16 

Cl ( I , J ) =A1 ( I • J ) •‘■Bl ( I » 'J ) 
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C 
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RETURN 

END 


SUBROUTINE MTADRC( Al , B1 , Cl , CON ) 


COMPLEX Bl.Cl 

DIMENSION A1(16,16),B1(16,16).C1(16,16) 

DO 100 1=1,16 

DO 100 J=l,16 

CKI,J)=A1(I,J)+B1(I,J)*C0N 

RETURN 

END 



SUBROUTINE MTADCR(Al,Bl,Ci,CON) 


COMPLEX Al, Cl ^ 

DIMENSION Al ( 16 , 16 ) , B1 ( 16 , 16 ) , Cl ( 16 , 16 ) 
DO 100 1=1,16 
DO 100 J=l,16 

100 C1(I.J)=A1(I,J)+B1(I.J)*C0N 

RETURN 
END 

SUBROUTINE MTMLRRCAl ,B1 .Cl ) 

DIMENSION Al ( 16 , 16 ) , B1 ( 16 , 16 ) , Cl ( 16 , 16 ) 
DO 100 1=1,16 
DO 100 J=l,16 
CKI, J)=0.0 
DO 100 K=l,16 

C1(I, J)=C1(I, J)+A1(I,K)*B1(K, J) 

100 CONTINUE 
RETURN 

END 

SUBROUTINE MTMLCC(A1 ,B1 ,C1 ) 

COMPLEX Al(16,16) ,B1(16,16),C1(16,16) 
DO 100 1=1,16 
DO 100 J=l,16 
CKI, J) = (0. 0,0.0) 

C1(I?J)=C1(I.J)+A1(I,K)*B1{K,J) 

100 CONTINUE 
RETURN 
END 



